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ABSTRACT 


This  study  deals  with  the  low  mass-density  ratio  hydroelastic 
behavior  of  fully  wetted  lifting  surfaces  of  sufficiently  high  aspect  ratio 
to  justify  the  fundamental  assumptions  of  two-dimensional,  unsteady  hy¬ 
drodynamic  fluid  forces  and  the  concept  of  a  rigid  streamwise  section. 
Hence,  this  study  is  particularly  applicable  to  lifting  surfaces  in  marine 
usage.  In  particular,  the  strut  supports  of  hydrofoil  boats  which  operate 
at  sub -cavita ting  speeds  may  be  regarded  as  typical  of  the  hydroelastic 
systems  toward  which  this  investigation  is  directed. 

Part  I  of  this  investigation  is  specifically  concerned  with  the  low 
mass-density  ratio  flutter  behavior  of  the  two -degree -of -freedom  "typical 
section."  Part  I  shows  that  there  exists  a  "critical"  mass-density  ratio 

(U  )  below  which  the  classic  "typical  section"  analysis  will  predict  no 
a 

dynamic  instabilities.  However,  some  extremely  sensitive  cause  and 

effect  relationships  are  noted  and  discussed  in  Part  I. 

Part  II  of  this  investigation  deals  with  the  low  mass -density  ratio 

hydroelastic  behavior  of  fully  wetted,  uniform,  cantilevered  struts  or 

other  lifting  surfaces,  both  with  and  without  an  angle  of  sweepback.  As 

a  result  of  this  study,  it  has  been  ascertained  that  the  low  mass -density 

ratio  (u  <  u  )  flutter  behavior  of  such  structures  can  only  be  properly 
a 

predicted  if  the  spanwise  character  of  the  strut  is  explicitly  included  in 
the  mathematical  model  used  in  the  analysis.  Further,  it  is  shown  in 
Part  II  that  a  positive  non -zero  sweep  angle  is  essential  for  low  mass- 
density  ratio  bending -tor s ion  flutter  of  uniform  cantilevered  struts  to 
occur.  In  addition,  Part  II  shows  what  structural  parametric  combina¬ 
tions  may  be  expected  to  have  a  low  mass -density  ratio  flutter  problem. 
Finally,  a  physical  explanation  of  the  phenomenon  is  offered. 

Viewed  as  a  whole,  this  study  shows  that  the  low  mass-density 
ratio  region  is  one  of  extreme  sensitivity  where  small  parametric  varia¬ 
tions  can  result  in  large  changes  in  the  hydroelastic  behavior  pattern  of 
a  lifting  surface. 
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GENERAL  INTRODUCTION 

Recent  interest  in  hydrofoils  has  quite  naturally  resulted  in 
numerous  investigations  of  the  engineering  problems  associated  with  the 
design  of  such  craft.  Because  of  the  obvious  similarities  between  hydro¬ 
foils  and  airfoils,  it  had  been  assumed  that  many  of  the  problems  encoun¬ 
tered  in  aircraft  design  could  well  be  expected  to  arise  with  hydrofoils. 
However,  several  investigations  conducted  during  the  last  decade  seemed 
to  indicate  that  the  problem  list  might  well  have  one  important  exception. 
In  particular,  these  investigations  (Refs.  2,  3, 4,  5)  indicated  that  flutter, 
a  problem  of  considerable  magnitude  in  aircraft  design,  might  not  occur 
with  hydrofoils.  The  fact  that  the  lifting  surfaces  of  some  marine  craft 
had  long  been  subjected  to  dynamic  pressures  higher  than  those  attained 
by  the  fastest  aircraft  and  had  not  encountered  flutter  strengthened  the 
notion. 

Specifically,  the  investigators  referenced  above  demonstrated,  by 
simplified  analytical  arguments  based  upon  a  simplified  representation  of 
the  problem,  that  when  the  ratio  of  the  mass  of  a  foil  section  of  unit  span 
to  the  mass  of  a  cylinder  of  surrounding  fluid  (^j)  fell  below  some  "critical" 
value,  classical  flutter  was  unlikely.  Further,  a  limited  amount  of  ex¬ 
perimental  evidence  seemed  to  confirm  what  the  theoretical  analyses  in¬ 
dicated.  Accordingly,  by  I960,  the  majority  of  the  literature  available 
on  the  subject  of  low  mass -density  ratio  flutter  concluded  that  the  occur¬ 
rence  of  the  phenomena  was  highly  unlikely.  However,  in  spite  of  this 
pleasant  outlook,  several  reports  purported  to  document  actual  observed 
cases  of  hydrofoil  or  other  very  low  mass-density  ratio  flutter,  (e.g.  , 

Ref.  7).  These  reports  were  not  numerous  and  the  flutter  occurrences 
usually  not  well  documented.  Nevertheless,  the  possible  consequences 
of  ignoring  these  reports  necessitated  a  thorough  investigation  of  the 
subject.  To  complicate  matters  even  further,  Herr  (5)  had  obtained 
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unstable  flutter  solutions  for  certain  discrete  swept  foils  in  spite  of  the 
fact  that  his  own  simplified,  "typical  section"  theory  had  indicated  that 
such  solutions  should  be  non-existent. 

Obviously,  the  contradictory  evidence  outlined  above  necessitated 
that  other  more  exhaustive  investigations  be  undertaken.  One  of  these 
investigations  was  the  one  undertaken  by  Baird,  Squires,  and  Caporali, 
(Ref.  10).  In  the  course  of  that  investigation,  it  was  demonstrated  both 
analytically  and  experimentally  that  coupled  bending -tor sion  flutter  could, 
in  general,  occur  regardless  of  the  mass-density  ratio.  There  were, 
however,  several  disquieting  aspects,  particularly  with  regard  to  the 
theoretical  portion  of  the  study. 

First,  the  analytical  method  which  successfully  predicted  the  ex¬ 
istence  of  a  flutter  boundary  was  the  standard  "production"  modal  analysis 
technique  in  which  the  flutter  mode  was  assumed  to  be  a  linear  combina¬ 
tion  of  a  finite  number  of  the  natural  modes  of  the  structure.  The  disturb¬ 
ing  aspect  of  these  analyses  was  that  a  preliminary  check  based  upon  the 
method  outlined  in  Ref.  (5)  (derived  from  the  "typical  section"  analysis) 
showed  the  mass -density  ratio  below  which  no  flutter  should  occur  to  be 
well  above  the  mass -density  ratio  associated  with  the  test  struts.  Quite 
naturally,  the  first  thought  was  that  the  section  analysis  was  an  over¬ 
simplification  and  that  the  modal  analysis  succeeded  because  the  character 
of  the  structure  was  more  accurately  represented  by  the  inclusion  of  more 
modes.  This  theory,  however,  was  shown  to  be  less  than  entirely  correct 
when  it  was  demonstrated  that  flutter  could  be  predicted  by  the  modal 
technique  when  as  few  as  two  modes  were  used.  Admittedly,  the  predic¬ 
tions  were  not  accurate.  Nevertheless,  the  discrepancies  were  obviously 
a  matter  of  degree  and  not  the  result  of  fundamental  errors  in  principle. 
Thus,  it  became  apparent  that  the  difference  between  the  two  methods  of 
analysis  was  much  more  subtle.  This  was  indeed  unfortunate  as  the  modal 
technique,  although  well  suited  to  production  type  work,  involved  a  degree 
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of  numerical  complexity  that  made  a  true  physical  understanding  of  the 
phenomena  difficult  to  obtain.  Inasmuch  as  the  possibility  of  classical 
flutter  had  been  established,  the  extrapolation  of  the  results  for  use  as 
a  design”tool”  required  precisely  this  physical  understanding.  Histori¬ 
cally,  this  type  of  understanding  had  always  been  obtained  from  the  more 
simple  but  well  proven  two -degree -of -freedom  ’’typical  section”  analysis, 
such  as  outlined  in  Ref.  (1).  Since  conclusions  based  upon  results  ob¬ 
tained  from  the  typical  section  analysis  had  been  proven  wrong,  two  pos¬ 
sibilities  existed.  Either  the  erroneous  conclusions  had  been  the  result 
of  mis -interpretations  of  the  simplified  analysis,  or  the  fundamental 
analytical  technique  had  been  inadequate.  The  study  of  the  two  possibili¬ 
ties  and  the  conclusions  resulting  therefrom  is  precisely  the  subject 
matter  contained  herein. 

It  would  seem  that  the  possibility  of  a  mis -inter pretation  should 
have  been  easy  to  discard.  However,  preliminary  work  by  Brunelle 
(Refs.  6  and  8)  and  Squires  (Ref.  12)  showed  the  possibility  definitely 
existed  that  all  of  the  pertinent  information  had  not  been  extracted  from 
the  simplified  analysis.  Essentially,  this  preliminary  work  implied  that 
it  might  be  possible  to  show  that  the  ’’critical”  mass -density  ratio  did  not 
in  fact  establish  a  true  limit. 

The  second  possibility  was  also  quickly  confirmed.  The  author  and 
others  (Ref.  10)  noted  that  the  unstable  flutter  solutions  obtained  by  the 
modal  technique  were  always  obtained  for  swept,  cantilevered  lifting  sur¬ 
faces.  Thus,  the  angle  of  sweep,  and  possibly  the  cantilevered  mounting, 
could  be  singled  out  as  the  possible  differences  between  the  two  methods 
inasmuch  as  the  analysis  based  upon  the  two -degree -of -freedom  section 
had  not  explicitly  included  sweep  or  the  spanwise  characteristics  of  the 
cantilevered  structures  as  distinct  parameters. 

Thus,  this  investigation  was  undertaken  in  an  effort  to  explain  low 
mass -density  ratio  classical  flutter.  This  investigation  can  be  described 
as  consisting  of  two  related  parts. 
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1)  A  comprehensive  investigation  of  the  two  “degree  -of  “freedom 
section  analysis  using  both  the  simplified  quasi-steady  hydro¬ 
dynamics  and  the  full  unsteady  hydrodynamics;  the  purpose  of 
this  phase  being  to  determine  whether  finite -valued  flutter 
boundaries  could  exist  at  very  low  mass-density  ratios  for 
unswept,  fully  wetted  sections,  (i.e.  ,  values  of  p  below  the 
'’critical11  p  of  a  given  section.  ) 

2)  A  comprehensive  investigation  of  the  characteristics  peculiar 
to  a  swept,  cantilevered  lifting  surface  to  determine  the  effects 
of  such  characteristics  upon  the  low  mass -density  ratio  flutter 
boundaries  of  such  structures. 

Principally,  this  documentation  of  the  investigation  outlined  above 
shows  that: 

1.  For  any  given  structure  assumed  to  be  represented  by  a  "typical 
section"  there  exists  a  mass-density  ratio  below  which,  for  all 
practical  purposes,  no  dynamic  hydroelastic  instabilities  can 

be  predicted,  (i.  e.  ,  for  u  <  p  .  ) 

3. 

2,  The  low  mass -density  ratio  flutter  behavior  of  a  swept,  canti¬ 
levered  lifting  surface  is  directly  a  function  of  the  structural 
characteristics  peculiar  to  such  a  structure;  namely  the  span- 
wise  variation  in  the  local  pitch  to  bending  ratio  and  the  sweep 
of  the  elastic  axis. 

In  addition,  related  items  of  interest  uncovered  in  the  course  of 
determining  the  above  are  included  herein. 
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DISCUSSION  AND  RESULTS 

PART  I  -  THE  TWO-DIMENSIONAL,  TWO -DEGREE -OF -FREEDOM 
TYPICAL  SECTION 

I -A  -  Background 

Flutter  analysts  in  the  aircraft  field  have  long  used  the  technique 
of  examining  the  behavior  of  a  representative  two -degree -of -freedom,  two- 
dimensional  section  of  a  structure  to  qualitatively  determine  the  flutter 
characteristics  of  the  entire  structure.  Since  even  this  simplified  analysis 
leads  one  to  two  transcendental  equations  which  must  be  simultaneously 
satisfied  when  the  complete  unsteady  hydrodynamic  theory  is  used,  the 
usual  practice,  at  least  when  only  qualitative  results  have  been  desired, 
has  been  to  use  quasi-steady  hydrodynamic  coefficients  obtained  by  setting 
the  ordinarily  complex  Theodorsen  function,  C(k),  equal  to  either  1  or 
(See  Ref.  (1)  or  Ref.  (13)  for  a  complete  explanation  of  the  Theodorsen 
function.  )  Briefly,  these  two  values  of  C(k)  are  the  limiting  values  ob¬ 
tained  by  letting  the  reduced  frequency  k  “  0  and  k  -4  “respectively.  Ap¬ 
pendix  A  shows  a  plot  of  the  dependence  of  the  Theodorsen  function  upon  k. 
(In  order  to  differentiate  between  the  coefficients,  the  set  obtained  with 
C(k)  =  ^  will  be  referred  to  as  the  qua  si -unsteady  set). 

Essentially,  the  conclusion  that  hydrofoil  or  other  low  mass- 
density  ratio  flutter  was  unlikely  to  occur  was  largely  the  result  of  ex¬ 
tending  the  simplified  techniques  discussed  above  to  the  low  mass -density 
ratio  regime.  Several  investigators  (e.g.  ,  Ref.  5)  demonstrated  that  the 
flutter  boundary  of  some  physically  plausible  two -degree -of -freedom, 
two-dimensional  sections  could  be  expressed,  as  a  function  of  the  mass 
density  ratio,  as  shown  in  Fig.  1. 

The  normal  procedure  followed  was  to  show  that  under  the  quasi- 

uf 

steady  assumption  the  flutter  speed  coefficient  (— —  )  could  be  expressed 

bu 

a 

as  the  ratio  of  two  polynomials  in  u.  By  setting  the  denominator  of  the 
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polynomial  equal  to  zero,  it  was  possible  to  demonstrate  that  for  some 
’’critical"  value  of  the  mass -density  ratio,  the  flutter  speed  would  be 
infinite.  This  value  has  been  termed  u  -  asymptote  (u  ).  The  equation 

cl 

from  which  has  been  determined  is  reproduced  here  to  reveal  its 
general  character.  This  equation  is  discussed  in  considerable  detail 
elsewhere  in  this  report. 


(jJ  2  2 

1  h  1  a8  a  1 

2+(u  '  r  3  (  2  ”  8  +  1 6 

a  a 


Ma  = 


(l-l) 


id  3 

T  +  a  +  x  -  (  )  — —  [x  (a2  +  a  -  j)  +  r  ®  ( y  +  a )] 

Z  cvcj  p  a  4  a  £ 

a  r 
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It  was,  of  course,  well  known  that  finite  flutter  speeds  existed  for  mass- 
density  ratios  larger  than  u  •  However,  the  notion  that  no  finite  flutter 

cl 

boundaries  existed  for  mass -density  ratios  less  than  the  ratio  which 
yielded  an  infinite  flutter  speed  was  extrapolated  from  the  fact  that  no 
boundaries  were  ever  found  for  the  particular  combinations  of  system 
parameters  used  in  the  calculations.  Thus,  the  rather  widely  accepted 
idea  of  a  minimum  mass -density  ratio  below  which  no  unstable  solutions 
existed  was  really  an  assumption  taken  to  be  valid  until  (and  if)  some 
negating  calculations  could  be  shown.  To  express  the  issue  somewhat 
more  precisely,  it  had  not  been  established,  even  for  the  simplified 
quasi-steady  and  qua  s  i  -unsteady  cases,  that  no  physically  plausible  com¬ 
bination  of  parameters  existed  which  could  result  in  a  re -appearance  of 
real  solutions  for  some  mass -density  ratio  below  the  so-called  "critical11 
value.  Mathematically  stated,  it  had  not  been  demonstrated  that  the 
ratio  of  polynomials  in  u  could  not  undergo  a  second  sign  reversal  for 
some  positive  mass -density  ratio  less  than  \j  .  Therefore,  a  careful 

cl 

re-examination  of  the  simplified  analyses  described  above  became  the 
starting  point  of  investigations  attempting  to  determine  why  the  two- 
degree -of -freedom  section  analysis  had  apparently  led  to  erroneous 
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conclusions.  The  results  of  several  of  these  re-examinations  have  been 
documented  in  Refs.  6,  8,  and  12.  In  addition,  the  author  has  also  studied 
the  problem  and  certain  aspects  of  the  results  are  discussed  herein.  To 
summarize,  these  investigations  indicated  that,  under  specia  1  c ircum  - 
stances,  certain  unswept,  fully  wetted  lifting  surfaces  apparently  had 
finite  valued  flutter  boundaries  for  mass -density  ratios  well  below  the 
’’critical'1  value.  However,  the  results  were  not  conclusive. 

Using  the  quasi-steady  hydrodynamic  theory  associated  with  a  re¬ 
duced  frequency  (k)  of  zero,  it  was  established  that  at  least  some  physi¬ 
cally  plausible  configurations  had  finite  valued  flutter  boundaries  at  very 
low  mass-density  ratios.  For  the  special  case  of  an  "in  vacuum"  un¬ 
coupled  plunging  to  uncoupled  pitching  frequency  ratio  of  zero,  Squires 
(Ref.  12)  established  that  the  existence  of  solutions  for  u  <  u  was  entirely 

a  function  of  the  center  of  gravity  location  (x  ).  Brunelle  also  included 

m 

several  examples  in  Refs.  6  and  8.  However,  Brunelle  also  noted  that 

the  reduced  frequencies  calculated  for  points  on  the  flutter  boundaries 

(u  <  U  )  which  were  obtained  from  the  quasi-steady  hydrodynamics  tended 
a 

to  invalidate  the  fundamental  quasi -steady  assumption  of  a  low  reduced 
frequency.  Attempts  to  obtain  finite -valued  flutter  boundaries  for  u  <  U 
using  the  qua  si -unsteady  hydrodynamic  theory  associated  with  a  reduced 
frequency  of  infinity,  for  the  very  same  configurations,  were  completely 
unsuccessful.  Thus,  the  unfortunate  and  unexpected  strong  dependence  of 
the  low  mass -density  ratio  solutions  on  the  reduced  frequency  limit  as¬ 
sumed  was  discovered.  The  term  unexpected  is  used  because  at  the  mass- 
density  ratios  normally  encountered  in  aircraft  flutter  analyses,  the  flutter 
boundary  is  normally  only  weakly  dependent  upon  the  reduced  frequency. 
The  discovery  is  termed  unfortunate  because  it  implied  that  the  determina¬ 
tion  of  the  validity  of  the  two -degree -of -freedom  section  analysis  for  low 
mass -density  ratios  would  require  the  inclusion  of  the  complete  complex 
value  of  the  Theodor  sen  function,  C(k),  in  the  analysis.  This  step  would 
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have  increased  the  complexity  of  the  problem  by  at  least  an  order  of  magni¬ 
tude.  The  transcendental  nature  of  the  resulting  equations  would  have  re- 
quired  that  a  lengthy  numerical  technique  be  used  to  obtain  solutions  with 
the  result  that  general  statements  regarding  classes  of  configurations 
would  have  been  virtually  impossible  to  make  without  a  prohibitive  amount 
of  computational  effort.  Accordingly,  a  systematic  method  for  evaluating 
the  two -deg  ree -of -freedom  section  analysis  was  devised.  The  method  and 
the  conclusions  resulting  from  its  application  are  the  subject  treated  in 
Part  I.  In  effect,  Part  I  of  this  report  might  be  regarded  as  "an  analysis 
of  an  analysis .  M 

I-B  -  Preliminary  Evaluation 

The  difficulties  and  inconsistencies  encountered  when  attempting 
to  study  numerically  the  flutter  behavior  of  a  two -degree -of -freedom  typical 
section  while  incorporating  only  the  approximate  quasi-steady  or  quasi - 
unsteady  hydrodynamic  theories  have  been  discussed.  Clearly,  these  in¬ 
consistencies  imply  that  the  full  two-dimensional,  incompressible,  un¬ 
steady  theory  must  be  employed  to  assure  the  validity  of  any  conclusions. 

The  development  of  the  complex  "flutter  determinant"  by  using  the 
"exact"  complex  Theodorsen  function,  C(k),  is  a  classical  derivation.  Vir¬ 
tually  any  text  on  the  subject  of  flutter  (e.g.  ,  Refs.  13,  14,  15)  will  contain 
the  development  in  some  form.  For  completeness,  an  abbreviated  develop¬ 
ment  is  included  herein  as  Appendix  A.  However,  for  purposes  of  this  dis¬ 
cussion,  it  is  convenient  to  show  several  of  the  steps  here. 

By  using  Lagrange's  equation,  the  equations  of  motion  for  the 
typical  section  of  Figure  Z  may  be  written  as: 

mh  +  S  a  +  mu  2  h  =  -  L 
a  h 

S  h  +  I  a  +  I  w  2a  =  M 
a  a  a  a 


(1-2) 
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where  L  is  the  lift  per  unit  span  and  M  the  moment  about  the  elastic  axis 
per  unit  span.  The  assumption  of  simple  harmonic  motion  at  the  inception 
of  the  flutter  condition  changes  equations  (1-2)  into: 


-cj  2  mh  -  u  2  S  a  +  u  2mh  =  -L. 
f  f  a  h  ha  r 


-cj  2  S  h  -  cj  2  I  a  +  cj  2  I  cl 
fa  f  a  a  cl 


M, 


har 


(1-3) 


The  terms  L  and  M  are  dependent  upon  the  Mach-number  range 
considered.  Since  this  study  is  fundamentally  concerned  with  the  behavior 
of  foils  in  water,  the  expressions  used  to  replace  L  and  M  are  obtained 
from  incompressible  theory.  This  allows  L  and  M  to  be  put  into  the  nota¬ 
tion  of  Ref.  16.  The  resultant  homogeneous  equations  constitute  an  alge¬ 
braic  eigenvalue  problem  with  meaningful  solutions  occurring  at  those  com¬ 
binations  of  speed  and  frequency  for  which  the  characteristic  determinant 
vanishes.  This  is  called  the  flutter  determinant  and  may  be  written  as: 


CJ  2  CJ  2 

CuCl  tux  +  [  L  -L  (i  +a)]} 

u  2  u*  h  °  a  h 

a  f 


CJ 


(uxa  +Mh-Lh(l  +a)l  turo3[l  +Mo 


-(L  +M  )<|+a)  +  L  (£+a)3} 
ct  h  n 


=  0 


(1-4) 
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where  the  L  ,  L  ,  and  M  are  functions  of  the  reduced  frequency  (k)  only 
ha  a  J  7 

and  M  is  equal  to  §  for  the  incompressible  case,  (see  Appendix  A).  Ex- 
h 

animation  of  the  determinant  reveals  that  the  flutter  frequency  ratio 
and  the  reduced  frequency  (k)  are  dependent  upon  the  values  of  five  dimen¬ 
sionless  parameters.  These  parameters,  which  are  listed  in  the  table  of 
symbols  are: 


a  -  the  elastic  axis  location 

(J  / (J  ^  -  the  plunging -pitching  frequency  ratio 

x  -  the  nondimensional  static  unbalance 
a 

r^  -  the  nondimensional  radius  of  gyration 
U  -  the  mass-density  ratio. 


It  is  thus  apparent  that  a  reasonably  comprehensive  parameter 
study  is  not  a  trivial  undertaking.  This  fact  is  made  even  more  apparent 
when  one  briefly  considers  the  actual  details  of  solving  Eq.  (1-4).  The 
complex  form  of  the  aerodynamic  coefficients  causes  the  determinant  to 
expand  into  real  and  imaginary  equations.  Because  of  the  nature  of  the 
hydrodynamic  terms,  both  equations  include  complicated  functions  of  k. 
Thus  to  solve  for  the  flutter  frequency  ratio  and  the  reduced  frequency 
with  the  above  physical  properties  known  requires  a  lengthy  trial  and 
error  process.  Hence,  it  may  be  concluded  that  a  comprehensive  study 
of  the  behavior  of  the  two -degree -of-freedom  analysis  as  a  function  of 
the  above  noted  system  parameters  presents  practical  computational  and 
data  processing  difficulties  of  such  a  magnitude  that  the  time  required  for 
completion  would  be  prohibitive;  even  with  the  use  of  high  speed  digital 
equipment  (IBM  7090)  and  the  range  of  interest  for  each  individual  para¬ 
meter  confined  by  practical  considerations.  Furthermore,  past  exper¬ 
ience  (Refs.  5,  10)  has  indicated  that  the  low  mass -density  ratio  regime 
may  be  characterized  by  an  extreme  sensitivity  to  small  changes  in  the 
values  of  the  system  parameters  -  thus  implying  that  unusually  small 
step  changes  might  be  required  in  order  to  correctly  ascertain  trends 
over  all  ranges  of  interest. 
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In  view  of  the  anticipated  computational  difficulties,  it  is  apparent 
that  some  scheme  which  serves  to  reduce  the  dimensional  size  of  the  prob¬ 
lem  is  highly  desirable  -  at  least  for  the  start  of  the  investigation.  A 
scheme  outlined  by  Brunelle  in  Ref.  6  is  promising  as  it  offers  the  pos¬ 
sibility  of  ignoring  two  of  the  system  parameters.  Fundamentally,  the 
analysis  involves  a  study  of  the  special  case  of  a  zero  mass -density  ratio. 

I-C  -  The  Zero  Mass -Density  Ratio  Analysis 

Solutions  for  the  Limiting  Zero  Mass -Density  Ratio  Case 

An  analysis  based  upon  a  zero  mass -density  ratio  has  several 
attractive  features.  One  is  concerned  with  the  fact  that  two  of  the  system 
parameters  do  not  have  to  be  considered.  In  particular,  the  non-dimensional 
mass  unbalance  and  the  non-dimensional  radius  of  gyration  do  not  appear  in¬ 
dependently  in  the  flutter  equations.  Physically  this  could  be  interpreted  as 
the  result  of  the  inertial  effects  of  the  structure  becoming  negligible  when 
compared  to  the  virtual  or  apparent  inertial  effects  of  the  surrounding  fluid. 

A  second  attractive  feature  is  that  the  zero  mass -density  ratio  case  is  a 
logical  limiting  case  for  the  low  mass-density  ratio  regime.  Thus,  one 
might  intuitively  expect  that  any  peculiarities  resulting  from  a  low  mass- 
density  ratio  would  be  most  pronounced  for  the  limiting  condition.  Un¬ 
fortunately,  a  study  of  the  zero  mass -density  ratio  case  cannot  be  under¬ 
taken  without  introducing  a  different  set  of  non-dimensional  system  para¬ 
meters. 

It  may  be  seen  from  equations  (1-2)  and  (1-3)  that  although  the  equa¬ 
tions  describe  a  two -degree -of -freedom  spring-mass  system,  the  spring 

constants  k  and  k  do  not  appear  explicitly.  Rather,  the  spring  constants 
h  a 

have  been  replaced  by  functions  involving  the  uncoupled  natural  frequencies 

in  plunging  and  in  pitching  through  use  of  the  relations  k^  =  mu^3  and 

k  =  I  u  3.  As  was  pointed  out  by  Brunelle  (Ref.  6),  regarding  the  spring 
a  a  a 

constants  as  being  fixed  requires  that  during  the  limiting  process  of 
m  —  y  -♦  0,  u  3  and  w  a  both  approach  infinity  in  such  a  fashion  as  to  keep 
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and  at  their  predetermined  values.  It  was  further  noted  by  Brunelle 

that  as  u^2  approaches  infinity,  the  flutter  speed  coefficient  (U^/bw^)  has 

to  approach  zero  even  though  the  flutter  speed  might  in  reality  have  a  finite 

non-zero  value.  Hence,  it  is  apparent  that  the  substitutions  for  k.,  and  k 

ha 

noted  above  are  poorly  suited  for  investigating  the  case  of  the  mass,  and 
subsequently  the  mass -density  ratio,  approaching  zero,  since  all  flutter 
points  collapse  onto  the  origin.  It  is  then  obvious  that  the  study  of  the  zero 
mass -density  case  requires  that  k^  and  k^  be  expressed  in  terms  of  gene¬ 
ralized  masses  and  frequencies  that  do  not  become  zero  or  infinite  as  the 
mass -density  ratio  approaches  zero.  Following  Ref.  (6),  it  is  possible  to 
obtain  a  more  useful  set  of  generalized  masses  and  frequencies  by  writing 
the  equations  for  a  "typical  section"  (Fig.  2)  which  is  submerged  in  a 
stationary  fluid.  If  all  the  coupling  terms  are  suppressed,  the  equations 
may  be  written: 


[  m  +  tt  pb2  ]  h  +  k^h  =  0 


[  I  +  TT  pb4  (  77  +  a 2  )  ]  a  +  k  a  =0 
a  8  a 


Thus,  the  artificially  uncoupled  "in  fluid"  frequencies  become: 


m  +  pbs 


and  CJ  *2  = 
a 


la  +7rPb4<t  +  a2) 


If  the  zero  mass -density  ratio  is  taken  to  be  the  limit  reached  as 
the  mass  is  made  to  approach  zero,  rather  than  the  limit  achieved  as  the 
fluid  density  is  made  to  approach  infinity,  the  artificially  uncoupled  water 
frequencies  are  seen  to  approach  u*a  =  k^/frpb2  and  u  J3  =  k^/V  pb4  (  + as  ) . 

Clearly  this  implies  that  the  generalized  mass  and  the  generalized  moment 
of  inertia  -  for  the  zero  mass -density  ratio  case  -  are  determined  by  the 
density  of  the  fluid. 
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As  one  would  expect,  the  "in  fluid"  and  the  "in  vacuum"  frequencies 
are  directly  related,  the  exact  relationships  being  a  function  of  the  mass- 
density  ratio.  For  purposes  of  comparison  the  following  equations  are 
given: 


If  the  above  relations  are  used  in  conjunction  with  the  unsteady, 
two  dimensional,  incompressible  hydrodynamic  theory  exemplified  by  an 
arbitrary  C(k),  the  equations  of  harmonic  motion  of  the  typical  section  are 
given  by: 

(  1  +u)£  +  2C(k)l  +  (u  +  1  )(bu*/u)3II3  4  +  ( uxQ  -a  )a 


+  2[|  +  C(k)(£  -  a)>  +  2C(k)a  =  0 
(Ux^  -  a )l  -  2(a+|)C(k)£  +  [|jrQ3  +  (|+  a3)> 


(1-8) 


+  <*-a)[l  -2C(k)(|+a)]a  +  {far*  +  (  i-  +  aa  )](bu  */  U)a 
-  2(  a  +  i)C(k)}ft  =  0 
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For  purposes  of  suitably  non-dimensionalizing  the  above,  the  dif¬ 
ferentiations  indicated  are  with  respect  to  a  non-dimensional  time.  That 
da  b  3  ds 

is  (  )  indicates  -  =  (— )  -  ,  etc.  ,  as  is  indicated  in  the  list  of 

i  2  U  .2 

ds3  dt 

symbols . 

At  this  point,  Brunelle  solved  the  equations  of  motion  for  the  limit¬ 
ing  case  of  a  zero  mass -density  ratio  using  both  the  approximate  hydro¬ 
dynamics  associated  with  a  reduced  frequency  of  zero  and  of  infinity.  This 
solution  was  based  upon  the  customary  harmonic  substitution  dn/dsn  = 

(jo o  )ne^GJ  S  and  further  noting  that  u  =  (oj^/oj  *)/(uf/b  OJ  J).  The  results  are 
presented  below: 

For(j  =  0  and  k  =  0 


bo; 

a 


! q2  +  (1  +  8aa  ) [ II  ,a  ( 1  -  q)  -  q] 
2 /T  V  q  +  (2a  +  1)  n,s 


( I —9a ) 


GJ 

— =  /q  where  q  =  4a(Za  -  1)11,, 2  +  (1  +  8aa  )  (l-9b) 

OJ 

a 


For  =  0  and  k  =  oo 


bcj 

a 


+  8aa)[F;a(l  -q>)  -  q] 
q  +  (2a  +  1  )IT, a 


(I-10a) 


Id 


— —  =  / q  where  q  =  —  [  2(2a  -  1  )a  FI,3  +  ( 1  +  8aa  )] 

3  * 


c j 


(I-10b) 


a 


-  16  - 


These  results  are  seen  to  be  quite  simple  and  allow  for  a  direct 
calculation  of  the  frequency  and  flutter  speed  coefficient  as  a  function  of 
n,a  and  (a).  Unfortunately,  the  results  obtained  by  using  the  above  equa¬ 
tions  suffered  from  the  inconsistencies  described  in  the  previous  section. 
Namely,  when  equations  (1-9)  were  used,  many  practical  combinations  of 
n,3  and  (a)  were  found  which  yielded  real,  positive  values  for  U^/bcj^*  and 
However,  the  resultant  reduced  frequency  obtained  by  combining 
the  frequency  and  flutter  speed  coefficients  was  invariably  much  too  large 
to  justify  the  choice  of  the  quasi-steady  hydrodynamic  theory.  Rather, 
the  values  of  reduced  frequency  which  were  obtained  indicated  that  equa¬ 
tions  (1-10)  should  be  used.  However,  the  use  of  equations  (1-10)  failed 

to  yield  any  real,  positive  values  for  U  /bio  *  and  u  ,/cj  *  for  any  of  the 
3  f  a  {’  a  7 

n,3  and  (a)  combinations  checked.  Thus  the  use  of  the  "exact"  theory  is 

implied. 

To  obtain  the  equations  associated  with  the  "exact"  unsteady  hydro - 
dynamic  theory  and  a  mass -density  ratio  of  zero  requires  that  the  equa¬ 
tions  of  motion,  Eqs.  (1-8),  be  solved  for  the  zero -mas  s -density  ratio 
case  with  the  complex  form  of  the  Theodorsen  function,  (C(k)  =  F  +  jg), 
employed.  Again  the  solution  is  based  upon  the  customary  harmonic 
substitution  and  follows  closely  the  development  shown  in  Appendix  A. 

By  making  the  appropriate  substitutions,  and  after  a  considerable 
amount  of  algebraic  manipulation,  one  form  of  the  resultant  flutter 
equations  may  be  written  as  shown  on  the  following  page  as  Eq.  ( I  - 1  1 ) . 

The  symbols  used  in  Eq.  (I  —  1 1 )  are  all  included  in  the  list  of 


symbols , 
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U  s  n»(l  +  8as  ) 

(— '*>  =  - 1 - 

bw^  ka(l  +  8a3)(  1  +  Il,f  )+2Gk(l  +  8as  -211*  +  8a3  IT  *  )  +  8IT  3  F(  1  +2a) 


[k(^  -  j)  -  G][na(l  +  8as  )] 

{k[2F(i-  +  a®  )  +  (^-|^-a  +2asF)n?]  -  Ilf  G(1  +  2a)}  {k2  -2Gk  -  8F} 

o  Z  Z  V  V 

(I-lla) 


U 


_ Gk  - ka  <f  -  r1 _ 

bu  *  k[2F(^+aa  )+(|--~  -  a  +  2aaF)na]  -IlaG(l+2a) 


(I-llb) 


To  solve  Eqs.  (1-11),  a  trial  and  error  process  must  be  used. 
Briefly,  this  involves  assuming  a  value  of  the  reduced  frequency  k, 
which  in  turn  prescribes  the  F  and  G  values,  solving  for  the  flutter  speed 
and  flutter  frequency  coefficients,  and  then  combining  the  coefficients  so 
as  to  obtain  a  computed  value  for  the  reduced  frequency.  A  valid  solution 
is  obtained  when  the  calculated  value  of  the  reduced  frequency  is  the  same 
as  the  assumed  value. 

An  examination  of  Eq.  ( I  —  1 1 )  shows  that  the  non-dimensional  static 

unbalance,  x  ,  and  the  non-dimensional  radius  of  gyration,  r  ,  do  not 

a  a 

appear.  The  exclusion  of  these  two  non-dimensional  system  parameters, 
coupled  with  the  mass -density  ratio  being  fixed  at  zero,  allowed  a  rather 
extensive  parameter  study  to  be  undertaken.  Even  though  the  computational 
process  remained  lengthy,  it  was  possible,  by  using  high  speed  digital 
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equipment,  to  check  literally  hundreds  of  combinations  of  II  2  and  (a)  over 
a  range  of  reduced  frequencies  which  effectively  covered  the  entire  unsteady 
region. 

Practical  considerations  dictated  that  limits  be  placed  upon  the 
values  of  Ilf  and  (a)  to  be  investigated.  Clearly,  only  values  of  Ilf  ^  0 
have  any  physical  meaning,  accordingly  only  positive  values  were  studied. 
With  regard  to  the  elastic  axis  location,  (a)  one  might  conclude  that  only 
values  of  (a)  between  plus  and  minus  one  need  be  considered.  However, 
for  reasons  which  will  be  explained  later,  the  limits  chosen  were 
-4  £  a  ^  3.  Actually,  no  lengthy  discussion  of  parametric  combinations 
is  necessary  at  this  time  -  it  being  sufficient  to  remark  that  none  of  the 
combinations  investigated  yielded  real  finite  flutter  points  However,  it 
was  noted  that  certain  of  the  combinations  of  IT  2  and  (a)  did  appear  to  be 
approaching  flutter  solutions  as  the  reduced  frequency  was  increased,  the 
trend  continuing  even  for  extremely  large  values,  ^on  the  order  of  1  03  ). 

For  any  practical  physical  system  b  and  must  be  real,  finite,  and 
non  zero.  Hence,  the  above  noted  behavior  implied  that  certain  of  the 
combinations  appeared  to  have  a  zero  flutter  speed.  Since  a  zero  flutter 
speed  is  associated  with  a  reduced  frequency  of  infinity,  it  was  concluded 
that  the  combinations  yielding  a  zero  flutter  speed  could  be  deduced  from 
the  equations  resulting  from  the  use  of  the  approximate  qua  si -unsteady 
aerodynamics;  Eq.  (I  —  1  0 ) . 

Examination  of  Eq.  (I-10a)  reveals  that  a  zero  flutter  speed  results 
if: 

qa  +  (  1  +  8a3)  [n*(  1  -  q)  -  q]  =  0 

where  q  =  1-  [2  (2a  -  1  )3  IT3  +  ( 1  +  8a“  )]  (1-12) 

Substituting  for  q  and  re-arranging  leads  to  an  expression  which 
can  be  identified  as  a  quadratic  in  Ila.  The  resultant  expression  may  be 
written  as : 
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Obviously  only  positive  values  of  the  frequency  ratio  Ilf  are  of 
practical  interest;  hence,  it  may  be  seen  that  positive  values  of  Ilf  will 
only  result  for  -  <  a  <  Thus,  for  any  elastic  axis  location  between 

the  three-eighth  and  the  three-quarter  chord  positions,  there  exists  a 
frequency  ratio  Ilf  which  will  result  in  a  zero  flutter  speed  solution  at 
a  mass -density  ratio  of  zero.  Figure  3  is  a  plot  of  all  such  combinations. 
As  a  check,  it  was  verified  that  the  combinations  obtained  by  using 
Eq.  (1-15)  included  those  that  were  indicated  by  the  numerical  analyses 
using  the  "exact11  unsteady  hydrodynamics.  In  effect,  the  "exact"  analy¬ 
sis  showed  that  the  flutter  points  obtained  by  using  the  quasi -steady  theory 
associated  with  a  reduced  frequency  of  zero  could  not  be  accepted  as  valid 
while  the  quasi-unsteady  theory  (k  -*  ®)  had  indeed  been  correct  in  not  pre¬ 
dicting  any  finite,  real,  non -zero  solutions.  Some  insight  as  to  why  the 
quasi -unsteady  theory  should  prove  more  reliable  than  the  quasi-steady 
(k  -*  0)  theory  may  be  gained  by  examining  Eq.  ( I  - 1 1 ) . 
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A  glance  at  Eq.  ( I  —  1 1 )  reveals  that  the  expressions  for  both  the 
flutter  speed  coefficient  and  the  flutter  frequency  coefficient  include  the 
reduced  frequency  (k).  It  may  be  further  noted  that  the  highest  power  of 
k  to  be  found  in  each  expression  is  never  combined  with  the  term  G,  G 
being  the  imaginary  portion  of  the  Theodorsen  function. 

The  behavior  of  the  Theodorsen  function  as  a  function  of  the  re¬ 
duced  frequency  is  well  known,  a  simple  plot  being  included  in  Appendix  A. 
Briefly,  for  the  quasi-steady  case  of  k  -*  0  F  =  1  and  G  =  0.  As  k  is  in¬ 
creased,  F  tends  toward  a  value  of  -g  while  G  approaches  its  largest 
(absolute)  value  of  approximately  -.2  for  a  k  of  approximately  .2.  Fur¬ 
ther  increases  in  k  cause  G  to  again  rapidly  approach  zero  in  such  a 
fashion  that  for  a  k  of  6.  0,  the  value  of  G  is  approximately  02. 

Returning  again  to  Eq.  (1-11),  it  is  clear  that  for  large  values  of 
the  reduced  frequency  (kb  the  terms  including  the  higher  powers  of  k  will 
be  dominant.  Since  these  higher  powers  are  never  combined  with  the  terms 
containing  G,  neglecting  the  already  small  values  of  G  by  setting  G  equal  to 
zero  could  logically  be  expected  to  have  little  effect  upon  the  character  of 
the  resultant  solutions.  On  the  other  hand,  for  very  small  values  of  the 
reduced  frequency,  the  higher  powers  of  k  tend  to  become  insignificant. 
Thus  the  inclusion  of  the  correct  value  of  GP  however  small,  might  be  ex¬ 
pected  to  have  a  correspondingly  greater  effect.  Hence  one  would  expect 
to  find  the  results  based  upon  the  quasi-steady  theory  unreliable  for  even 
small  departures  from  zero  k  values;  while  the  quasi -unsteady  results 
would  be  expected  to  remain  valid  over  a  much  larger  range  of  reduced 
fr equenc  ie  s  . 

The  results  discussed  so  far  may  be  summarized  on  an  extremely 
simple  plot.  On  Figure  4,  the  =  0  case  appears  as  the  ordinate  of  the 
set  of  axes . 

By  an  extensive  numerical  search  it  has  been  established  that  no 

reasonable  or  practical  combination  of  U  3  and  (a)  exists  which  will  yield 

/  * 

a  finite,  positive,  non-zero  U^/bu  &  .  Hence  no  flutter  points  exist  along 
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the  ordinate  between  zero  and  infinity.  Similarly,  the  numerical  search 

coupled  with  a  study  of  the  quasi -unsteady  flutter  equations  have  shown 

that  certain  combinations  do  yield  a  U,/bu  *  of  zero.  Thus  solutions  exist 

f  a 

at  the  origin.  To  complete  the  study  requires  only  that  the  possibility  of 

a  U  /b(J  of  infinity  be  considered;  that  is,  the  possibility  of  the  zero  mass- 
f  a 

density  ratio  ordinate  being  an  asymptote.  This  possibility  can  be  investi¬ 
gated  by  examining  Eqs.  (1-9),  which  are  derived  from  the  quasi-steady 
approximation,  inasmuch  as  these  equations  become  exact  as  the  flutter 
speed  coefficient  approaches  infinity. 

It  may  be  seen  from  Eq.  (I-9a)  that  a  flutter  speed  coefficient  of 
infinity  requires  that: 


g  +  (2a  +  1  )FI  *  -  0 


(I-I6a) 


where:  g  =  4a(2a  -  1  )IT®  +  (1  +  8as) 

Substituting  and  collecting  terms  yields  a  quadratic  in  (a)  of  the 

form: 


as(  1  +  Il*)8  -  2all*  +  (1  +  II*)  =  0  (I-I6b) 

No  real  values  of  (a)  will  satisfy  Eq.  (I-l6b)  if  II®  is  restricted  to 
positive  values.  Thus  no  physically  possible  combination  of  IT®  and  (a) 
can  be  expected  to  asymptotically  approach  a  flutter  solution  along  the 
[j  =  0  ordinate . 

Expressed  in  the  terms  of  Figure  4,  the  only  valid  flutter  solutions 
found  for  the  zero  mass -density  case  have  been  found  to  be  at  the  origin  of 
the  axes,  and  these  have  required  certain  carefully  prescribed  combina¬ 
tions  of  n?  and  (a). 

'fi 
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Solution  Gradients  as  the  Mass -Density  Ratio  Becomes  Finite 

Clearly,  to  be  of  any  practical  significance,  the  zero  flutter  speed 
solutions  discussed  in  the  previous  section  would  have  to  be  the  terminal 
points  of  valid  flutter  boundaries.  To  express  the  situation  graphically, 
consider  Figures  5a  and  5b.  Note  that  to  avoid  imaginary  terms,  the 
ordinate  has  been  designated  the  square  of  the  flutter  speed  coefficient. 

If  the  zero  flutter  speed  solutions  are  associated  with  flutter 
boundaries  of  the  types  depicted  by  Fig.  5a,  they  are  of  practical  interest, 
Conversely,  boundaries  of  the  type  depicted  by  Fig.  5b  are  of  little  sig¬ 
nificance.  To  determine  whether  any  of  the  configurations  depicted  by 
Fig.  3  are  associated  with  the  types  of  boundaries  shown  in  Fig.  5a 
requires  that  at  least  a  representative  sample  of  the  Ilf  and  (a)  com¬ 
binations  be  investigated  at  mass -density  ratios  other  than  zero. 

Unfortunately,  investigations  conducted  for  non-zero  mass-density 
ratios  require  that  all  of  the  non-dimensional  system  parameters  be  in¬ 
cluded.  Thus  the  computational  problem  would  once  again  appear  to  be  a 
major  obstacle.  However,  since  the  anticipated  boundaries  should  depart 
from  the  zero  flutter  speed  point,  it  is  reasonable  to  expect  that  for  at 
least  a  limited  range  of  mass -density  ratios  the  quasi -unsteady  equations 
associated  with  a  k  of  infinity  will  satisfactorily  predict  the  trend  of  the 
solutions.  In  addition,  for  very  low  mass-density  ratios,  the  effects  of 
the  structural  mass  and  inertial  parameters  will  be  minimal.  Thus  the 
initial  departure  of  the  solutions  from  the  origin  of  the  axes  should  be 
only  slightly  dependent  upon  the  values  of  x^  and  r^  ,  so  long  as  the 
numerical  values  are  physically  reasonable.  Clearly  then,  the  intended 
method  of  investigation  requires  that  the  flutter  equations  based  upon 
first  the  quasi-unsteady  and  then  the  "exact"  unsteady  hydrodynamics 
be  developed.  Both  can  be  easily  obtained  from  Eqs.  (1-8),  the  funda¬ 
mental  equations  of  motion. 
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By  replacing  the  C(k)  terms  in  Eq.  (1-8)  by  ^  and  then  solving  the 
associated  flutter  determinant,  the  quasi -unsteady  flutter  equations  take 
the  following  form: 


U  3  +  2u  ( 1  -  4a  +  4a3  +  8ax  -4x  +  4  r  3  ) 

a  a  a  a 

CJfs  1  +  8a3  +  2n  3(1  -  4a  4-  4a3)(l  +  jj)  +  g  8  r^3 

(1-17) 


(1  +  |i)(l  +  8as  +U8roS)[  II  * 


(1  +  n.,2)] 


u. 


<b-*)= 


O' 


CJ  4 

+  f  1  +U  [1  +  8as  +  1  6ax  -u8x  3  +8r  S(1  +  u  )]} 

cj  Ci  M  a  a 

a 


4nf(l  +U)(1  +2a)  +  4  — T  1  -u(l  +2a  +  2x  )] 

'r>  *  a  a 

LJ 

a 


In  a  similar  fashion,  the  flutter  equations  incorporating  the  "exact” 
unsteady  theory  may  be  obtained  by  substituting  the  complex  Theodorsen 
function  (F  +  jG)  into  Eq.  (1-8).  These  results  may  be  written  as  shown 
on  the  following  page,  (Eqs.  1-18). 


k[2 -F  +2U(1  -F -2a  +4a2Fx  +8aFx  -  2x  ,+  4Fr  2  )]  +  4G[  1  -  U 1  +  2a  +  2x  )] 


(q8I-l)  PUB  (B8I-I)  -sb3 
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If  Eqs.  (1-17  and  18)  are  compared  with  their  respective  counterparts  for 
the  zero  mass-density  ratio  case,  Eqs.  ( I  - 1 0  and  11),  the  increase  in 
complexity  is  seen  to  be  considerable.  However,  the  remarks  made  for 
the  u  =  0  case,  Eq.  (I - 1 1 ) ,  regarding  the  grouping  of  the  k  and  G  terms 
are  seen  to  remain  valid  for  the  more  general  case  represented  by 
Eq.  (1-18).  Thus  the  conclusions  reached  regarding  the  reliability  of 
the  solutions  obtained  from  Eq.  (1-10)  still  apply  to  Eq.  (1-17). 

Having  obtained  the  necessary  relationships,  the  investigative 
procedure  consisted  of  taking  a  number  of  the  combinations  depicted  by 
Figure  3  (chosen  so  as  to  cover  the  complete  range  of  possible  elastic 
axis  locations),  substituting  these  values  into  Eqs.  (1-17),  and  noting 
the  resultant  flutter  frequency  and  flutter  speed  coefficients.  Values 
of  x^  and  r^  were  chosen  so  as  to  more  than  cover  the  range  normally 
encountered,  in  spite  of  the  fact  that  the  expected  effect  of  these  para¬ 
meters  at  very  low  mass -density  ratios  was  expected  to  be  slight; 
while  the  range  of  values  investigated  went  from  zero  to  one  in  such 
a  fashion  that  very  small  step  sizes  were  taken  near  zero  with  the  step 
sizes  increasing  as  jj  increased.  Obviously,  the  computational  effort 
involved  was  considerable  even  though  high  speed  digital  equipment 
was  used.  Again  no  detailed  discussion  of  the  results  is  required  as 
no  valid  flutter  speed  solutions  were  found.  In  general,  the  flutter 
frequency  coefficients  were  positive  and  on  the  order  of  unity.  How¬ 
ever,  the  resultant  squares  of  the  flutter  speed  coefficients  were  all 
negative  -  indicative  of  imaginary  flutter  speeds.  Graphically  the  re¬ 
sults  could  be  portrayed  by  a  plot  similar  to  Figure  5b. 

It  might  be  of  academic  interest  to  mention  that  the  configura¬ 
tion  consisting  of  Ilf  =  1  and  x  =  a  =  0  was  found  to  yield  a  zero 
flutter  speed  solution  for  all  jj. 

Having  obtained  the  results  based  upon  the  quasi-unsteady  theory, 
a  number  of  the  configurations  were  checked  by  the  "exact"  theory  as 
depicted  by  Eqs.  (1-18).  Because  of  the  increased  complexity  of  the 
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calculations,  the  number  of  configurations  checked  was  restricted. 
Further,  because  of  the  transcendental  nature  of  the  equations,  imaginary 
solutions  could  not  be  assumed  and  thus  the  "flutter  boundaries"  obtained 
from  the  approximate  theory  could  not  be  duplicated  However,  it  is 
noteworthy  that  no  positive  finite  flutter  points  were  obtained. 

To  summarize,  the  non-zero  mass  -density  ratio  investigations 
described  in  this  section  may  be  regarded  as  the  final  portion  of  the  work 
undertaken  in  an  effort  to  ascertain  what  type  of  flutter  behavior  is  to  be 
expected  near  a  zero  mass-density  ratio.  The  work  described  in  this 
section,  in  conjunction  with  that  described  in  the  previous  section  effec¬ 
tively  shows  that  for  a  two-dimensional,  fully  wetted,  unswept  lifting 
surface,  no  flutter  points  of  practical  significance  may  be  expected  for 
a  mass -density  ratio  of  zero.  Further,  it  has  been  shown  that  the  zero 
flutter  speed  solutions  obtained  for  the  limiting  case  are  not  the  terminal 
points  of  real,  finite  flutter  boundaries;  but  rather  are  associated  with  a 
family  of  imaginary  solutions  and,  hence,  of  only  academic  interest. 

I-D  -  The  Non-Zero  Ma s s  -Dens lty  Ratio  Analysis 
The  Fundamental  Approach 

The  fact  that  the  two -degree -of -freedom  typical  section  analysis 
for  the  limiting  condition  of  a  zero  mass -density  ratio  will  not  yield  any 
flutter  points  of  practical  interest  may  be  regarded  as  a  valuable  por¬ 
tion  of  the  complete  low  mass -density  ratio  investigation  as  it  allows 
the  conclusion  that  no  finite  flutter  speed  boundaries,  plotted  as  a  func¬ 
tion  of  the  mass -density  ratio  i x >  may  either  approach  or  intersect  the 
=  0  ordinate  of  Figure  1.  This  simplifies  considerably  the  remaining 
investigation. 

Obviously,  the  absence  of  a  finite  flutter  speed  at  =  0  for 
physically  plausible  configurations  does  not  preclude  the  existence  of  a 
flutter  boundary  for  small  but  non-zero  values  of  the  mass -density 
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ratio.  Unfortunately,  a  parametric  study  for  finite  values  of  the  mass- 
density  ratio  is  considerably  more  involved  than  for  the  (j  =  0  case  simply 
because  the  mass  and  inertial  properties  of  the  structure  must  be  con¬ 
sidered.  This  meant  that  every  effort  had  to  be  made  to  decrease  the 
computational  requirements  without  sacrificing  the  general  applicability 
of  any  conclusions  resulting  from  the  study.  Accordingly,  it  was  decided 
to  follow  a  program  designed  to  keep  the  computational  work  to  a  minimum 
by  changing  what  were  to  be  regarded  as  the  fundamental  non-dimensional 
system  parameters  and  also  by  obtaining  as  much  information  as  possible 
from  the  simplified  quasi -steady  and  qua  si -unsteady  flutter  equations. 

The  several  changes  in  the  system  of  non-dimensional  parameters  will 
be  discussed  first. 

Historically,  the  flutter  speed  coefficient  (U^/bu^)  and  flutter  fre¬ 
quency  coefficient  (u  J  u  )  are  considered  to  be  functions  of  the  five  sys¬ 
tem  parameters  u,  FI2,  a,  x  ,  and  r  .  For  the  non-zero  mass  density 

a  a 

ratio  portion  of  this  investigation,  the  non-dimensional  parameters  chosen 
a  re : 


(j  and  n2 


a 


x 

m 


r 

c 


the  mass -density  ratio  and  frequency  ratio 
squared  respectively;  same  definitions  as 
previously  given. 

the  non-dimensional  elastic  axis  location 

for  a  "typical  section";  more  generally  the 

location  of  the  axis  of  zero  static  coupling. 

Positive  for  axis  aft  of  midchord  . 

the  non-dimensional  distance  between  the 

midchord  and  the  center  of  gravity  =  x^  +  a. 

Positive  for  c.g.  aft  of  midchord. 

non-dimensional  radius  of  gyration  about  the 

1  T 

/  S  S  \"?  C  O’ 

center  of  gravity  =  (r  -  x  )  or  — *  5* 

0  "  mb2 


2 
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The  use  of  the  mass-density  ratio  g  requires  no  discussion.  How¬ 
ever,  the  very  fact  that  it  is  desired  to  keep  the  effects  of  the  mass-density 
ratio  entirely  separate  leads  to  the  choice  of  the  "in  vacuum"  frequency 
ratio  squared,  IIs,  and  subsequently  to  the  use  of  the  "in  vacuum"  pitching 
frequency  (u^)  in  the  expressions  defining  the  flutter  speed  and  flutter  fre¬ 
quency  coefficients. 

It  may  be  recalled  that  the  analysis  performed  for  the  zero  mass- 
density  ratio  required  that  the  frequencies  be  computed  for  the  structure 
submerged  in  the  surrounding  fluid,  (11^,  u  *,  and  (j*)*  No  disadvantage 
was  encountered  by  the  inclusion  of  the  mass -density  ratio  in  the  frequency 
parameters  as  the  entire  investigation  was  conducted  for  one  specific  value 
of  jj.  Conversely,  it  may  be  noted  that  requiring  that  the  "in  fluid"  fre¬ 
quency  ratio  remain  constant  over  a  range  of  mass -density  ratios  would 
be  to  imply  a  change  in  the  structure  itself;  the  change  being  a  function  of 
the  mass-density  ratio  change  involved.  Thus,  although  any  results  of 
such  an  investigation  would  certainly  be  mathematically  correct,  a  physical 
interpretation  might  be  somewhat  difficult  to  obtain.  On  the  other  hand,  by 
using  the  "in  vacuum"  frequencies,  and  thereby  keeping  the  mass-density 
ratio  entirely  separate,  it  is  a  simple  matter  to  trace  the  flutter  behavior 
of  any  reasonable  physical  configuration  through  any  range  of  values  of 
the  mass -density  ratio;  the  simple  physical  interpretation  being  that  of  a 
given  structure  submerged  in  fluids  of  varying  density.  The  two  approaches 
are,  of  course,  entirely  equivalent  everywhere  except  at  g  =  0.  Thus,  al¬ 
though  it  must  be  admitted  that  the  "in  fluid"  frequencies  are  probably 
more  representative  of  the  fluid -structure  system  taken  as  a  whole,  the 
desirability  of  completely  isolating  any  effects  due  to  changes  in  the  mass- 
density  ratio  requires  the  use  of  parameters  which  are  not  explicitly  de¬ 
pendent  upon  the  mass -density  ratio  in  any  way. 

The  remaining  three  dimensionless  parameters  are  completely 
determined  by  the  configuration  of  the  structure  which  the  section  is  as¬ 
sumed  to  represent.  The  parameter  (a)  is  the  classical  non-dimensional 
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location  of  the  axis  of  zero  static  coupling,  as  denoted  in  Figure  2,  and  is 
only  discussed  because  this  axis  does  not  necessarily  coincide  with  the 
elastic  axis  of  the  section  itself  Historically,  the  typical  section  has 
been  assumed  to  be  a  section  located  at  some  specified  spanwise  station 
of  a  straight  wing  which  could  be  analyzed  and  would  give  results  valid 
for  the  entire  wing.  Clearly,  the  correct  location  of  the  axis  of  zero 
static  coupling  for  such  a  representation  is  indeed  the  elastic  axis  of  the 
section  However,  there  are  several  important  classes  of  structures 
for  which  the  elastic  axis  of  the  section  may  not  be  regarded  as  the  axis 
of  zero  static  coupling.  For  example,  it  is  not  difficult  to  visualize 
certain  completely  movable  control  surfaces  for  which  the  axis  of  zero 
static  coupling  might  be  more  closely  related  to  the  pivot  point  than  to 
the  elastic  axis  of  a  similar  section  Further,  and  of  even  more  im¬ 
portance  insofar  as  this  investigation  is  concerned,  one  of  the  effects 
of  sweeping  a  cantilever,  when  a  is  taken  in  a  streamwise  plane,  is  to 
locate  the  axis  of  zero  static  coupling  at  other  than  the  section  elastic 
axis;  the  exact  location  being  a  function  of  the  spanwise  location  of  the 
representative  section.  Figure  6  shows  this  comparison  Hence,  it 
is  logical  to  expect  that  part  of  the  static  mechanical  effects  of  sweep¬ 
ing  a  structure  can  be  inserted  into  the  basic  section  analysis  by  the 
proper  choice  of  an  (a)  value.  This  is  precisely  why  the  range  of  (a) 
values  investigated  was  much  larger  than  that  required  to  cover  the 
section  proper,  (i  e  ,  range  for  numerical  work  was  -4  £  a  £  3 .  )  It 
was  felt  that  the  investigation  of  the  swept  cantilever  effects  (Part  II) 
could  be  simplified  and  the  results  clarified  by  the  resultant  separation 
of  the  causes  and  effects.  In  what  follows,  the  term  "elastic  axis" 
means  the  axis  of  zero  static  coupling  and  not  the  section  elastic  axis, 
unless  it  is  so  specified. 

The  decision  to  regard  the  elastic  axis  location  (a)  as  essen¬ 
tially  independent  of  the  shape  of  the  section  makes  it  desirable  to  sepa¬ 
rate  the  inertial  characteristics  from  the  elastic  axis  location.  Accord¬ 
ingly,  the  last  two  non-dimensional  parameters  adopted  are  x^  and  r^ 
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rather  than  their  classical  counterparts,  x  and  r  ,  since  the  latter  set 

a  a 

explicitly  includes  the  location  of  the  elastic  axis  (a).  Conversely,  the 

location  of  the  center  of  gravity  with  respect  to  the  midchord  (x  )  and 

m 

the  value  of  the  radius  of  gyration  about  the  center  of  gravity  (r  )  are 
obviously  independent  of  any  artificially  contrived  axis  of  zero  static 
coupling.  Again,  it  must  be  realized  that  either  system  could  be  used; 
the  choice  being  strictly  a  matter  of  convenience.  However,  a  clear 
understanding  of  the  physical  changes  in  the  system  implied  by  the  vari¬ 
ous  parametric  variations  is  not  ea  sily  obtained  if  the  parameters  are 
inter  -related 

From  a  numerical  standpoint,  the  value  of  r  8  may  be  regarded 

as  a  constant  for  most  reasonable  section  shapes,  a  value  of  .25  being 

representative.  (For  physically  sensible  section  shapes,  r  3  varies 

c 

between  approximately  .  20  and  .  33  depending  on  the  exact  shape  of  the 

section.  )  Since  previous  investigations  have  shown  the  flutter  behavior 

of  a  section  to  be  comparatively  insensitive  to  small  changes  in  the  radius 

of  gyration,  it  is  apparent  that  the  numerical  work  using,  m  general,  the 

single  value  of  r  3  may  be  regarded  as  representative,  thus  reducing  the 
c 

computational  effort.  However,  certain  qualitative  effects  of  varying  r  3 
are  discussed  later  in  some  detail 

Contrary  to  what  has  been  said  regarding  r  3,  the  location  of  the 
center  of  gravity  (x  )  is  a  strong  function  of  the  shape  of  the  section. 
Hence  the  choice  of  a  particular  value  for  x  ,  (c.g.  confined  to  section 
proper  or  -1  £  x  £  l)v  qualitatively  determines  the  actual  shape  of  the 
structure  being  investigated. 

Development  of  the  Equations  Used  for  the  Non -Zero 
Mass -Density  Ratio  Analyses 

As  previously  noted,  much  information  may  be  obtained  from  the 
approximate  quasi-steady  and  quasi-unsteady  results.  In  order  to  show 
how  the  approximate  theories  may  be  useful  in  determining  areas  where 
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finite  flutter  boundaries  exist,  the  flutter  equations  resulting  from  both 
the  approximate  and  the  '’exact”  unsteady  theories  must  be  presented. 

The  derivation  of  the  "exact"  flutter  equations  follows  naturally 
from  the  expansion  of  the  flutter  determinant  (Eq.  1-4)  as  illustrated  in 
Appendix  A.  The  approximate  flutter  equations  can  be  obtained  from  the 
"exact"  results  or  directly  from  Eq.  1-4  by  using  the  suitable  unsteady 
hydrodynamic  derivatives.  These  derivatives  are  listed  below. 


Regardless  of  the  method,  the  equations  to  be  presented  can  be 
obtained  without  any  conceptual  difficulty,  although  the  algebraic  effort 
required  is  considerable. 
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Case  1  C(k)  =  ^ 


id  8  li  (a8  -  a  +  -7  +  r  8  +  2ax  -  x  )  +  7 

a  ^3  _ _ 4 _ a _ a  a  8 

Ufa  uraa+un*(aa  -  a  +  i) 


(1-19) 


$Ai  (r  8  +  ur  8  +  8na  4~  +  a8na)  -  ft4|j8nar 
a  a  a  8  c 


U  a  -|j(ra+|ara-nxa  +  2ax  =  aa  +  ^-)  -  !• 

^  j  _xa  _  M  a  _a _ a _ 8  8 

b%  fiarua(i  -  r2®na) +%(i -!T2ana)+ ux  -  £] 

^  W 


To  incorporate  the  system  parameters  chosen  for  this  investiga 
tion,  the  following  substitutions  are  made: 

x  -  x  -  a 
a  m 


ra  =  ra+(x  -  a)a 
a  c  m 


This  results  in  Eqs,  (1-19)  taking  the  form: 


^rca  +  (Xm-^]+8 
|i  [rca  +  (xm  -a)a+  IIa(a  -  I)3] 


(1-2 Oa) 
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nan  [(r  a  +  [x  -  a  ]3)(  i  +  n  +  un3 )+  na(r+  a)] 
c  m  o 

-  nVna[r  a-(x  -  a)2-  (i  f[  r  3  +  (x  -a)»][l  +  u] 
cm  cm 


X3  = 


-  |i  (x  -  a)3  +  2ax  -  a3  +  h  -  i- 

m _ m _ o  8 _ 

S7a[ya(l  -r^n3) +-^-(1  -f2®n3)  +  yx 


(1-2  Ob) 


Since  X3  does  not  appear  in  Eq.  (I-20a)  it  follows  that  substitu¬ 
tion  of  Eq  (I -20a)  into  Eq.  (I -2 Ob)  will  result  in  an  explicit  expression 
for  the  flutter  speed  coefficient.  For  present  purposes,  it  is  only  neces¬ 
sary  to  remark  that  the  resultant  expression  may  be  written  as  a  ratio 
of  polynomials  in  jj  of  the  general  form 


X3  = 


£nn[f.  (r  3 *  x  ,  a,  n3)] 


£  M  F  Tg  (r  3  t  x  ,  a,  n3  )] 

o  pc  m 


(1-21) 


where,  for  any  given  typical  section,  the  coefficients  of  the  may 

be  regarded  as  constants 

Case  2.  C(k)  =  1 

Proceeding  as  in  Case  1,  the  following  quasi-steady  flutter 
equations  may  be  obtained. 
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r  3  -  a(|  -  a)  -  x  (|  -  2a)  +  i-u 

57s  =  — - - - - - 

r  3  -  a(i  -  a)II3 
a  * 

(1-22) 

nVra3n3  +  (nrff3  +|+a3)(u  +1  -  flPWn»  ) 

-  *(H+  1)  -  (|ix  -a)3 

^3  — 

2nVn3(|  +  a)  -  20?  [-|  +  M  (4  +  a  +  xft  )] 


By  replacing  r^3  and  with  r  3  and  x  ,  Eqs.  (1-22)  may  be 


re-written  in  the  form: 


U[r  3  +  x  (x  -4)  +£- 
c _ mm  8 

|j  [r3  +  (x  -a)3  +aIT3  (a  ~4)] 


(I-23a) 


&U9Jp[r  9  +  (x  -  a)3]  + 

[|J+  1  -  ^n3]  fU  [rca  +  (xm-a)3]+a3  +  i} 

-  ^®u(u+  l)[r  3+  (x  -a)3]-  (y(x  -a)  -a]3 

X3= - 2 - 22 - 22 -  (I-23b) 

2n4ulf(4+a)  203  [y(xm  +  4)  -  4] 

As  in  the  previous  case,  substitution  of  Eq.  (I-23a)  into 
Eq.  (I-23b)  results  in  an  explicit  expression  for  the  flutter  speed  co¬ 
efficient  which  may  be  regarded  as  a  ratio  of  polynomials  in  p(Eq.  1-21). 
Thus  the  flutter  equations  obtained  from  the  approximate  quasi-steady 
and  quasi -unsteady  hydrodynamics  are  seen  to  be  similar  in  form. 
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Case  3,  C(k)  =  F(k)  +  jG(k) 

For  the  case  where  the  "exact"  unsteady  theory  is  used,  the 

flutter  equations  may  be  written  as  shown  below,  (Eqs.  1-24).  Once 

again,  r  3  and  x  may  be  replaced  by  r  3  +  (x  -  a)3  and  (x  -  a), 
a  a  c  m  m 

respectively,  to  yield  Eqs.  (1-25),  as  shown  on  page  36. 


Eqs.  (I-24a)  and  (I-24b) 


+  a3  +  u(r  3  +  [  x  -a]3)}{k[2 -F  + u2(l  -2x  +  F[4x  3  +4r  3  -  1  ])] 

^  c  m  m  m  c 


36 


Eqs.  (I-25a)  and  (I-25b) 
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As  with  the  approximate  equations,  it  may  be  seen  that  X3  does 
not  appear  in  Eq.  (I-25a).  Thus,  in  principle,  Eq.  (I-25a)  may  be  sub¬ 
stituted  into  Eq.  ( I -2 5 b )  to  obtain  an  expression  for  the  flutter  speed  co¬ 
efficient  as  a  function  of  the  five  system  parameters,  ( IP  ,  a,  u,  x  ,  r  3  ), 

me 

the  reduced  frequency  k,  and  the  real  and  complex  parts  of  the  Theodorsen 
function,  (F  and  G  respectively),  which  are  themselves  functions  of  k. 
Obviously,  such  an  expression  may  only  be  regarded  as  valid  for  those 
combinations  of  the  system  parameters  and  the  reduced  frequency  for 
which  the  solutions  to  Eq.  I-25a  and  I-25b  can  be  combined  to  yield 

.  With  this  restriction  in  mind,  it  may  be  noted  that  the  re¬ 
sultant  expression  for  the  flutter  speed  coefficient  may  be  written  as  a 
ratio  of  polynomials  in  jj  of  the  general  form. 

r  [f  (r  a,  X  ,  a,  If,  k,  F(k),  G(k))] 

X3^ - 2__£ - 22 - -  (1-26) 

£uPrg  (r3,  x  ,  a,  II3,  k,  (F(k),  G(k))] 
o  P  cm 

Now,  unlike  the  equations  resulting  from  the  approximate  theories, 
the  coefficients  of  the  ^  cannot  be  regarded  as  constants  for  any  given 
typical  section  since,  in  general,  the  value  of  the  reduced  frequency  k 
obtained  from  solving  Eqs.  I-25a  and  I-25b  changes  as  the  mass  density 
ratio  is  changed.  Further,  it  is  possible  for  more  than  one  value  of  k  to 
satisfy  Eqs.  I-25a  and  I-25b  for  certain  sets  of  system  parameters.  Thus 
a  given  typical  section  operating  at  some  fixed  mass -density  ratio  could 
mathematically  have  more  than  one  flutter  speed  coefficient.  The  im¬ 
plications  of  this  last  possibility  will  be  discussed  in  greater  detail  later. 

All  of  the  equations  presented  in  this  section  are  useful  in  deter¬ 
mining  the  general  character  of  the  flutter  solutions  in  the  low  mass- 
density  ratio  region  and  will  be  discussed  when  their  specific  contribu¬ 
tions  are  considered.  For  present  purposes  it  is  sufficient  to  discuss 
the  general  forms  of  the  expressions  for  the  flutter  speed  coefficient 


38 


squared  in  order  to  justify  the  reasoning  used  in  limiting  the  types  of 
flutter  boundaries  which  might  possibly  occur  in  the  low  mass -ratio 
region.  The  reasoning  only  requires  that  for  real,  finite  values  of  the 
system  parameters,  the  square  of  the  flutter  speed  coefficient  be  a 
continuous  function  of  (j  everywhere  except  that  at  certain  discrete 
values  of  jj  the  denominator  may  vanish 

For  the  approximate  unsteady  hydrodynamic  theories,  the  flutter 
speed  coefficient  squared  has  been  shown  to  be  expressible  as  a  ratio  of 
polynominals  in  jj,  (Eq.  I  21);  where  for  any  given  typical  section  the 
coefficients  of  the  polynomials  may  be  regarded  as  constants.  Accord¬ 
ingly,  one  may  expect  the  expression  for  Xs  as  a  function  of  jj  to  be 
smooth  and  well  behaved  everywhere  except  at  those  points  for  which 
the  denominator  vanishes.  Hence,  a  plot  of  X2  versus  jj  may  be  ex¬ 
pected  to  be  smooth  and  continuous  with  the  possible  exception  of  certain 
discrete  values  of  u  for  which  Xs  may  be  infinite  Clearly  only  positive 
values  of  Xs  have  any  practical  physical  significance,  but  the  smooth  be¬ 
havior  of  the  expression  may  be  expected  regardless  of  the  practical 
significance  of  the  result 

For  the  "exact"  theory,  the  square  of  the  flutter  speed  coefficient 

has  also  been  shown  to  be  expressible  as  a  ratio  of  polynomials  in  , 

(Eq.  1-26).  However,  the  coefficients  of  the  ^  terms  may  not  be  regarded 

as  constant  for  fixed  values  of  If  ,  a,  x  ,  and  r  3;  but  rather  are  also 

m  c 

dependent  upon  the  reduced  frequency  k.  Nevertheless,  considering  the 
fact  that  F  and  G  are  continuous  functions  of  k  and  considering  the  form 
of  the  expression,  it  is  not  unreasonable  to  expect  a  plot  of  Xs  versus 
to  be  smooth;  again  excepting  those  points  for  which  the  denominator 
vanishes.  There  are  reasons  other  than  the  form  of  the  expression  for 
expecting  a  smooth  behavior.  First,  it  is  well  known  that  the  results 
for  moderate  values  of  the  mass -density  ratio  have  always  exhibited  a 
smooth  dependence  of  X  upon  ,  as  depicted  in  Figure  1  Second,  as 
has  been  shown,  the  results  obtainable  from  the  quasi-steady  and 
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quasi -unsteady  theories  must,  in  general,  be  smooth  and  well  behaved. 
These  approximate  theories  are  the  opposite  extremes  of  the  "exact11 
theory.  It  is  mentioned  again,  however,  that  it  is  possible,  for  certain 
combinations  of  the  system  parameters,  for  more  than  one  value  of  the 
reduced  frequency  to  satisfy  the  "exact"  flutter  equations.  This  implies 
that  for  certain  values  of  the  mass-density  ratio,  the  same  typical  sec¬ 
tion  may  yield  more  than  one  flutter  speed  coefficient.  Clearly,  from  a 
practical  standpoint,  only  the  lower  value  is  of  interest.  However,  the 
types  of  flutter  boundaries  which  might  result  because  of  these  multiple 
values  have  been  considered  and  are  discussed  in  later  sections. 

Possible  Types  of  Low  Mass-Density  Ratio  Flutter  Boundaries 

Basically,  the  arguments  of  the  previous  section  have  been  pre¬ 
sented  to  show  that  real,  finite,  flutter  boundaries  of  physically  plausible 
structures  (e.  g.  ,  Figure  1)  cannot  be  segments  of  curves  with  the  ter¬ 
minal  points  having  a  finite,  non -zero  value  for  either  the  flutter  speed 
coefficient  or  the  mass -density  ratio.  Further,  the  zero  mass -density 
ratio  analysis  has  shown  that  the  flutter  boundaries  can  neither  terminate 
at  nor  approach  asymptotically  the  line  depicting  a  zero  mass -density 
ratio,  (i.e.  ,  the  ordinate  of  Figure  1  or  Figure  7)  It  will  be  shown  that 
the  combination  of  these  constraints  serves  to  restrict  the  form  of  flutter 
boundaries  that  might  exist  in  the  low  mass -density  ratio  regime  to  a 
relatively  small  number.  Just  what  forms  these  boundaries  might  take 
can  best  be  illustrated  by  a  few  simple  plots. 

The  flutter  boundary  shown  in  Figure  7  is  representative  of  the 

type  of  boundary  known  to  result  from  a  two -degree -of -freedom  analysis 

of  many  aircraft  wing  "typical  sections.  "  The  line  labeled  jj  represents 

a 

the  value  for  which  the  quasi -steady  hydrodynamic  theory  predicts  an 

infinite  flutter  speed,  (Eq.pl).  Other  investigators  have  shown  (e.  g.  , 

Refs.  3,5),  that  for  so-called  "reasonable"  sections,  the  value  of  ^  is 

a 
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of  the  order  of  one.  Since  the  existence  of  finite  flutter  boundaries  for 

u  >  is  an  acknowledged  fact,  generally  the  region  of  greatest  interest 
a 

herein  is  jj  >  jj  ;  the  shaded  region  of  Figure  7  However,  one  possi- 
a 

bility  which  must  be  considered  is  that  the  boundary  depicted  in  Figure  7 

might,  for  certain  as  yet  uninvestigated  configurations,  be  displaced  in 

such  a  manner  that  the  flutter  free  region  might  be  extremely  small.  In 

view  of  the  works  of  previous  investigators  <’Refs.  3,5),  this  is  deemed 

unlikely.  However,  because  of  the  rather  extreme  section  configurations 

(e.  g.  ,  center  of  gravity  near  trailing  edge)  sometimes  used  for  hydrofoil 

struts  or  lifting  surfaces,  the  dependence  of  the  value  of  jj  upon  the  other 

a 

system  parameters  will  be  discussed  in  detail  in  the  next  section. 

Another  possibility  for  finite  low  mass -density  ratio  flutter 
boundaries  is  that  completely  separate  and  as  yet  undiscovered  boundaries 
exist  for  certain  reasonable  combinations  of  system  parameters.  In  view 
of  the  restrictions  previously  noted,  it  may  be  seen  that  any  "newM 
boundaries  confined  to  the  U  <  Ua  reglon  would  have  to  have  one  of  the 
four  general  forms  depicted  in  Figures  8a  ~d. 

Figures  8a  and  8d  illustrate  the  possibility  of  real  solutions  re¬ 
appearing  from  plus  infinity  and  either  tending  toward  a  zero  flutter  speed 
or  again  exhibiting  an  asymptotic  behavior.  Figures  8b  and  8c  show  solu¬ 
tions  reappearing  from  zero.  These  latter  boundaries  would  be  the  posi¬ 
tive  extensions  of  flutter  speed  coefficient  squared  (U^/bu^)3  solutions 
approaching  from  minus  infinity.  Further,  the  possibility  of  multiple 
solutions  implies  that  combinations  of  the  above  could  exist. 

Now  it  may  be  noted  that  since  the  jj  =  0  line  has  been  ruled  out 
as  a  possible  asymptote,  the  existence  of  boundaries  8a  and  8c  would  re¬ 
quire  that  the  solutions  approach  infinity  for  some  finite  value  of  the  mass- 

density  ratio  less  than  u  In  short,  for  the  boundaries  8a  and  8c  to  be 

a 

possible,  it  would  be  necessary  to  show  that  for  some  given  set  of  system 

parameters,  at  least  a  second  asymptotic  solution  (i.  e.  ,  a  second  u  ) 

a 

could  exist  for  some  value  of  0  <  |j  <  ^  .  This  means  that  the  requisite 

a 
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information  as  to  whether  such  boundaries  might  exist  may  be  obtained 
by  a  careful  study  of  the  quasi -steady ,  (k  -*  0),  flutter  equations. 

On  the  other  hand,  the  existence  of  boundaries  such  as  those  de¬ 
picted  in  Figures  8b  and  8d  would  require  that  valid  solutions  exist  on 
the  abscissa;  i.e  ,  valid  zero  flutter  speed  solutions  Hence,  these 
possibilities  may  be  checked  by  a  careful  study  of  the  quasi -unsteady, 

(k  -►<»),  flutter  equations. 

Thus,  a  detailed  study  of  the  approximate  flutter  equations  will 

reveal  whether  any  physically  plausible  configurations  may  be  expected 

to  yield  separate  and  distinct  flutter  boundaries  in  the  very  low  mass- 

density  ratio  region  -  i.  e  ,  for  some  jj  less  than  the  |j  predictable  for 

a 

that  configuration.  Accordingly,  the  approximate  flutter  equations  have 
been  carefully  investigated  These  investigations  are  discussed  in  de¬ 
tail  in  the  next  two  sections. 

The  possibility  of  the  "exact"  section  analysis  predicting  several 
flutter  solutions  for  some  range  of  mass -density  ratios  has  been  men¬ 
tioned.  Clearly  this  suggests  several  additional  types  of  flutter  boun¬ 
daries  that  might  exist.  Two  possibilities  are  depicted  in  Figure  9, 
where  the  "new"  boundaries  have  been  superimposed  upon  what  might 
be  termed  the  typical  representation  of  Figure  7 

In  view  of  the  results  of  the  modal  analyses  discussed  in 
Refs.  10,  11,  and  12,  the  boundary  tending  toward  a  flutter  speed  of 
zero  as  the  mass -density  ratio  approaches  zero  might  be  viewed  as  a 
particularly  attractive  possibility.  Nevertheless,  it  may  be  noted  that 
the  existence  of  these  "new"  boundaries  may  also  be  checked  by  exam¬ 
ining  the  approximate  flutter  equations  as  the  remarks  made  for 
Figures  8  again  apply. 

The  last,  and  probably  least  attractive  possibility  is  shown  in 
Figure  10. 

This  is  termed  least  attractive  for  several  reasons.  First,  it 
is  noted  that  previous  investigators  (e  g.  ,  Refs.  3  and  5)  have  shown 
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the  value  of  u  to  be  a  good  indicator  of  the  limit  of  the  "typical"  boundary 
for  all  the  cases  considered,  although  some  valid  solutions  have  been  ob¬ 
tained  for  values  of  jjvery  slightly  less  than  jj  .  Second,  it  is  apparent 
that  the  above  behavior  could  only  be  uncovered  by  using  the  "exact” 
flutter  equations  for  a  section.  Third,  the  resultant  boundary  could  not 
agree  with  the  low  mass -density  ratio  results  of  the  modal  solutions  per¬ 
formed  on  actual  three  dimensional  swept  structures  (e.g,  ,  Refs.  10,  11, 
and  12).  Nevertheless,  this  possibility  has  been  considered  and  is  eventu¬ 
ally  discussed  herein.  Essentially,  it  will  be  shown  that  all  the  evidence 
indicates  that  it  is  extremely  unlikely  that  any  boundary  such  as  the  one 
depicted  in  Figure  10  can  exist  in  the  low  mass -density  ratio  region  for 
reasonable  two  dimensional,  unswept  configurations. 

To  summarize,  it  has  been  shown  that  only  a  limited  number  of 
possibilities  have  to  be  examined  to  determine  whether  very  low  mass- 
density  ratio  flutter  may  be  predicted  by  a  two-degree -of -freedom  repre¬ 
sentative  section  analysis  Further,  it  has  been  noted  that  most  of  the 
possibilities  may  be  either  substantiated  or  discarded  by  a  careful  ap¬ 
praisal  of  the  approximate,  (k  -+  0;  »),  flutter  equations.  Accordingly, 
the  following  is  an  outline  of  the  work  tc  follow 

1.  The  approximate  flutter  equations  associated  with  the  quasi - 
steady  hydrodynamic  theory  (k  0)  are  discussed  in  detail. 

This  study  reveals  whether  any  physically  plausible  configura¬ 
tions  exhibit  an  asymptotic  approach  to  an  infinite  flutter 
speed  for  more  than  one  value  of  ^  and  hence  reveals  whether 
certain  of  the  flutter  boundaries  discussed  above  exist.  In 
addition,  the  dependence  of  |j  upon  the  various  structural 
system  parameters  is  determined  in  order  to  ascertain 
practical  minimum  mass -density  ratios  for  which  the 
"typical"  flutter  boundary  of  Figure  7  may  be  expected. 

2.  The  approximate  flutter  equations  associated  with  the  qua  si - 
unsteady  hydrodynamic  theory,  (k  -*  «)  are  discussed  in  detail. 

The  purpose  of  this  study  is  to  determine  whether  (and  if  so, 
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under  what  conditions)  zero  flutter  speed  solutions  can 
occur.  Thus,  the  existence  of  the  types  of  flutter 
boundaries  terminating  on  the  abscissa,  Figures  8b, 
c,  d,  and  9,  is  checked. 

3.  A  careful  comparison  of  the  asymptotic  tendencies  of 

the  quasi-steady  and  quasi -unsteady  solutions,  (jj  and 

a 

U  *  ),  is  made  -  with  particular  emphasis  on  the  low 
a 

mass -density  ratio  region.  This  comparison  indicates 

whether  flutter  boundaries  of  the  type  shown  in  Figure  10 

can  exist.  (A  flutter  boundary  of  the  type  being  considered 

logically  requires  that  ^  ’  be  appreciably  less  than  ^  ; 

a  a 

and  further  that  flutter  solutions  with  a  relatively  high 

reduced  frequency  be  expected  for  some  values  of 

U  <  U  •  )  Some  combinations  are  also  checked  using  the 
a 

"exact"  theory. 

The  discussions  of  the  above  investigations  and  the  conclusions 
resulting  therefrom  are  treated  in  the  remaining  portions  of  Part  I. 

Constraints  on  the  Flutter  Boundary  of  a  Section  Which  are  Imposed  by 
the  Character  of  the  Asymptotic  Flutter  Speed  Solutions 

Previous  sections  have  served  to  show  that  the  typical  section 

analysis  will  admit  to,  at  most,  a  limited  number  of  low  mass -density 

ratio  flutter  boundaries.  This  section  shows  how  a  careful  investigation 

of  the  dependence  of  the  ^  solutions  upon  the  various  non-dimensional 

a 

system  parameters  further  restricts  the  nature  and  the  number  of  pos¬ 
sible  boundaries;  the  fundamental  approach  being  a  study  of  the  qua  si - 
steady  flutter  equations. 

The  methods  of  obtaining  the  flutter  equations  associated  with 
a  reduced  frequency  (k)  of  zero  have  already  been  discussed  with  the  re¬ 
sultant  expressions  presented  as  Eqs.  1-22  and  1-23.  For  convenience, 
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the  numerator  and  the  denominator  of  Eq.  I -23b  are  divided  by  fi3  .  This 
leads  to  the  following  set  of  equations : 


n* 


U[r  s  +  X  (x  -  £)]  +  - 
_  c _ mm  8 _ 

|j[rs  +  (x  -af  +  aH^a-i)] 
K  c  m  *  J 


(I-27a) 


ft3,,8!!3  [r  s  +  (x  -  a)3  1  + 
M  L  c  m  J 


[ 


-  una]{U[rc8  +  (x 


m 


•H  +  a"  +  fj 


X3  = 


-u(u+  D[rca  +  -7^ru(xm-a)  -ai8 

Z&U&it  +a)  -  2[|i(xm  +#)  -*] 


(I -27b) 


For  a  physically  plausible  set  of  system  parameters,  the  above  equations 
become  "exact"  as  the  flutter  speed  becomes  infinite  or  the  flutter  fre¬ 
quency  zero.  Hence,  it  is  possible  to  obtain  from  the  above  equations  the 
value  or  values  of  the  mass -density  ratio,  expressed  as  a  function  of  the 
other  system  parameters,  for  which  the  flutter  speed  increases  without 
limit. 

Since  in  any  real  physical  system  an  infinite  flutter  speed  co¬ 
efficient  implies  an  infinite  flutter  speed,  it  is  obvious  that  an  infinite 
flutter  speed  results  if  either  the  numerator  of  Eq.  (I-27b)  is  infinite  or 
the  denominator  zero.  Clearly,  the  former  possibility  is  only  of  limited 
interest  as  it  immediately  requires  that  at  least  one  of  the  system  para¬ 
meters  be  infinite.  To  consider  any  of  the  system  parameters  as  infinite 
is,  in  general,  of  no  practical  interest  -  with  the  exception  of  the  section 
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frequency  ratio,  II2.  However,  it  may  be  seen  that  Eq.  (I-27b)  becomes  an 
indeterminate  expression  if  the  frequency  ratio  is  taken  to  approach  infinity. 
Hence,  the  flutter  speed  coefficient  becomes,  in  such  a  case,  a  function  of 
the  other  system  parameters  and  is  not,  in  general,  infinite  unless  the 
structure  is  regarded  as  having  no  torsional  stiffness,  (a)  =  0).  For  this 

latter  case,  the  resultant  infinite  flutter  speed  coefficient  is  seen  to  exist 
regardless  of  the  values  of  the  other  system  parameters.  This  may  be  re¬ 
garded  as  an  unfortunate  consequence  of  using  the  uncoupled  torsional  fre¬ 
quency  as  the  non-dimensionalizing  quantity.  Hence  it  may  be  concluded 
that  the  only  infinite  flutter  speed  coefficients  of  practical  interest  are 
those  that  result  from  the  vanishing  of  the  denominator  of  Eq.  (I-27b). 

If  Eq.  (I-27a)  is  substituted  into  the  denominator  of  Eq.  (I-27b) 
and  the  resulting  expression  is  set  equal  to  zero,  an  equation  involving 
the  mass -density  ratio,  jj,  to  only  the  first  power  results.  Solving  for  ^ 
in  terms  of  the  other  system  parameters  yields  the  value  of  |Jl  for  which 
the  flutter  speed  is  infinite,  \Jl  .  The  resultant  expression  is,  of  course, 

cl 

the  equivalent  of  Eq.  (1-1)  and  may  be  written. 


*(*m  -ap:[xmti.] 


(1-28) 


- 


m 


a)(aa+ a  -  j)  +  (^-  +  a)  [r  s  +(x  -a)8]} 

4  2  c  m  J  J 


Since  the  parept  expression  to  the  above  equation  contains  jjt  to 

only  the  first  power,  it  is  evident  that  for  any  single,  physically  possible 

combination  of  system  parameters,  only  one  U  will  exist.  Hence,  all  of 

a 

the  flutter  boundaries  discussed  in  the  previous  section  which  would  have 
required  at  least  a  second  asymptotic  solution  may  immediately  be  dis¬ 
carded  as  valid  possibilities.  Thus  the  boundaries  depicted  in  Figures  8a, 
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8c,  and  the  boundary  of  Figure  9,  which  shows  an  infinite  solution  for  some 
U  <  |i  may  all  be  eliminated  as  possibilities.  Further,  any  of  the  other 
types  of  boundaries  which  are  depicted  as  approaching  an  asymptote  must 
approach  the  one  determined  from  Eq.  (1-28). 

Since  it  has  been  shown  in  previous  sections  that  no  valid  flutter 
points  exist  for,  or  approach  asymptotically,  the  zero  mass -density  ratio 
line,  it  is  clear  that  if  some  combinations  of  system  parameters  can  be 
found  which  yield  a  negative  value  for  U^,  the  only  valid  finite  flutter  boun¬ 
daries  that  might  possibly  exist  are  those  shown  to  require  one  or  more 
crossings  of  the  abscissa.  In  the  next  section,  this  possibility  will  be  dis¬ 
carded.  Thus,  although  it  has  not  as  yet  been  definitely  proven,  it  will  be 
stated  at  this  time  that  an  infinite  or  a  negative  indicates  that  no  physi¬ 
cally  meaningful  flutter  solutions  exist  for  any  positive  mass -density  ratio. 
This  is  in  accordance  with  what  was  suggested  by  Herr,  (Ref,  5),  but  which 
unfortunately  was  not  proven.  This  entire  matter  will  be  discussed  in  de¬ 
tail  in  the  next  section,  which  also  includes  the  evidence  supporting  the 
above  statement. 

Eq.  (1-28)  may  be  studied  to  determine  the  dependence  of  upon 

the  various  structural  system  parameters.  That  cannot  be  zero  has 

already  been  established  and  discussed.  Hence  the  determination  of  what 

minimum  positive  values  may  be  expected  for  physically  reasonable 

systems  becomes  of  particular  interest.  In  addition,  the  dependence  and 

sensitivity  of  u  to  changes  in  the  various  structural  parameters  is  desired, 
a 

The  physical  constraints  imposed  for  this  investigation  are  listed 
and  discussed  below. 

1.  II8  must  be  equal  to  or  greater  than  zero.  Clearly  no 
meaningful  physical  systems  can  possess  negative  natural 
frequencies . 

2.  r  9  must  be  positive  and  will  be  restricted  to  values  less 

c 

than  or  at  most  equal  to  one  as  this  includes  all  meaningful 
physical  systems. 
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3.  The  center  of  gravity  will  be  constrained  such  that 

-1  £  x  £  1.  This  is  somewhat  arbitrary,  but  includes 
m 

all  systems  with  the  c.g.  on  the  section  proper. 

4.  The  limits  imposed  on  the  location  of  the  axis  of  zero 
static  coupling  are  again  arbitrary.  However,  for 
reasons  which  have  already  been  presented,  the  numeri¬ 
cal  analysis  will  consider  values  between  -4  £  a  £  3. 

The  dependence  of  the  calculated  values  upon  the  individual  non- 

dimensional  system  parameters  may  be  readily  established  by  investigating 

the  behavior  of  Eq.  (1-28)  as  all  but  one  of  the  structural  parameters 

(say  x  )  are  regarded  as  being  fixed.  Clearly,  this  poses  no  problems 
m 

However,  the  determination  of  a  positive  minimum  could  become  quite  in¬ 
volved  if  it  were  approached  in  the  classical  manner  whereby  the  above 
stated  constraints  would  be  imbedded  into  Eq.  (1-28)  and  the  partial  deriva¬ 
tives  of  the  combined  expression  examined.  In  addition,  the  boundaries  of 
the  mathematical  space  represented  by  the  function  and  its  constraints  and 
the  points  for  which  Eq.  (1-28)  does  not  possess  a  derivative  or  where  the 
function  is  not  continuous  would  have  to  be  examined.  Fortunately,  the 
classical  approach  is  not  required  as  the  positive  minimum  may  be  de¬ 
duced  by  combining  the  results  of  the  individual  dependence  analyses  with 
results  which  have  already  been  discussed.  That  the  deductions  are  cor¬ 
rect  is  proven  by  numerical  examples,  some  of  which  are  included  herein. 

A  glance  at  Eq.  (1-28)  will  show  that  it  may  be  put  into  the  follow¬ 
ing  form: 


A  +  IIs  B 
Ua  "  C  +  n8D 


(1-29) 


where  A,  B,  C,  and  D  are  functions  of  x  ,  r  a,  and  (a).  For  the  present 

m  c 

time,  these  may  be  regarded  as  constants.  To  determine  whether  any  IIs 

values  will  result  in  a  stationary  point,  the  derivative  of  U  with  respect  to 

a 

II8  is  set  equal  to  zero. 
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Accordingly,  the  following  may  be  obtained. 

BC  -  AD 

- —  =  -  =  0  (1-30) 

d(n3)  (c  +  ipd)3 

Therefore : 


a(ns) 

if  BC  =  AD  (for  all  finite  IT3  if 

DIP  +  C  does  not  equal  zero) 

or  if  IP  =  oo 

Thus,  there  exist  no  finite  values  of  II8  for  which  /g(IP)  is 

a 

zero  although  the  limiting  value  of  an  infinite  frequency  ratio  does  appear 

as  a  possibility.  The  first  solution  noted  above  is,  however,  of  some 

interest  as  it  yields  a  condition  for  which  the  value  of  jj^  is  independent  of 

IP  .  This  is  discussed  in  detail  later  in  this  section. 

Eq.  (I  -30)  also  shows  that  djj  / 3 ( IIs  )  approaches  infinity  as 

a 

DIP  +  C  approaches  zero.  This  would  yield  combinations  of  parameters 

for  which  (i  is  very  sensitive  to  changes  in  II8  . 
a 

With  regard  to  the  immediate  problem,  it  follows  from  the  above 

that  if  some  combination  of  (a),  r  8,  and  x  can  be  associated  with  a 

c  m 

minimum  positive  |j  ,  this  value  must  occur  for  one  of  the  limiting  values 
a 

of  IP  or  else  must  be  determined  by  the  behavior  of  the  function,  (Eq.  1-28), 
as  it  approaches  infinity.  This  latter  possibility  can  be  discarded  as  of 
no  practical  interest  since  an  asymptotic  behavior  of  U  as  a  function  of 

cl 

IP  could  only  yield  a  positive  minimum  if  a  curve  of  jj  versus  II3  passed 
through  zero  for  some  value  of  II8  -  in  which  case  the  minimum  U  would 
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be  zero.  However,  it  has  been  shown  that  cannot  be  zero  for  any  sys¬ 
tem  of  practical  interest.  Thus,  only  the  former  possibility  need  be  given 
serious  consideration. 

Clearly,  negative  values  of  IIs  are  of  no  practical  interest.  Thus 
the  limiting  values  of  IIs  that  need  be  considered  are  zero  and  plus  infinity 
The  case  of  IIs  equal  to  zero  is  easily  discerned  from  Eq.  (1-28)  and  seen 
to  yield  the  following  result. 


U 

a 

ns  =  o 


i 

1  +  2x 

m 


(1-31) 


This  is  a  particularly  simple  expression  showing  the  u  for 

a 

ns  =  0  to  be  entirely  dependent  upon  the  location  of  the  center  of  gravity. 
Now  to  ascertain  whether  some  positive  can  exist  which  is  less  than 
what  could  be  obtained  from  Eq.  (1-31),  the  other  limiting  value  of  II3 
equal  to  infinity  will  be  considered. 

If  Eq.  (1-28)  is  investigated  as  II8  approaches  infinity,  an  inde¬ 
terminate  expression  results  which,  when  evaluated,  leads  to  the  follow¬ 
ing: 


U 

a 

ns 


(a  -  x  )(as  +a  -■r)-(T+a)[rS  +  (x  -  a)3  ] 
m  4  2  L  c  m  J 


(1-32) 


Clearly,  the  above  expression  is  considerably  more  involved  than 
that  shown  in  Eq.  (1-31).  Since  it  contains  all  three  of  the  remaining  non- 
dimensional  parameters,  it  follows  that  the  minimum  value  which  Eq.  (1-32) 
might  achieve  is  a  function  not  only  of  the  center  of  gravity,  (x  ),  but  also 
of  the  location  of  the  elastic  axis  (a),  and  the  radius  of  gyration  of  the  sec¬ 
tion  (r  *).  Thus,  to  determine  the  extremums  of  Eq.  (1-32),  these  other 
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parameters  must  be  considered.  A  fruitful  approach  is  to  consider 
Eq  (I-3Z)  as  a  polynomial  in  (a)  and  arrange  the  expression  accordingly: 


U 

a 


a®(x 

m 


a®  a  1 

T  "  8  +T6 


-  a(x  a  +  r  a  +t)  +  (^ 
m  c  4  4 


r  s 
c 


Z 


(1-33) 


Now  the  extremums  of  U  as  a  function  of  (a)  will  be  sought. 

ana  — *  oo 

Again  the  extremums  may  occur  where  the  derivative  of  the  function  to  be 
extremized,  (U  ),  taken  with  respect  to  (a),  vanishes;  for  the  limiting 

ana-.ro 

values  of  (a);  or  may  be  determined  by  the  behavior  of  Eq.  (1-33)  as  it 
approaches  infinity.  Since  the  desired  quantity,  if  it  exists,  is  to  be  a 
minimum  positive  value,  the  third  possibility  may  be  neglected  for  the 
same  reasons  given  when  discussing  Eq.  (I-Z9). 

Since  (a)  is  only  arbitrarily  restricted,  the  limiting  value(s)  of 
(a)  may  be  taken  as  infinity.  Clearly  then,  the  second  possibility  would 
lead  to  an  indeterminate  form  for  Eq.  (1-35),  which,  when  evaluated,  will 
yield  the  same  expression  as  Eq.  (1-31).  Therefore,  for  very  large  values 
of  (a),  U  approaches  U  .  The  typical  numerical  solutions  of 

an®-a>  an®-o 

Eq.  (I-Z8)  which  are  included  herein  as  Figures  11a  through  1  lg  clearly 
show  this  trend.  Hence,  Eq.  (1-31)  is  seen  to  apply  for  still  another 
limiting  case.  Thus,  the  first  possibility  is  the  only  one  of  interest  from 
the  standpoint  of  any  "new"  positive  minimum 

Taking  the  required  derivative  of  Eq.  (1-33)  and  setting  it  equal 
to  zero  yields 


=  {—  -^-+  4z)(x  8  +  r  8  +  J- 2a*  -a)  +(a -h[aa(x  +h 

2  8  16  m  c  4  m  8Lm2 


■a(xm8+rc3  +?)+(Hr  -"T 


)]  =  0 


2 
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The  above  may  be  expanded  and  rearranged  to  yield  the  following 


expre  s  sion: 


x  x  8  r  8  x  x  8  r  a 

as/J2 L  JUL  JZL  c 

8  ~  2  2"  16'  8  ”  2  2 


x  x  8  r  3 

JZL__£ - L)  =  o 

4 '  8  2  2  16' 


which  may  be  expressed  in  particularly  simple  -  and  fortunately  separated 
form  as: 


X  x  ® 


(1-34) 


For  physically  plausible  systems,  the  expression  involving  the 
center  of  gravity  location  and  the  radius  of  gyration  cannot  be  zero.  Thus, 
the  required  extremizing  condition  may  be  obtained  by  setting: 


a8  +  a  -  —  =  0 
4 


(1-35) 


1  fz 

The  roots  to  Eq.  (1-35)  are  -  ~  ±  ~  and  represent  two  values  of 
(a)  for  which  Eq.  (1-35)  is  stationary.  The  result  is  particularly  interest¬ 
ing  because  it  is  independent  of  x  and  r  8.  Numerically  these  roots  are 
approximately  -1.2  and  +  .2. 

If  the  roots  of  Eq.  (1-35)  are  substituted  into  Eq.  (1-33),  the 
following  two  equations  result. 


1  /2 


*■  a  =  ‘2  -  — 


_ 

T  Xma  +  Xm(1  t;T)+ 1  +  !/2  +  7rc’ 


(I-36a) 
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(1-3  6b ) 


The  numerators  of  both  of  the  above  expressions  are  positive. 

However,  the  denominator  of  Eq  (I-36a)  cannot  be  negative  for  any  real 

x  and  any  r  3  >  0  while,  for  the  same  conditions,  the  denominator  of 
m  c 

Eq.  (1-3 6b )  can  never  be  positive.  (This  statement  may  be  checked  by 

regarding  the  denominators  as  quadratics  in  x  and  checking  the  dis  - 

m 

criminant  of  each.  )  Accordingly,  one  may  conclude  that  only  Eq.  (I-36a) 

can  yield  a  positive  and  thus  is  really  the  only  expression  of  interest. 

a 

The  second  expression  yields  a  minimum  negative  u  as  II3  -*<». 

a 

These  conclusions  may  be  confirmed  by  referring  to  Figures  11a 

through  llg,  from  which  it  may  be  noted  that  no  positive  values  exist 

a 

for  high  values  of  IIs  and  an  (a)  of  approximately  +  .2.  It  may  be  further 
noted  by  examining  Figures  11a  -  lie  that  the  positive  jj  obtained 

ana  » 

from  Eq.  (I-36a)  is  indeed  lower,  for  the  same  x  <  -i,  than  the  value 

m  * 

obtained  from  Eq.  (1-31)  for  the  limiting  IIs  =  0  case.  It  is  also  interest¬ 
ing  to  note  that  although  no  positive  U  could  be  obtained  for  the  IIs  =  0 

a 

case  with  the  center  of  gravity  ahead  of  the  quarter  chord,  (x  £  -i), 

m  ^ 

there  are  no  center  of  gravity  locations  for  which  the  extremum  of  the 

IP  -+oo  case,  (Eq.  I-36a),  will  not  yield  a  positive  U  -  provided  only  that 

a 

r  3  is  greater  than  zero.  Figures  1  If  and  llg  illustrate  this  fact. 

For  purposes  of  identification,  the  result  obtained  from  Eq.  (I- 36a) 

will  be  referred  to  as  (li  )  as  it  may  be  confirmed  from  the 

aip  _a)mm 

numerical  evidence  presented  in  Figures  1  la  through  llg  that  the  value 
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is  indeed  the  minimum  which  may  be  obtained  for  any  given  physically 

sensible  x  and  r  3.  By  examining  Eq.  (I-36a),  the  following  observa- 
m  c 

tions  regarding  the  dependence  of  this  minimum  u  upon  x  and  r  3  may 

a  m  c 

be  made. 

1.  (U  )  is  lowered  by  increasing  r  a.  However, 

a_  mm  7  s  c 


it  may  be  seen  that  for  the  range  of  r  3  values  normally 
encounte  red,  (.20  <r  3  <  .33),  the  effect  may  be  ex  - 


pected  to  be  slight. 


2.  (U 


)  is  lowered  by  increasing  x  positively, 
min  m  J 


3.  (U 


)  is  raised  by  increasing  x  negatively 
mm  m 


until  /Z  x  becomes  greater  than  x  (  1  4-/2/  2). 
m  m 

Beyond  this  point,  increasing  x  negatively  will  again 

decrease  the  extremum  value  in  such  a  fashion  that  the 

extremum  value  approaches  zero  for  very  large  positive 

or  negative  x  .  For  the  practical  range  of  -1  £  x  £  1 , 
m  m 

the  effect  of  center  of  gravity  location  is  monotonic  with 
a  forward  movement  of  the  center  of  gravity  tending  to 
increase  the  flutter  free  region 

The  extremum  value  can  only  be  made  to  approach  zero 

by  having  x  or  r  3  approach  infinity, 
m  c 

For  the  practical  limiting  value  of  the  center  of  gravity 
at  the  trailing  edge,  (x  =1),  and  with  an  r  3  equal  to 


25,  (U 


)  is  approximately  .  26. 
min 


Virtually  all  of  the  above  statements  are  shown  by  the  numerical 
results  presented  in  Figures  lla-llg.  However,  before  proceeding  with 
a  discussion  of  those  results,  it  is  interesting  to  return  to  Eq.  (1-30).  If 
the  BC-AD  of  Eq.  (1-30)  is  expanded,  and  rearranged,  the  following  may 
be  obtained: 
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For  physically  possible  systems,  the  first  factor  can  never  be 
zero.  Thus,  the  second  factor  must  vanish  for  the  above  equation  to  be 
valid. 

The  second  factor  may  be  written  as  a  polynomial  in  (a)  and  set 
equal  to  zero.  The  resultant  expression  may  be  written  as 


x  8  - 
m 


rcs  -  5>  *  ?<- 


x  S  "  r  3  ~  o) 

m  c  8 


(1-38) 


For  any  physically  possible  system,  the  expression  containing 

x  and  r  8  cannot  be  zero.  Thus,  the  only  valid  solution  to  Eq.  (1-43) 
m  c 

is:  (a)  =  -  ^  Hence,  an  (a)  of  -  ^  (the  elastic  axis  located  at  the  quarter 

chord)  represents  a  condition  for  which  u  is  independent  of  the  frequency 

a 

ratio,  IP  .  From  Figures  lla-llg,  the  validity  of  the  above  conclusion  is 
apparent. 

Figures  lla-llg  show,  as  has  been  stated  earlier,  the  results 

of  solving  Eq.  (1-28)  numerically.  Because  of  the  comparatively  weak 

dependence  of  LI  upon  r  8  for  typical  r  8  values,  (.  20  <  r  8  <  .33),  all 
r  a  c  c  c 

of  the  calculations  shown  are  based  upon  an  r  8  of  .  25. 

c 

It  is  easy  to  see  that  the  trends  shown  by  Figures  lla-llg  are  in 
accordance  with  the  analytical  results  which  have  been  discussed.  Further, 
the  combination  of  the  numerical  and  theoretical  results  allows  the  follow 
ing  statements  to  be  made.  It  is  again  emphasized  that  these  statements 
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are  based  upon  the  two-dimensional  typical  section  analysis  and  hence  may 
very  well  not  apply  for  planforms  which  cannot  be  accurately  represented 
in  such  a  manner. 

1.  For  configurations  which  have  the  elastic  axis,  (a),  forward 

of  the  quarter  chord,  u  will  decrease  as  IIs  increases.  For 

a 

(a)  aft  of  the  quarter  chord,  the  opposite  is  true. 

2.  For  a  given  c.g.  location  and  a  given  non-dimensional  radius 

-1  -/2 

of  gyration,  the  minimum  U  will  occur  with  an  (a)  of  — - - 

a 

and  IIs  of  infinity.  However,  it  may  be  seen  from  the  numeri¬ 
cal  results  that  the  (a)  required  for  a  minimum  U  is  only 

a 

slightly  dependent  upon  II®. 

3.  With  II®  =  0,  no  unstable  solutions,  regardless  of  the  mass- 
density  ratio,  may  be  obtained  if  the  center  of  gravity  is  for¬ 
ward  of  the  quarter  chord.  However,  for  the  case  of  II®  -*<», 
some  unstable  solution  can  always  be  obtained  if  the  elastic 
axis  is  near  the  leading  edge.  A  comparison  is  shown  in 
Figure  12. 

4  In  general,  aft  c.g.  locations  lower  |j  ,  (Figure  12).  How- 

a 

ever,  there  are  specific  exceptions  associated  with  the  maxi¬ 
mum  U  values  which  occur  for  moderate  values  of  II®  and 


a 

aft  elastic  axis  locations. 

5.  For  center  of  gravity  locations  forward  of  the  quarter  chord, 

(x  <  -.5),  no  positive  U  values  exist  for  (a)  locations  aft 
m  a 

of  the  quarter  chord  regardless  of  II®.  Thus,  it  is  extremely 
likely  that  configurations  with  a  c.g.  forward  of  the  quarter 
chord  and  an  elastic  axis  aft  are  flutter  free  for  all  U.  How¬ 
ever,  these  configurations  might  be  particularly  susceptible 
to  static  divergence. 

6.  Considering  only  the  flutter  aspects  of  the  problem,  it  is 
apparent  that  for  all  practical  values  of  x^,  certain  con¬ 
figurations  associated  with  comparatively  large  frequency 
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ratios,  (IT3  >1),  may  be  made  flutter  free  for  all  |i  by 
locating  the  elastic  axis  in  an  aft  position. 

The  above  statement  regarding  the  dependence  of  jj  upon  the 

a 

center  of  gravity  has  been  made  by  other  investigators  (e.g  ,  Refs.  3,  5, 

12).  Ref  5  also  included  a  figure  depicting  the  dependence  of  u  upon  x 

a  m 

and  (a),  for  a  moderate  II3  of  .  64  Unfortunately,  the  calculations  were 
performed  using  a  constant  r^3  with  the  result  that  many  of  the  configura¬ 
tions  depicted  were  physically  impossible,  thus  showing  a  questionable 

inte r -dependence  of  |i  upon  x  and  (a). 

a  m 

In  effect,  the  last  of  the  above  statements  says  that  for  a  given  x 

m 

and  r^3  ,  certain  combinations  of  II3  and  (a)  yield  a  negative  ua  which*  f°r 

reasons  which  have  been  discussed,  is  an  extremely  strong  indicator  that 

the  configuration  will  have  no  valid  flutter  points  for  any  positive  U.  Hence, 

it  is  desirable  to  show  for  what  parametric  combinations  a  negative 

might  be  expected  Obviously,  what  is  required  is  to  determine  what 

physically  possible  and  plausible  configurations  can  cause  |i  as  expressed 

a 

by  Eq.  (1-28)  to  change  sign.  Since  the  numerator  of  Eq  (1-28)  must  re¬ 
main  positive  for  meaningful  systems,  the  sign  of  U  must  be  determined 

a 

by  the  denominator;  the  boundary  being  those  cases  for  which  the  de¬ 
nominator  is  equal  to  zero. 

If  the  denominator  of  Eq  (1-28;  is  set  equal  to  zero,  the  following 
expression  may  be  obtained 

x  r  3  x  3 

na  [aa(x  +  j)  -  a(x  a  +  r  a  +  j)  +  -  -j—  ~  y— )} 

'•m2  mc44Z  2  J 


+  [r  a  +  (x 
c  m 


a)8  ]  [x  +  T]  =  0 

J  L  m  2J 


d-39) 


across  x 
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For  computational  purposes,  Eq.  (1-39)  may  be  re  -written  as: 


IP  = 


f r  8  +  (x 
L  c  m 


+  7l 


(1-40) 


-a8(x  +  —•)  +  a  ( x  8  +  r  8  4-  j)  - 

m  Z  m  c  4  4  2 


x  8 
m 


Now  if  r^s,  a,  and  x^  are  specified,  a  valid  boundary  point  will 
exist  if  the  II8  calculated  using  the  above  is  positive. 

Some  typical  boundaries  are  shown  in  Figures  13a~13i.  Fig¬ 
ures  13a -13i,  show  that  the  location  of  the  elastic  axis  can  be  very  im¬ 
portant.  Also  of  interest  is  the  complete  change  of  character  of  the 
boundaries  as  the  center  of  gravity  is  moved  forward  For  c.g.  loca¬ 
tions  aft  of  the  quarter  chord,  valid  flutter  solutions  may,  in  general, 
be  expected  for  some  positive  U;  the  exceptions  being  fairly  well  con¬ 
fined  to  configurations  with  an  aft  elastic  axis  and  a  frequency  ratio 
greater  than  one.  Further,  the  combination  of  a  c.g.  aft  of  and  an 
elastic  axis  forward  of  the  quarter  chord  will  always  result  in  un¬ 
stable  solutions  for  somejj.  Conversely,  if  the  c.g.  is  forward  of  the 
quarter  chord,  valid  flutter  solutions  may  not,  in  general,  be  expected; 
the  exceptions  being  confined  to  configurations  with  a  forward  elastic 
axis  and  a  non -zero  II8  .  Further,  the  combination  of  a  c  g.  forward 
of,  and  an  elastic  axis  aft  of  the  quarter  chord  is  always  stable  for  all 
mass-density  ratios.  These  conclusions  are,  of  course,  in  agreement 
with  what  was  indicated  in  Figures  l2a-12g.  It  is  of  interest  to  observe 
that  certain  all  movable  control  surfaces  have  a  comparatively  high  II8 
and  may  be  characterized  by  a  center  of  gravity  and  an  elastic  axis 
forward  of  the  quarter  chord. 

As  has  been  previously  stated,  the  primary  reason  for  studying 

the  expression  for  U  ,  (Eq.  1-28),  has  been  the  desire  to  obtain  the 
a 

practical  limits  of  the  "typical"  flutter  boundary  depicted  by  Figure  7. 
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However,  many  of  the  remarks  concerning  the  connection  between  the 

value  obtained  for  u  and  the  existence  of  valid  flutter  solutions  have  im  - 
a 

plied  that  the  validity  of  the  flutter  speed  coefficients  as  determined  from 

Eq.  (I-27b)  is  solely  dependent  upon  the  algebraic  sign  of  the  denominator 

Clearly,  if  the  numerator  were  to  change  algebraic  sign  for  some  U  <  U  , 

a 

real  solutions  would  again  result  The  possibility  of  the  qua  si -steady 

equations  predicting  a  "new11  flutter  boundary  for  some  u<  |i  does  not 

a 

need  to  be  discussed  in  detail  as  other  investigators,  (e  g,  ,  Refs.  6,  12) 

have  shown  that  solutions  do  indeed  reappear  However,  as  has  been 

stated  earlier,  the  solutions  which  result  are  such  that  they  invalidate 

the  fundamental  quasi -steady  assumption.  Further  the  reappearance  of 

valid  solutions  for  U  <  |J  implies  that  one  of  the  boundaries  sketched  in 

a 

Figure  8  must  exist  The  proof  that  none  of  these  boundaries  is  possible 
will  be  completed  in  the  next  section  To  clarify  these  last  statements, 
an  analysis  similar  to  that  contained  m  Ref  12  will  be  presented  For 
computational  ease,  the  analysis  is  based  upon  the  assumption  of  II3  =0 
which  is  fairly  representative  of  many  structures  Further,  the  points 
discussed  above  are  very  well  illustrated  by  this  example. 

If  IIs  is  taken  equal  to  zero,  Eqs.  (I ’27a  and  27b)  may  be  written 


U[r  8  +  x  (x  -  I>]  +  I 
u  L  c _ mm  i  _  8 

u  [  r  *  +  (x  -  a)a  ] 
c  m  J 


(1-4  la) 


Xa 


fi8  fu(|J4  U[rca  +  (xm 


a>*]} 


(1-4  lb) 
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Substituting  Eq.  (I-41a)  into  the  numerator  of  Eq.  (I  —4 lb)  yields: 

NUM  =Wxmk(xm  -£)  -4]  (1-42) 

From  Eq.  (1-42)  it  may  be  seen  that  for  g.  greater  than  (^^(x  -  i), 

the  algebraic  sign  of  the  numerator  will  be  the  same  as  the  sign  of  x 

m 

Conversely,  for  U  less  than  (ijr)/(x^-  J),  the  sign  of  the  numerator  will  be 
opposite  to  the  sign  of  x 

m 

Since  must  be  positive  for  a  meaningful  solution,  the  sign  of 
the  denominator  may  be  seen  to  be  governed  by  the  now  familiar  expres¬ 
sion: 


na=o 


2x  4  1 
m 


(1-31) 


Thus  the  denominator  is  greater  than  zero  for  and  less  than 

cL 

zero  for  U  <  M-  . 

a 

The  foregoing  information  allows  a  diagram  (Figure  14)  of  the 
regions  of  real  solution  to  be  drawn.  In  Figure  14,  the  shaded  portions 
indicate  the  regions  of  real  solutions. 

The  interesting  thing  to  be  noted  about  Figure  14  is  that  real  solu¬ 
tions  are  predicted  for  |i  <  if  the  center  of  gravity  is  aft  of  midchord. 
Largely  as  a  result  of  this  fact,  it  was  concluded  in  Ref.  14  that  unswept 
sections  with  aft  centers  of  gravity  might  have  valid  flutter  solutions  for 
g  <  g  and  that  these  solutions  would  probably  be  characterized  by  very 

3. 

low  reduced  frequencies.  This  was  indeed  an  attractive  possibility  as 
hydrofoils  are  characterized  by  aft  centers  of  gravity  and  further  the  low 
mass -density  modal  flutter  solutions  obtained  by  the  author  and  others, 
(Ref.  10),  were  characterized  by  very  low  reduced  frequencies  Unfor¬ 
tunately,  however,  the  conclusions  are  not  valid.  If  Eqs  (I-41a  and  41b) 

are  solved  for  x  >0  and  g  <  g  ,  positive  values  will  result.  However, 
m  a 


-  60  - 


iff23  and  Xs  are  then  combined  to  find  the  resultant  reduced  frequency,  it 
will  be  found  that  the  quasi-steady  (k  -»  0)  assumption  is  very  poor  every¬ 
where  except  for  jj  very  close  to  U  .  Thus,  the  results  can  only  be  re - 

a 

garded  as  another  example  of  the  inadequacy  of  the  quasi -steady  equa¬ 
tions  with  regard  to  the  prediction  of  valid  flutter  boundaries  in  the  low 
mass -density  ratio  region.  That  the  solutions  indicated  by  Figure  14  for 
U  <  cannot  exist  will  be  definitely  proven  in  the  next  section.  Further, 
it  will  be  shown  that  the  correct  trend  can  be  predicted  by  the  quasi¬ 
unsteady  (k  -►  co)  flutter  equations. 

Constraints  on  the  Flutter  Boundary  of  a  Section  Which  are  Imposed  by 
the  Character  of  the  Zero  Flutter  Speed  Solutions 

It  is  shown  in  this  section  that  a  careful  investigation  of  the  depen¬ 
dence  of  the  (mas s -density  ratio  for  which  a  system  has  a  zero  flutter 

speed  coefficient)  solutions  upon  the  various  non-dimensional  system 
parameters  completes  the  proof  that  a  given  section  has  only  one  flutter 
boundary,  the  "typical"  boundary  of  Figure  7.  This  requires  that  the 
quasi-unsteady  flutter  equations  be  investigated. 

The  methods  of  obtaining  the  flutter  equations  associated  with  a 
reduced  frequency  (k)  of  infinity  have  been  discussed  earlier  in  this  re¬ 
port  with  the  resultant  expressions  presented  as  Eqs.  1-19  and  1-20.  It 
is  convenient,  from  the  standpoint  of  the  work  to  follow,  to  divide  the 
numerator  and  the  denominator  of  Eq.  (I-20b)  by  ft3.  This  leads  to  the 
following  set  of  equations: 


Is  3 

Lt[r  a  +  (x  -  j)  ]  +  — 

c _ m  2  J  8 _ 

u[r  s+  (x  -  aP  +  IIa (a  -  h8] 

K  L  c  m  c  J 


(I-43a) 
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u[(rc3  +  [xm_  a3a  H1  +u  +una)  +  n»(I  +  a)  ] 


-n»^n«[rc*  +  (xm-a)»]  {[rc8  +  (xm-aP][l+u] 

Q 


1  ,  1 


X3  = 


-u(xm  -a)«+2axm  -  a*  +  - }  - -Q 


ua(l  -  fisns)  +  j(l  +  Ux  -ya-i 


(I-43b) 


For  a  physically  meaningful  set  of  system  parameters,  the  above 
equations  become  "exact"  as  the  flutter  speed  approaches  zero.  Hence, 
it  is  possible  to  use  the  above  equations  to  obtain  the  value  or  values  of 
the  mass -density  ratio,  expressed  as  a  function  of  the  other  system  para¬ 
meters,  for  which  the  flutter  speed  becomes  zero.  It  is  well  to  note  that 
in  the  above  system  of  equations,  the  IT2  8  expression  (Eq.  I-43a)  can  never 
be  negative  for  meaningful  physical  systems  as  (j,  r  3,  and  113  cannot  be 
negative.  Therefore,  the  validity  of  any  flutter  solutions  is  entirely  de¬ 
pendent  upon  the  algebraic  sign  of  the  X3  expression  (Eq.  I-43b). 

For  any  meaningful  physical  system,  a  zero  flutter  speed  co¬ 
efficient  is  coincident  with  a  zero  flutter  speed.  Strictly  speaking,  a 
zero  flutter  speed  coefficient  results  if  either  the  denominator  of 
Eq.  (1—4 3b )  is  infinite  or  the  numerator  zero.  Clearly  the  former  possi¬ 
bility  is  of  no  practical  interest  for  reasons  which  have  been  discussed 
in  the  previous  section  when  comparable  situations  were  encountered  with 
the  quasi-steady  equations.  Thus,  any  meaningful  zero  flutter  speed 
solutions  must  be  coincident  with  the  vanishing  of  the  numerator. 

It  may  be  recalled  that  certain  of  the  flutter  boundaries  previously 
presented  as  possibilities,  (Figures  8  and  9),  were  shown  to  originate  or 
terminate,  or  both,  with  a  zero  flutter  speed  coefficient.  In  effect,  the 
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existence  of  any  of  these  boundaries  would  require  that  for  certain  specific 
geometrical  configurations,  X8  ,  (Eq.  I-43b),  be  able  to  change  smoothly 
from  positive  to  negative  by  varying  |i.  Thus,  the  behavior  of  Xs  as  a 
function  of  \Jit  near  those  points  for  which  Xs  is  zero,  is  as  important  as 
the  determination  of  the  parametric  combinations  yielding  a  zero  flutter 
speed 

Before  proceeding  with  an  examination  of  the  general  quasi - 
unsteady  flutter  equations,  it  is  instructive  to  again  consider  the  case  of 
ns  =  0,  This  is  admittedly  a  special  case  and  indeed  certain  of  the 
specific  conclusions  are  not  applicable  to  the  more  general  case  of  a  non¬ 
zero  frequency  ratio.  Still,  the  character  of  the  flutter  boundaries  ad¬ 
mitted  by  the  limiting  n8  =  0  case  is,  as  will  be  shown  later,  quite  typical. 

If  IIs  is  assumed  to  approach  zero,  Eqs.  I-43a  and  43b  may  be 
written  as: 


a 


s  - 


1  s  3 

n[r  8  +  (x  -  j)  ]  +  -r 
^L_c _ m  Z  J  8 

n[r  8  +  (x  -  a)8  ] 
c  m  J 


(I-44a) 


X8  = 


fi8[n(rc8  +  )<1  +w)]  -w  (  rc8  +  x^)  - 


a  a  U  1 

r  3  -  —  -  — 

C  8  8 


n*Cu<*m4»-7] 


(I-44b) 


If  Eq.  (I -44a)  is  substituted  into  the  numerator  of  Eq.  (I -44b), 
the  following  expression  may  be  obtained. 


NUM  =  m8(x 

M  m 


1  8 


-U<x 


m 


1_ 

4 


(I-45a) 
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The  above  expression  may  be  recognized  as  a  perfect  square 
which  may  be  written  as. 


NUM  =  [u(x 

m 


(I-45b) 


From  this  expression  it  may  be  seen  that  the  numerator  has  a 
single  zero  which  is  obtained  when: 


U  -  U 


o  Zx  -  1 

IIs-  0  m 


(1-46) 


What  is  of  even  greater  importance  is  that  for  any  meaningful 
system,  the  numerator  of  Eq.  (I-44b)  can  never  be  negative.  Thus,  the 
algebraic  sign  of  Xs,  and  hence,  the  significance  of  the  result,  is  entirely 
determined  by  the  denominator  expression.  This  of  course  implies  that, 
at  least  for  the  II2  -  0  case,  none  of  the  flutter  boundaries  which  are  de¬ 
picted  as  originating  or  terminating  with  a  zero  flutter  speed  solution  are 
possible,  (Figures  8b,c,d,  and  9),  unless  the  change  of  sign  of  the  de¬ 
nominator  of  Eq.  (I -5 Ob)  is  coincident  with  the  vanishing  of  the  numerator. 

With  regard  to  Eq.  (1-46),  the  following  observations  may  be  made: 


U  >0  for  x 
o  m 


Lt  =0  for  x 
o  m 


U  <0  for  x 
o  m 


>  h 


-  l 
”  2 


<  i 


Therefore,  no  meaningful  zero  flutter  speed  solutions  can  exist 
with  a  center  of  gravity  forward  of  the  three  quarter  chord. 


-  64  - 


Now  the  denominator  of  Eq.  (I-44b)  may  be  considered.  Since  for 
all  meaningful  systems  the  flutter  frequency  coefficient  squared,  fi2  ,  must 
be  positive  and  non-zero,  the  zero,  and  hence  the  algebraic  sign  of  the 
denominator  must  be  determined  by  the  expression: 

U(x  +  £)-£  (1-47) 

This  is  seen  to  be  identical  to  the  expression  obtained  by  the 
quasi-steady  (k  -►  0)  theory  and  is  seen  to  yield  a  denominator  zero  -  or 
infinite  Xs  -  when: 


u  =  14 


ns=o 


2x  4  1 
m 


(1-31  and  48) 


Further, 


^a  >  0  for  x  > 
m 

1 1  -♦  ao  f or  x  — 


1  <  0  for  x  < 
m 


-i 

2 


,i 

2 


1 

2 


Thus  the  denominator  is  positive  for  u  >  U  1  and  negative  for 

cL 

U  <  U  1 .  In  addition,  it  may  be  seen  that  the  numerator  and  the  denominator 

cL 

will  not  vanish  simultaneously,  thus  assuring  that  none  of  the  "new"  flutter 
boundaries  previously  discussed,  (Figures  8b,8c,8d,  and  9),  exist. 

As  a  summary  of  the  examination  of  the  quasi -unsteady  (k  -►  ao) 
flutter  equations  for  the  limiting  case  of  II3  =0,  a  diagram  analogous  to 
that  shown  in  Figure  14  may  be  constructed.  This  is  included  herein  as 
Figure  15  The  shaded  portion  of  Figure  15  is  the  region  where  real  solu¬ 
tions  exist. 
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From  Figure  1  5  it  may  be  noted  that,  unlike  Figure  14,  no  valid 

solutions  are  predicted  for  U  <  U  .  Further,  since  all  separate  and  "new" 

a 

flutter  boundaries  have  been  shown  to  be  impossible,  it  is  obvious  that 

the  correct  trend  must  be  shown  by  Figure  15  and  not  by  Figure  14  In 

addition,  if  Eqs.  (I-44a  and  44b)  are  solved  for  x  >  -i-  and  y  >  u  1  ,  the 

m  a 

resultant  reduced  frequencies  will  tend  to  confirm  the  quasi -unsteady 

assumption  except  for  a  comparatively  small  U  range  near  U  f.  Since 

a 

U  1  =  U  ,  this  is  no  restriction  but  only  indicates  that  the  asymptotic  ap- 
a  a 

proach  to  an  infinite  flutter  speed  indeed  occurs.  As  might  be  expected, 

the  quasi-unsteady  assumption  is  particularly  accurate  for  x  >  ^  and 

m  * 

for  mass -density  ratios  approximately  equal  to  U. 

Since  all  separate  and  "new"  flutter  boundaries  have  been  dis¬ 
missed  for  the  case  where  IIs  -  0,  it  is  concluded  that  the  flutter  boundary 

is  of  the  "typical11  form  depicted  in  Figure  7.  If  x  >  A,  a  valid  positive 

m  * 

zero  will  exist  and  the  boundary  will  approach,  touch,  and  depart  the  zero 

flutter  speed  coefficient  abscissa  of  Figure  7  as  M  is  decreased.  Further 

decreasing  jj,  will  cause  the  boundary  to  approach  an  infinite  flutter  speed 

coefficient  for  jj  =  [j,  =  |i  1  .  No  real  valid  solutions  can  be  found  for 

a  a 

U  <  U  . 
a 

It  might  be  here  noted  that  the  inclusion  of  structural  damping  into 
the  analyses  would  prevent  the  calculation  of  a  zero  flutter  speed  which  is, 
of  course,  physically  impossible.  Hence,  the  zero  flutter  speed  points 
could  be  more  correctly  viewed  (from  a  physical  standpoint)  as  regions 
where  the  system  may  be  expected  to  have  a  very  low  flutter  speed. 

The  procedure  to  be  followed  in  analyzing  the  more  general  case  of 
a  positive,  non -zero  frequency  ratio  is  quite  similar  to  what  has  been  fol¬ 
lowed  for  the  II8  =  0  case.  Further,  the  conclusions  regarding  the  non¬ 
existence  of  any  "new"  flutter  boundaries  will  be  shown  to  be  valid  in  the 
general  case.  However,  as  with  the  quasi-steady  case,  the  introduction 
of  a  non -zero  frequency  ratio  results  in  expressions  involving  all  of  the 
structural  system  parameters. 
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To  obtain  the  conditions  which  result  in  a  zero  flutter  speed,  the 
numerator  of  Eq.  (I-43b)  may  be  set  equal  to  zero.  If  Eq.  (I-43a)  is 
substituted  into  the  numerator  of  Eq.  (I-43b)  and  the  result  equated  to 
zero,  the  following  equation  may  be  obtained 


[u3y  +  g  us  HucpU  +  ns)  +  (cp  +  £  )}{cp  +  e] 

-[/ nscp[usy*  +  4  uy  +  £■£■}  ~  u4rcs  (cp  +  e)s  (1-49) 

us 

-  (a3  (cp  +  6  )(cp  +  e  )s  ~  -g-(cp  +  e)8  =  0 


where: 

Cp  =  r  ^  ( x  -  ap 

v  c.  m 

J3  =  ns  (a8  +  i-) 

y  =  r  8  +  (x 
'  c  m  Z 

1  S 

e  =  n8  (a  -  ^ ) 


6  = 


2a  x  -  a3 
m 


It  may  be  noted  that  all  of  the  factors  in  Eq.  (1-49)  contain  the 
mass -density  ratio  jj  to  at  least  the  second  power.  Since  the  case  of 
|i  =  0  has  been  discussed  and  has  been  shown  to  be  of  no  interest  in  this 
particular  problem,  a  |j3  may  be  factored  out  of  the  above  equation.  The 
remaining  expression  can  be  arranged  as  a  quadratic  in  jj  as  follows: 
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Us{  (cp  +  c)(l  +  ns)5py  -  IlSySup  -  r^s  (jp  +  e)3} 

+  u((cp  +  e)[y(cp  +  ^  )+ g'cp(  1  +  ns  )]  -|-  ns(py  -  (cp  +  5 )(cp  +  e)3 } 

+  f-  (cp  +  £  Kcp  +  0)  -  ns^  _  I  (cp  +  e)s  =  0 

(1-50) 

Equation  (1-5  0)  was  solved  for  literally  hundreds  of  combinations 
of  the  non-dimensional  structural  parameters  and  some  of  the  results 
are  included  herein  as  Figures  l6a-l6o  and  will  be  discussed  later.  Of 
particular  importance  is  the  fact  that  the  numerical  solutions  virtually 
proved  that  Eq.  (1-5  0)  is  a  perfect  square  of  the  form 

[V  -  f(r  2,  x  ,  ns,  a)]  2  =  0  (1-51) 

u  c  m  J 

Because  of  the  numerical  work,  the  algebraic  proof  was  under¬ 
taken  with  the  result  that  after  what  can  most  conservatively  be  described 
as  a  great  deal  of  algebraic  manipulation,  it  was  shown  that  the  discrim¬ 
inant  of  Eq.  (I“50)  is  identically  zero.  Thus,  the  numerator  of  the  X2 
expression,  (Eq.  I-43b),  can  never  be  negative  for  any  physically  possible 
combination  of  the  system  parameters.  Hence,  the  algebraic  sign  of  the 
X2  solution  is  entirely  determined  by  the  sign  of  the  denominator  Further, 
a  given  combination  of  structural  parameters  will  only  admit  to  a  zero 
flutter  speed  solution  for  one  value  of  the  mass -density  ratio.  Since  the 
form  of  the  numerator  will  not  allow  the  X2  solutions  to  proceed  smoothly 
from  positive,  through  zero,  to  negative  -  or  vice  versa  -  it  may  be  con¬ 
cluded,  as  it  was  for  the  II2  =  0  case,  that  none  of  the  boundaries  requir¬ 
ing  a  zero  flutter  speed  coefficient  terminal  point  are  possible,  (Fig¬ 
ures  8b,  c,  d,  and  9).  Since  at  most  one  meaningful  ju  and  one  meaning¬ 
ful  fJi  can  exist,  at  most  one  physically  valid  flutter  boundary  can  exist 
o 

for  any  given  two-dimensional  section.  Thus,  it  must  be  concluded  that 
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this  boundary  must  be  of  the  form  depicted  in  Figure  7,  or  possibly 

Figure  10,  except  that  in  certain  cases  the  boundary  may  touch  the 

abscissa  for  some  discrete  value  of  the  mass -density  ratio. 

The  possibility  of  dual  flutter  speed  solutions  for  /i  <  /J  -  thus 

a 

resulting  in  a  flutter  boundary  of  the  type  shown  in  Figure  10  -  will  be 

discussed  later  in  this  and  in  the  following  section.  In  general  it  will  be 

shown  that  although  some  valid  solutions  are  possible  for  /i  <  /i  ,  it  is 

a 

very  likely  that  the  "undershoot,  11  at  least  for  a  two-dimensional  section, 
is  always  very  slight,  and  thus  of  very  little  practical  significance.  In¬ 
deed,  analyses  using  the  "exact”  theory  have  produced  no  examples  to 
the  contrary. 

Because  of  the  complexity  of  Eq.  (1-50),  the  dependence  of  n  up- 

o 

on  the  various  parameters  was  determined  numerically.  However,  before 

discussing  the  numerical  results,  it  is  well  to  remember  that  only  those 

H  solutions  associated  with  a  real  flutter  boundary  (X2  ^  0)  are  of  any 
o 

practical  interest.  Hence  the  denominator  of  Eq  (I-43b)  must  be  con¬ 
sidered  because  for  a  given  section,  the  can  only  be  a  point  on  the 
flutter  boundary  if  for  =  fj,  the  denominator  expression  is  greater  than 
zero,  (i.  e.  ,  for  ^  >£Z  *  ).  Clearly  then,  the  values  for  which  the  de- 

O  3. 

nominator  vanishes  are  of  interest. 

If  the  expression  for  ,  (Eq.  I~43a),  is  substituted  into  the  de¬ 
nominator  of  Eq.  (I-43b)  and  the  result  set  equal  to  zero,  the  following 
expression  may  be  obtained: 

M  (a  +  4-)fM  [r  8  +  (x  -  a)8  +  II8(a  -  — )  ] 

-  n«[n(rc»  +  im*  -*m  +  7>  +  5B 

+  M  (M  x  -  M a  -  i-)[r  8  +(x  -a)8  +n8(a  -  —)  ]  =0 

m  c.  u  c  m  c 


(1-52) 
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Since  a  \x  appears  as  a  factor  in  all  of  the  terms,  it  may  be 
factored  out  of  the  expression.  The  £t  factor  is  of  no  interest  as  it  may 
be  regarded  as  being  cancelled  by  a  similar  factor  in  the  numerator  ex¬ 
pression,  (see  Eq.  1-49).  The  remaining  denominator  expression  con¬ 
tains  \x  to  only  the  first  power,  thus  indicating  that  for  any  given  section, 
the  qua  si -unsteady  theory  will  only  admit  to  one  mass -density  ratio  for 

which  the  flutter  speed  is  infinite,  (ju  ’).  This  is  in  agreement  with  what 

a 

the  quasi-steady  theory  has  shown.  However,  it  will  be  shown  that,  in 

general,  ju  =/  ju 
a  a 

If  the  remaining  expression  is  solved  for  fx  ,  the  following  expres¬ 
sion  may  be  obtained: 


fX  -  [X  y  - 


r[rcs  +  (xm  -a)s] + 113 [■§““! +  76  ] 


a  [X  +  f][r  3  +  (x  -  a)3  ] 
L  m  2  J  L  c  m  J 


(1-53) 


For  meaningful  values  of  the  system  parameters,  Eq.  (1-52)  is 

greater  than  zero  for  fx  >  ix  '  and  less  than  zero  for  fx  <  fx  f  .  Thus, 

a  a 

since  the  numerator  is  a  perfect  square,  the  quasi -unsteady  theory  will 

admit  to  no  flutter  solutions  for  ix  <  fx  '  .  Further,  if  \x  >  fx  1  ,  it  is  a 

a  o  a 

point  on  a  valid  flutter  boundary  while  if  fx  <  M  1  it  is  associated  with 

o  a 

a  set  of  imaginary  solutions.  For  the  latter  case,  the  physically 
interesting  flutter  boundary  will  not  have  a  zero  flutter  speed  point  for 
any  mass -density  ratio. 

The  dependence  of  M  1  upon  the  other  system  parameters  may  be 

a 

determined  for  Eq.  (1-53)  by  following  the  same  procedure  as  that  which 
was  outlined  in  the  previous  section  and  used  to  analyze  the  ju  expres- 

2t 

sion,  (Eq.  1-28).  Because  of  the  complete  similarity  between  the  two 
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analyses,  only  the  results  that  may  be  obtained  by  a  study  of  Eq.  (1-53) 
will  be  presented  herein.  These  results  may  be  summarized  as  follows: 

1. p *  1  cannot  be  zero  for  any  meaningful  combination  of  non- 

a 

dimensional  structural  parameters. 

2.  If  the  elastic  axis  is  located  at  the  quarter  chord,  (a  =  -£), 

/i  1  is  independent  of  the  section  frequency  ratio,  IIs. 
a 

3.  The  minimum  positive  value  for  ju  1  occurs  as  the  frequency 

St 

ratio,  IIs  ,  approaches  infinity  and  for  an  elastic  axis  located 
one  quarter  chord  forward  of  the  leading  edge  (a  =  -1.5). 

4.  The  actual  minimum  numerical  value  taken  by  ju  1  is  depen- 

£L 

dent  upon  the  location  of  the  center  of  gravity  and  upon  the 
non-dimensional  radius  of  gyration.  The  dependence  is  ex¬ 
pressed  by  the  following  relationship. 

13 

Ma  =  - T1 -  (I'54) 

MIN  (x  +  f  r  +  r  3 
m2  c 

For  an  r^s  of  .  25  and  the  center  of  gravity  of  the  section  at 

the  trailing  edge,  the  above  expression  yields  a  ju  '  of 

a  min 

.25.  This  is  very  slightly  less  than  the  comparable  M  . 

a 

5.  Since  the  numerator  of  Eq.  (1-53)  cannot  be  negative  for 

any  physically  possible  section,  the  algebraic  sign  of /u  1 

a 

is  determined  by  the  denominator.  K  the  denominator  is 
set  equal  to  zero,  the  following  expression  may  be  obtained. 


(1-55) 
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The  positive  II8  values  obtained  from  the  above  for  combina¬ 
tions  of  a,  x  ,  and  r  8  define  a  boundary  for  which  u  1  is 
me  a 

infinite.  For  x  >  -4,  (c.  g.  aft  of  quarter  chord),  positive 
m  A 

fi  1  values  exist  for  frequency  ratios  less  than  the  frequency 
a 

ratio  for  which  H  1  is  infinite.  For  x  <  -4,  the  opposite  is 
a  m  ^ 

true . 


Strictly  speaking,  the  only  reason  for  computing  a  ^  1  is  to  deter  - 

a 

mine  whether  the  fj,  associated  with  the  given  section  is  associated  with 
the  real  flutter  boundary  of  the  section,  inasmuch  as  the  quasi-unsteady 
theory  (k  <»)  is,  in  general,  not  valid  as  U^/b^  approaches  infinity. 
However,  as  will  be  pointed  out  when  discussing  Figures  16a-l6o,  there 
are  cases  for  a  /u  1  <  n  where  the  flutter  boundary  may  tend  toward  the 

cl  cl 

lower  value  before  approaching  an  infinite  flutter  speed  at  A*  =  A*  . 

a 

Figures  16a-l6o  show  plots  of  fj,  ,  /u  ,  and  jj,  1  versus  the  elastic 

o  a  a 

axis  location  for  various  values  of  x  and  II^  .  For  all  plots,  the  value 

m 

of  r  8  is  considered  to  be  .25.  For  center  of  gravity  locations  aft  of  the 
c 

three  quarter  chord,  (x  >4),  and  low  values  of  the  frequency  ratio 

m  * 

(112  <  l)t  all  three  curves  behave  in  much  the  same  manner  with  /u  >  A*  1 

o  a 

and  fj,  for  most  elastic  axis  locations.  For  higher  112  values,  A*  tends 
a  o 

to  become  less  than  A*  1  and  M  for  many  of  the  important  elastic  axis 

a  a 

locations  between  plus  and  minus  one.  Other  characteristics  are  also 

evident.  It  can  be  seen  that  H  and  /i  1  both  equal  l/  (2x  +  1)  for  an 

a  a  m 

elastic  axis  at  the  quarter  chord  and  also  both  approach  the  same  value 

as  (a)  becomes  large  -  either  positively  or  negatively.  Likewise  H 

equals  l/(2x  -  1)  for  an  elastic  axis  location  at  the  three  quarter  chord 

m 

(a  =  and  seems  to  approach  the  same  value  as  (a)  becomes  large.  It 

is  interesting  to  note,  for  the  higher  frequency  ratios,  the  sensitivity  of 

H  to  elastic  axis  locations  in  the  vicinity  of  the  three  quarter  chord, 
o 

However,  it  must  be  mentioned  that  this  behavior  is  only  of  academic 

interest  as  it  may  be  seen  from  the  fJL  and  ^  1  curves  that  no  valid 

a  a 

flutter  solutions  exist  for  the  combinations  of  parameters  required  to 
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produce  such  an  extraordinary  behavior.  It  may  also  be  seen  from  the 

curves  that  the  minimum  positive  values  of  /Li  and  /Li  1  are  comparatively 

a  a 

weak  functions  of  the  frequency  ratio  and  occur  for  an  elastic  axis  loca¬ 
tion  which  can  be  generally  located  as  slightly  forward  of  the  leading 
edge.  On  the  other  hand,  the  minimum  value  of  /Li  is  a  much  stronger 
function  of  the  frequency  ratio  and  occurs  for  an  elastic  axis  location 
near  the  mid-chord.  Of  course,  the  minimum  /Li  is  only  of  practical 
interest  if  it  may  be  associated  with  a  family  of  real  flutter  solutions, 

{n  >  V  '  )• 

o  a 

For  centers  of  gravity  forward  of  the  three  quarter  chord  but  aft 

of  the  quarter  chord,  the  character  of  the  /Li  curves  is  seen  to  change 

while  the  curves  for  /Li  and  /Li  1  remain  essentially  the  same  as  for  the 

a  a 

more  aft  c.g.  locations.  For  a  IIs  of  zero,  /i  is  infinite  for  x  =  i. 

o  m  * 

However,  for  II3  >  0,  valid  positive  /Li  points  are  seen  to  exist  for  many 

elastic  axis  locations  forward  of  the  three  quarter  chord.  For  center  of 

gravity  locations  forward  of  the  three  quarter  chord,  the  frequency  ratio 

range  for  which  no  positive  /Li  values  may  be  obtained  is  seen  to  increase. 

As  an  example,  for  the  center  of  gravity  at  the  mid-chord  (x  =  0),  no 

m 

positive /Li  values  exist  for  II3  £  1.  Furthermore,  the  elastic  axis  range 
for  which  physically  interesting  /Li  values  may  be  obtained  is  seen  to  be 
extremely  limited;  the  range  being  generally  located  in  the  vicinity  of  the 
quarter  chord. 

If  the  center  of  gravity  is  considered  to  be  forward  of  the  quarter 

chord,  the  character  of  the  M  and  *  plots  changes  to  the  extent  that 

a  a 

valid  positive  values  only  exist  for  forward  (a  <  -£)  elastic  axis  loca¬ 
tions.  The  region  for  which  the  /Li  solutions  are  of  interest  remains 
limited  both  with  respect  to  the  frequency  ratio  and  the  location  of  the 
elastic  axis. 

Having  discussed  the  general  character  of  the  zero  and  asymptotic 
solutions  depicted  by  Figures  l6a-o,  it  is  of  some  interest  to  consider 
some  of  the  specific  information  that  may  be  obtained  from  the  curves. 
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For  /Li  >  /Li  1  ,  the  zero  flutter  speed  solution  is  a  valid  point  on  a  real 

O  cL 

flutter  boundary.  Hence,  a  configuration  may  be  expected  to  have  a  low 

flutter  speed  for  mass -density  ratios  in  the  vicinity  of  /Li  ,  f  /Li  >  /Li  1  . 

o  o  a 

For  /Li  just  slightly  less  than  /Li  1  ,  the  zero  flutter  speed  point  is 
o  a 

not  on  the  flutter  boundary.  However,  the  flutter  boundary  will  tend 

toward  the  zero  and  probably  attain  a  rather  low  flutter  speed  coefficient 

before  increasing  toward  an  infinite  flutter  speed  as  /Li  is  approached. 

a 

If  /Li  is  much  less  than  ^  \  it  is  likely  that  the  quasi -unsteady  theory 
will  be  of  little  value  in  accurately  predicting  the  flutter  behavior  of  the 
given  section.  Rather,  the  complete  '’exact*1  unsteady  theory  is  needed 
to  correctly  determine  the  flutter  behavior  at  most  mass-density  ratios 
with  the  practical  lower  limit  to  the  flutter  boundary  being  established 

by  ^  a  ’ 

Something  about  the  flutter  boundary  may  also  be  learned  by  noting 

how  much  larger  /Li  is  than  /Li  !  and  consequently  /Li  .  If  /Li  »  /Li  ,  the 
o  a  TTy.  t  a  o  a 


uf/ 

slope  of  the  flutter  boundary  ($  btt)  A/u  ^  ak°ut  /Li  may  be  expected  to  b 
moderate.  Thus  a  slight  difference  in  /Li  may  be  expected  to  have  only  a 


1  )  about  /u  may  be  expected  to  be 

dM  '  o 


slight  affect  upon  the  flutter  behavior  of  the  section.  Conversely,  if  /Li 

and  /u  are  nearly  equal,  an  extremely  sensitive  dependence  of  the  flutter 
a 

speed  upon  the  mass -density  ratio  may  be  expected  -  at  least  for  the  mass 

density  ratio  range  between  the  two  values.  Thus  flutter  predictions  made 

for  mass-density  ratios  near  /Li  and /Li  may  easily  be  in  considerable 

a  o 

error  unless  the  parameters  are  very  carefully  ascertained.  It  may  be 
noted  that  the  likelihood  of  such  an  extra  sensitive  behavior  increases  as 


the  frequency  ratio  is  increased. 

Little  has  as  yet  been  said  about  the  fact  that,  in  general,  /Li  =  /Li  ’ 

a  a 

Strictly  speaking,  the  asymptotic  approach  of  the  flutter  boundary  to  an 

infinite  flutter  speed  will  occur  at  /i  since  the  theory  from  which  /Li  1  is 

a  a 

derived  is  not  valid  as  k  -*•  0.  Hence,  if  /Li  ’  >  /Li  ,  the  approach  to  an 

a  a 

asymptotic  solution  at  /Li  may  be  expected  to  be  smooth  and  of  the  type 

a 

characterized  by  Figure  7.  In  this  case,  the  value  of /Li  1  is  unimportant. 

a 
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However,  if  /u  1  <  M  ,  and  further  if  /Li  is  sufficiently  close  to  /Li  and  /Li  1 
a  a  o  a  a 

that  the  quasi -unsteady  theory  may  be  expected  to  be  valid,  the  boundary 

may  tend  towaict  jJ.  1  before  approaching  infinity  for  /Li  =  /Li  .  Hence,  two 
a  a 

flutter  solutions  can  occur  for  some  range  of  /Li  < /Li  .  Checks  using  the 

a 

"exact"  unsteady  theory  have  confirmed  this  fact  and  an  example  is  in¬ 
cluded  herein  in  Figure  19.  However,  it  may  be  seen  from  Figure  19 
that  the  range  of  double  solutions  is  slight.  Further,  it  may  be  seen 

from  Figures  l6a-l6o  that  when  /Li  1  <  /lx  ,  it  is  generally  only  slightly 

a  a 

less;  this  statement  being  particularly  true  for  very  low  mass-density 

ratios.  Hence,  it  is  extremely  unlikely  that  real  flutter  solutions  will 

be  obtained  for  any  mass -density  ratios  significantly  lower  than  /Li  ^  - 

particularly  in  the  very  low  mass -density  ratio  region.  This  is  deemed 

to  be  another  indicator  that  the  type  of  flutter  boundary  depicted  by 

Figure  1  0  is  unlikely  In  addition,  as  has  been  previously  stated,  none 

of  the  many  "exact"  solutions  were  found  to  yield  any  counter-examples. 

Just  what  general  types  of  flutter  boundaries  may  be  expected  as 

a  result  of  the  relative  values  determined  for  /Li  ,  /Li  ,  and  /Li  *  may  be 

o  a  a 

shown  by  a  series  of  simple  diagrams,  (Figures  17a-17e). 

Combinations  other  than  those  depicted  by  Figures  17a-17e  are, 
of  course,  possible,  but  the  five  which  are  shown  include  most  of  the 
more  interesting  combinations. 

It  may  be  seen  from  Figures  l6a-l6o  that,  for  certain  completely 
realistic  systems,  the  mass -density  ratio  for  which  the  lowest  flutter 
speeds  may  be  expected  can  change  drastically  with  only  small  changes 
in  the  system  parameters.  This  statement  is  particularly  true  with 
regard  to  the  location  of  the  elastic  axis  for  all  but  the  very  lowest  fre¬ 
quency  ratios.  Hence,  it  is  obvious  that  general  statements  regarding 
the  dependence  of  the  flutter  velocity  upon  the  various  non-dimensional 
structural  parameters  cannot  be  made  unless  the  mass-density  ratio  is 
specified.  For  example,  it  was  stated  in  Ref.  3,  and  is  in  fact  an  almost 
universally  accepted  "rule  of  thumb,  "  that  a  rearward  movement  of  the 
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elastic  axis  is  generally  accompanied  by  a  lower  flutter  speed.  A  glance 

at  Figures  l6a-l6o  reveals  many  contradictions  to  the  above  statement. 

It  is  evident  that  the  flutter  boundary  of  a  given  section  is  not  only 

a  function  of  the  relative  values  of  M  ,  M  >  and  jx  but  is  also  dependent 

o  a  a 

upon  the  sensitivity  of  the  flutter  speed  coefficient  to  changes  in  the  mass- 
density  ratio,  (3X/3M  ).  Since  the  lowest  flutter  speed  coefficients  may 
logically  be  expected  to  occur  near  -  if  jU  >  jU  -  this  is  the  region 
that  is  discussed  herein.  Examples  of  some  flutter  boundary  calculations 
are  included  herein  as  Figures  18a -18h.  These  figures  are  all  based  upon 
the  approximate  quasi -unsteady  theory.  However,  only  those  solutions 
for  which  the  resultant  reduced  frequencies  are  calculated  to  be  greater 
than  six  are  included.  Hence,  the  departure  from  the  boundary  predict¬ 
able  from  the  "exact"  unsteady  theory  is  slight.  This  statement  has  been, 
in  fact,  confirmed  by  comparing  certain  of  the  boundaries  obtained  from 
the  "exact"  unsteady  theory  to  some  of  the  results  shown  in  Figures  18a- 
18h.  Figure  19  is  included  herein  as  a  typical  comparison. 

It  may  be  seen  that  almost  all  of  the  configurations  shown  in 
Figures  18a-18h  are  characterized  by  aft  locations  of  the  center  of  gravity. 
This  is  hardly  surprising  as  most  of  the  valid  n  points  are  associated 
with  aft  centers  of  gravity  -  particularly  for  structures  having  a  low  fre¬ 
quency  ratio,  113.  Figure  18h  is  the  only  exception  with  an  x^  of  zero. 

The  configurations  shown  in  Figure  18a  are  all  characterized  by 

a  zero  frequency  ratio.  Hence  jx  and  jx  are  independent  of  the  elastic 

o  a 

axis  location.  For  this  case,  it  may  be  seen  that  with  the  elastic  axis 

aft,  lower  flutter  speed  coefficients  may  be  expected  for  mass -density 

ratios  on  either  side  of  ^x  .  However,  as  11®  is  increased  to  even  the 

o 

nominal  value  of  .4,  the  dependence  of  the  flutter  coefficient  upon  the 
elastic  axis  location  is  seen  to  be  very  much  a  function  of  the  mass- 
density  ratio  being  considered.  Figures  18e-18g  show  the  same  trends 
for  a  slightly  more  forward  center  of  gravity. 
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Figure  18c  shows  a  set  of  boundaries  for  the  same  elastic  axis 
location  but  different  frequency  ratios.  It  is  particularly  interesting  to 
note  the  extreme  sensitivity  to  small  changes  in  the  mass -density  ratio 
that  is  exhibited  by  the  boundaries  of  the  higher  frequency  ratio  configura¬ 
tions.  Because  of  the  extreme  sensitivity  to  small  changes,  it  is  likely 
that  accurate  flutter  analyses  of  real  systems  characterized  by  high  fre¬ 
quency  ratios  would  be  difficult  to  perform  for  mass -density  ratios  near 
M  .  Admittedly  the  location  of  the  center  of  gravity  of  the  configuration 
in  Figure  18c  is  extreme  However,  an  unfortunate  choice  of  parameters 
could  result  in  a  similar  behavior  for  centers  of  gravity  considerably 
more  forward,  (e.g.  ,  Figure  18h).  Figure  18h  is  also  interesting  be¬ 
cause  it  shows  how  sensitive  a  low  mass -density  ratio  flutter  boundary 
can  be  to  small  changes  in  the  location  of  the  elastic  axis.  Again  the 
comments  regarding  the  analytical  difficulties  apply.  When  one  considers 
the  many  approximations  made  in  the  flutter  analysis  of  a  practical  struc¬ 
ture,  it  becomes  apparent  that  a  flutter  behavior  which  is  extraordinarily 
sensitive  to  small  changes  in  any  of  the  system  parameters  can  make  an 
accurate  prediction  of  the  flutter  speed  extremely  difficult  to  obtain.  It 
certainly  appears  that  some  all  movable  control  surfaces  might  exhibit 
just  such  a  sensitive  behavior.  This  fact  may  be  the  reason  why  theore¬ 
tical  attempts  to  analyze  observed  instabilities  have  generally  proven  to 
be  unsuccessful.  The  mere  fact  that  two-dimensional  unsteady  hydro - 
dynamic  coefficients  are  normally  used  immediately  suggests  a  potential 
source  of  considerable  error.  This  particular  possibility  will  be  dis¬ 
cussed  in  more  detail  in  a  section  to  follow. 

Results  Obtained  by  Utilizing  the  MExactM  Two-Dimensional  Unsteady 
Hydrodynamic  Theory 

The  methods  of  obtaining  the  flutter  equations  associated  with  an 
arbitrary  reduced  frequency  (k)  have  been  discussed  with  the  resultant 
expressions  presented  as  Eqs.  1-24  and  1-25.  Since  the  non-existence 
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of  any  "new"  flutter  boundaries  has  been  assured  by  the  results  obtained 
using  the  approximate  quasi -steady  and  quasi -unsteady  flutter  equations, 
the  use  of  the  "exact"  equations  is  confined  to  amplifying  the  conclusions 
drawn  from  the  aforementioned  analyses.  Because  of  the  transcendental 
nature  of  the  "exact"  equations,  it  is  obvious  that  only  a  reasonable  num¬ 
ber  of  numerical  checks  could  be  made.  It  may  be  stated  that  none  of  the 
numerical  checks  yielded  any  significant  contradictions  to  what  was  con¬ 
cluded  from  the  analyses  of  the  approximate  equations  -  particularly  with 
regard  to  the  flutter  free  regions.  In  addition,  it  may  be  stated  that  the 
"exact"  results  agreed  very  well  with  the  results  obtained  from  the  quasi¬ 
unsteady  theory  for  the  portion  of  the  flutter  boundary  characterized  by  a 
reduced  frequency  greater  than  six.  A  typical  comparison  of  the  flutter 
boundaries  shown  in  Figures  18a -18h  with  the  boundaries  predicted  by  the 
exact  theory  may  be  obtained  from  Figure  19. 

From  Figure  19  it  may  be  noted  that  the  approximate  theory  pre¬ 
dicts  a  more  conservative  flutter  boundary  everywhere  except  for  mass- 

density  ratios  very  near  (J  .  However,  it  may  also  be  noted  that  the 

3 

exact  theory  shows  that  two  flutter  solutions  are  predicted  for  some  mass- 

density  ratios  less  than  /u  but  that  the  "undershoot”  is  quite  small.  Cal- 

3. 

culations  using  Eqs.  1-28  and  1-53  show  that  for  the  combination  of  para¬ 
meters  considered  in  Figure  18,  ^  1  is  slightly  less  than  (X  . 

a  a 

It  may  be  recalled  that  the  quasi -steady  hydrodynamic  theory  in¬ 
dicated  that  for  aft  centers  of  gravity,  solutions  could  be  expected  for 

fj  <  ft  .  Indeed  Figure  19  shows  the  approach  of  the  "exact"  boundary  to 
a 

ti  to  be  from  the  side  of  lower  mass -density  ratios.  There  were  not 
a 

enough  boundaries  calculated  to  ascertain  whether  this  behavior  could 

always  be  expected  for  configurations  possessing  an  aft  c.g.  However, 

it  is  interesting  to  note  that  this  would  assure  compliance  with  the  quasi - 

steady  trends  predicted  to  be  applicable  for  mass -density  ratios  very 

near  without  requiring  the  existence  of  any  "new"  flutter  boundaries, 
a 


-  78  - 


Hence,  the  region  of  dual  solutions  may  be  expected  to  be  confined  to  a 

narrow  range  of  mass -density  ratios  near  M  .  Such  was,  in  fact,  the 

a 

case  for  all  the  conf igurations  actually  checked.  Thus  there  is  every 

indication  that,  for  all  practical  purposes,  /u  may  be  regarded  as  the 

a 

mass -density  ratio  below  which  the  two -degree -of -freedom  section 
analysis  will  not  predict  any  valid  flutter  solutions.  In  short,  it  is 
deemed  to  be  extremely  unlikely  that  any  physically  reasonable  com¬ 
binations  of  non-dimensional  structural  parameters  will  cause  the 
typical  section  analysis  to  yield  an  exaggerated  flutter  boundary  of  the 
type  shown  in  Figure  10. 

To  summarize,  it  may  be  stated  that  the  numerical  work  in¬ 
corporating  the  "exact1*  unsteady  flutter  equation  confirms  the  results 
of  the  approximate  theories.  Thus,  a  lengthy  discussion  of  many  of 
the  actual  numerical  results  is  unnecessary. 

On  the  Apparent  Contradictory  Aspects  of  the  Flutter  Behavior  of  a 
Typical  Section  as  the  Low  Mass -Density  Ratio  Region  is  Approached 

Results  included  herein  have  shown  that  some  of  the  lowest 

a 

values  -  and  fJL  values  -  are  associated  with  high  frequency  ratios  and 
forward  locations  of  the  elastic  axis.  This  would  seem  to  be  in  apparent 
contradiction  to  what  has  come  to  be  accepted  as  the  result  of  a  body  of 
aircraft  experience.  However,  a  simple  sketch  will  serve  to  show  that, 
in  fact,  no  contradiction  exists  and  the  configurations  conducive  to  flutter 
in  the  range  of  mass -density  ratios  peculiar  to  aircraft  can  logically  be 
extrapolated  to  yield  the  higher  fJ,  ^  and  /u  values.  This  fact  is  illustrated 
by  Figure  20. 

It  is  a  well  known  fact  that  the  flutter  speed  coefficient  will  vary 
approximately  as  the  square  root  of  the  mass -density  ratio  (constant 
dynamic  pressure  boundary  for  a  given  section)  in  the  range  peculiar  to 
aircraft.  Thus,  it  is  logical  to  expect  that  the  configuration  have  the 
lower  flutter  speed  in  the  aircraft  range  will  tend  toward  a  zero  or 
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minimum  flutter  speed  for  a  mass -density  ratio  higher  than  its  more 

stable  counterpart.  For  mass-density  ratios  less  than  JjL  ,  (or  the  M 

o 

which  corresponds  to  a  minimum  flutter  speed  if  /u  <  jU  *  ),  the  trends 

o  a 

may  be  reversed,  as  shown  in  Figure  20.  Thus,  the  seemingly  para¬ 
doxical  behavior  peculiar  to  the  low  mass-density  ratio  region  is  seen 
to  be,  qualitatively,  a  very  natural  extension  of  what  has  come  to  be 
accepted  on  the  basis  of  aircraft  experience.  However,  as  may  be  seen, 
it  is  an  extension  that  must  be  made  with  care  as  it  is  the  flutter  boun¬ 
daries  that  may  be  extrapolated  and  not  the  comparative  flutter  behavior 
at  some  fixed  mass -density  ratio. 

A  more  interesting  example  can  be  discussed  with  regard  to  the 
use  of  two-dimensional  unsteady  derivatives  in  a  flutter  analysis.  Air¬ 
craft  experience  has  shown  that,  in  general,  the  use  of  two-dimensional 
derivatives  results  in  a  predicted  flutter  speed  somewhat  lower  than  that 
predicted  if  three-dimensional  derivatives  are  used. 

Ref.  2  includes  some  experimental  and  theoretical  comparisons 
which  show  that  for  the  higher  mass -density  ratios,  the  theoretical  data 
are  definitely  conservative.  However,  for  the  lower  mass -density 

ratios  near  ju  ,  some  unconservative  predictions  are  noted.  In  particu- 
a 

lar,  what  could  be  regarded  as  the  experimental  /u  values  were  seen 

a 

to  occur  at  a  mass -density  ratio  less  than  the  /u  predictable  from  the 

cl 

two-dimensional  theory.  This  apparent  contradiction  to  what  has  come 
to  be  regarded  as  accepted  fact  can  be  explained  by  considering  the  basic 
flutter  determinant,  (Eq.  1-4). 

If  the  two  homogeneous  equations  from  which  the  flutter  deter¬ 
minant  is  derived  are  divided  by  the  mass -density  ratio  fJL  ,  the  flutter 
determinant  may  be  written  in  the  following  manner: 
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It  may  be  noted  that  in  the  above  expression,  all  of  the  unsteady 
derivatives  appear  divided  by  the  mass -density  ratio  and  further,  that  the 
mass -density  ratio  only  appears  in  combination  with  one  of  the  deriva¬ 
tives.  Hence,  if  it  is  assumed  that  each  of  the  four  derivatives  can  be 
modified,  to  account  for  finite  span  effects,  by  applying  the  same  cor¬ 
rection  factor  (<q),  it  is  clear  that  the  finite  span  effect  may  be  duplicated 
by  increasing  M  some  given  amount.  In  such  a  case,  the  effect  upon  the 
flutter  boundary  of  using  the  three-dimensional  coefficients  would  be  to 
displace  the  boundary  as  shown  in  Figure  21. 

Figure  21  shows  that  the  two-dimensional  theory  may  indeed  be 
expected  to  be  unconservative  for  certain  of  the  lower  mass -density 
ratios.  In  particular,  it  may  be  seen  that  the  /u  obtained  from  the  two- 
dimensional  theory  may  well  be  too  high  -  particularly  if  the  actual 
structure  has  a  fairly  low  aspect  ratio. 

It  is  interesting  to  consider  some  of  the  results  previously  noted 
from  the  standpoint  of  what  effects  a  finite  span  might  have  on  the  re¬ 
sults.  It  may  be  recalled  that  the  flutter  boundaries  of  certain  configura¬ 
tions  characterized  by  high  frequency  ratios  and  forward  elastic  axis 
locations  were  found  to  be  extremely  sensitive  to  mass -density  ratio 
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changes.  Since  it  has  been  shown  that  finite  span  effects  may  be  regarded 
qualitatively  as  mass -density  ratio  changes,  it  is  logical  to  expect  that 
systems  known  to  be  sensitive  to  small  mass -density  ratio  changes  may 
also  be  sensitive  to  finite  span  effects.  Thus  an  accurate  analysis  of 
such  systems  may  be  very  difficult  to  perform.  Since  the  extraordinarily 
sensitive  systems  are  more  characteristic  of  all  movable  control  surfaces 
than  main  support  or  lifting  surfaces,  these  remarks  are  even  more  per¬ 
tinent  as  control  surfaces  normally  have  a  comparatively  low  aspect  ratio. 
The  extraordinary  sensitivity  of  control  surface  type  structures,  coupled 
with  the  possible  finite  span  effects,  may  well  explain  why  certain  ob¬ 
served  instabilities  have  defied  theoretical  explanation.  An  obvious 
remedy  is  to  use  the  correct  three-dimensional  coefficients.  However, 
this  is  very  difficult  and  may  even  be  impossible  for  overly  sensitive 
systems.  Thus  a  series  of  analyses  covering  a  range  of  mass -density 
ratios  seems  to  be  in  order  for  control  surface  type  structures  designed 

to  operate  at  or  near  the  /i  of  the  system. 

a 

I-E  Conclusions  Resulting  from  Part  J. 

Before  proceeding  with  the  conclusions  resulting  from  Part  I  of 
this  investigation,  it  is  again  emphasized  that  the  statements  to  be  made 
apply  only  to  fully  wetted,  two-dimensional  sections  acted  upon  by  two- 
dimensional  unsteady  hydrodynamic  forces  and  moments  and  that  any 
extrapolation  of  the  conclusions  must  be  made  with  great  care.  Further, 

It  is  to  be  understood  that  the  structures  which  the  typical  section  are 
assumed  to  represent  are  unswept  in  the  sense  that  no  spanwise  varia¬ 
tions  in  the  vibratory  behavior  of  the  structure  are  considered.  That 
io,  t  n rW* In  of  the  inertial  -  or  mechanical  -  effects  of  a  swept  structure 
have  hcioii  considered  by  regarding  the  section  as  being  a  streamwise 
out  t  i 4 Ml  and  using  the  concept  of  an  axis  of  zero  static  coupling  rather 
titan  nt  a  sot  thm  elastic  axis.  However,  none  of  the  possible  effects 
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of  sweep  on  the  spanwise  variation  of  the  unsteady  hydrodynamic  forces 
and/or  the  spanwise  pitch  to  plunge  ratio  have  been  considered. 

It  may  be  recalled  that  other  investigators  have  indicated,  on  the 
basis  of  a  limited  number  of  experimental  and  numerical  results,  that 
lifting  surfaces  amenable  to  a  two-dimensional  section  representation 
appeared  to  have  a  minimum  mass -density  ratio  below  which  flutter 
could  not  occur,  (Refs.  3,5).  On  the  other  hand,  it  has  been  definitely 
established  that  certain  swept  cantilevered  struts  will  flutter  at  very  low 
mass -density  ratios  and  further  that  numerically  sophisticated  theoreti¬ 
cal  results  have  shown  that  no  low  fJ.  flutter  free  region  exists  for  such 
structures,  (Ref.  10)  This  fundamental  difference  is  in  contrast  to 
extrapolations  based  upon  aircraft  experience.  Hence,  this  investiga¬ 
tion  may  be  regarded  as  having  been  an  attempt  to  establish,  in  as 
general  a  sense  as  possible,  whether  the  swept  cantilevered  struts  and 
the  unswept  structures  which  may  be  represented  by  a  two-dimensional 
section  indeed  pose  separate  and  distinct  problems  with  regard  to  opera¬ 
tions  characterized  by  low  mass -density  ratios. 

In  the  light  of  the  above  statements,  undoubtedly  the  most  impor¬ 
tant  and  most  fundamental  conclusion  resulting  from  Part  I  of  this  in¬ 
vestigation  is  that  the  two-degree -of -freedom  section  analysis  will  not, 
for  all  practical  purposes,  admit  to  any  dynamically  unstable  solutions 
for  mass -density  ratios  below  what  may  be  termed  a  "critical"  value. 
This  "critical"  mass -density  ratio  may  be  closely  approximated  by 
determining  the  mass -density  ratio  for  which  the  approximate  quasi - 
steady  hydrodynamic  theory  associated  with  a  zero  reduced  frequency 

will  predict  an  infinite  flutter  speed,  (fJ.  ).  In  the  course  of  verifying 

a 

the  above  fact,  several  interesting  and  related  characteristics  of  the 
section  analysis  have  been  noted.  These  may  be  stated  as  follows; 

1.  For  any  physically  meaningful  geometrical  configuration, 
only  one  mass -density  ratio  exists  for  which  the  flutter 
boundary  will  asymptotically  approach  an  infinite  flutter 
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speed,  (m  ).  Further,  n  cannot  be  made  to  approach  zero 

3  3 

unless  or  r^8  3re  considered  to  3ppro3ch  infinity.  For 

sections  of  re3son3ble  geometry,  (me3ning  3n  r  8  of  between 

c 

.  2  3nd  .4  3nd  with  the  center  of  gr3vity  loc3ted  on  the  sec¬ 
tion  proper),  the  minimum  ju  is  3pproxim3tely  .25  3nd 

3 

occurs  with  3n  (3)  v3lue  of  approxim3tely  -1.2,  3nd  x  of 

m 

1.0,  3nd  for  very  high  values  of  IIs.  However,  the  minim  urn 
y  is  a  weak  function  of  the  frequency  ratio  n»  . 

St 

2.  The  numerator  of  the  expression  for  the  flutter  speed  co¬ 
efficient  (U^/bu^)  which  is  derivable  from  the  approximate 
theory  associated  with  an  infinite  reduced  frequency  (k  -4  ») 
may  be  written  as  the  product  of  the  mass -density  ratio  y 
and  a  quadratic  in  y  which  is  by  itself  a  perfect  square. 
Hence,  for  any  given  section  only  one  mass -density  ratio 
greater  than  zero  may  admit  to  a  real  zero  flutter  speed 
solution.  Further  no  meaningful  flutter  boundaries  may 
terminate  on  the  zero  flutter  speed  abscissa  (see  e.  g.  , 

Figure  7),  as  the  numerator  of  the  flutter  speed  expres¬ 
sion  cannot  change  algebraic  sign.  Hence,  the  physical 
validity  of  the  y  flutter  speed  solutions  is  determined  by 
the  denominator;  the  algebraic  sign  change  occurring  at 

that  value  of  y  for  which  the  denominator  is  zero,  (y  1  ). 

a 

In  the  low  mass -density  ratio  region,  y  1  and  y  are 

a  a 

nearly  equal  for  any  meaningful  combination  of  structural 

parameters.  Hence,  if  y  is  greater  than  y  ^  (and  thus  for 

practical  purposes  y  ),  it  is  associated  with  a  family  of 

a 

real  solutions  and  will  yield  a  mass -density  ratio  region 
for  which  a  given  section  can  be  expected  to  exhibit  a  very 
low  flutter  speed. 

3.  If  the  center  of  gravity  is  considered  to  be  confined  to  the 

section  proper,  y  ,  y  ',  and  y  are,  in  general,  all  lowered 
a  a  o 

as  the  center  of  gravity  is  moved  aft. 
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4.  If  the  frequency  ratio,  118,  is  zero,  fj.  ,  fJ.  ',  and  fj,  are 

a  a  o 

completely  independent  of  the  axis  of  zero  static  coupling, 

(a).  However,  the  value  of  (a)  quickly  becomes  important 
as  the  frequency  ratio,  n» ,  is  increased;  with  forward  axis 
locations  (negative  (a)  values)  tending  to  produce  the  lowest 
H  values.  It  may  be  noted  that  certain  systems  exhibit  an 

cL 

extreme  sensitivity  to  comparatively  small  changes  in  the 
axis  location.  As  an  example,  certain  configurations  em¬ 
ploying  a  forward  axis  location  and  moderate  or  high  fre¬ 
quency  ratios  (ns  approximately  equal  to  or  greater  than 
one)  will  yield  very  low  flutter  speeds  over  a  relatively 
restricted  range,  of  mass -density  ratios  while  the  same 
parametric  values  but  with  an  aft  axis  location  may  well 
not  yield  any  dynamically  unstable  solutions  for  any  mass- 
density  ratio. 

It  may  be  noted  that  statements  one  and  two  above  are  sufficient 
to  assure  that  for  any  meaningful  system,  "new"  or  distinct  flutter 
boundaries  cannot  exist  and  be  entirely  confined  within  the  mass -density 
ratio  region  below  M  a* 

The  above  conclusions  certainly  imply  that,  contrary  to  what 
aircraft  experience  would  dictate,  the  historically  insignificant  differ¬ 
ences  between  a  swept  cantilevered  lifting  surface  and  its  two  dimen¬ 
sional  "typical  section"  representation  are  indeed  extremely  important 
factors  in  determining  the  actual  low  mass -density  ratio  flutter  be¬ 
havior  of  the  surface  in  question. 
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PART  II  -  THE  LOW  MASS -DENSITY  RATIO  HYDROELASTIC  BEHAVIOR 
OF  A  UNIFORM  CANTILEVERED  LISTING  SURFACE 

II -A  Background 

It  was  shown  in  Part  I  of  this  investigation  that  the  simplified  two- 
degree -of -freedom  typical  section  analysis  will  not  admit  to  any  separate 
and/or  distinctive  flutter  boundaries  for  mass -density  ratios  significantly 
less  than  fJ.  .  Hence,  the  obvious  conclusion  which  may  be  drawn  is  that 

2L 

certain  of  the  low  mass -density  ratio  flutter  boundaries  previously  en¬ 
countered  both  theoretically  and  experimentally,  (see  e.g.  ,  Ref.  10), 
must  have  been  the  result  of  system  details  that  were  entirely  lost  in  the 
simplified  typical  section  representation.  Therefore,  even  though  the 
primary  purpose  of  this  investigation  is  to  provide  a  basic  understanding 
of  the  fundamental  mechanism  of  low  mass -density  ratio  flutter,  it  is 
clear  that  the  simplified  mathematical  model  required  for  such  a  purpose 
must  contain  more  of  the  salient  features  of  the  structure  and  the  result¬ 
ing  hydrodynamic  loadings  than  what  is  provided  by  the  typical  section 
representation. 

Before  continuing,  it  must  be  noted  that  the  above  statements  are 
not  to  be  interpreted  as  meaning  that  all  dynamic  instabilities  associated 
with  comparatively  low  mass -density  ratios  are  completely  beyond  the 
scope  of  the  simplified  theory.  Indeed  it  was  demonstrated  in  Part  I  of 
this  study  that  certain  classes  of  structures  which  can  be  well  repre¬ 
sented  by  a  typical  section  may  be  expected  to  encounter  dynamic  sta¬ 
bility  problems  for  comparatively  low  mass -density  ratios.  This  was 
shown  to  be  true  in  spite  of  the  fact  that  the  complete  flutter  boundary  of 
these  structures  could  be  described,  as  a  function  of  the  mass-density 
ratio,  by  the  "usual"  boundary  of  Figure  7.  As  was  indicated,  at  least 
some  of  the  analytical  difficulties  associated  with  these  structures,  for 
mass -density  ratios  close  to  fji  ,  may  well  have  been  due  to  the  extreme 
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sensitivity  of  the  flutter  solutions  to  small  changes  in  certain  parameters. 
The  attendant  results  were  that  the  effects  of  the  usual  approximations 
(e.g.  ,  two-dimensional  hydrodynamic  theory)  were  magnified  and  in  some 
cases,  certain  commonly  expected  trends  were  reversed. 

Nevertheless,  the  specific  low  mass -density  ratio  flutter  occur¬ 
rences  which  formed  the  background  and,  in  fact,  the  motivation  for  this 
entire  investigation  (see  e.g.  ,  Ref.  10)  could  not  be  explained  by  extend¬ 
ing  the  results  discussed  above.  Rather,  a  study  of  the  flutter  cases  of 
interest  indicated  that  some  could  be  associated  with  mass -density  ratios 
as  much  as  an  order  of  magnitude  below  the  respective  A*  of  the  struc- 

3L 

ture.  Further,  all  of  the  related  numerical  analyses  showed  that  the 
structural  parametric  combinations  associated  with  the  low  M  flutter 
cases  exhibited  comparable  flutter  points  for  mass -density  ratios  ap¬ 
proaching  zero.  In  essence,  this  meant  that  the  calculated  /I  value 
(Eq.  1-28)  of  the  structure  did  not  represent  a  lower  mass -density  ratio 
limit  to  the  flutter  boundary  of  the  specific  lifting  surfaces  under  con¬ 
sideration.  Rather,  some  form  of  boundary  completely  predictable  by 
production  type  modal  analyses  incorporating  the  usual  two-dimensional, 
unsteady,  strip  theory  existed  in  the  low  mass -density  region  well  below 

the  /i  value  predicted  by  the  simplified  analysis.  The  portion  of  the 
a 

flutter  boundary  for  mass -density  ratios  less  than  the  /i  of  the  structure 
which  cannot  be  predicted  by  the  simplified  typical  section  analysis  marks 
precisely  the  region  of  interest  of  Part  II  of  this  study.  Hence,  it  is  here 
noted  that  the  term  "low  mass -density  ratio"  as  used  in  Part  II  is  meant 
to  imply  an  association  with  the  mass -density  ratio  region  below  the 
characteristic  M  of  the  system  in  question  and  for  which  the  simplified 

cL 

section  analysis  has  been  proven  to  be  partially  inadequate. 

A  review  of  the  well  documented  experimentally  and/or  theoreti¬ 
cally  observed  low  mass -density  ratio  flutter  occurrences  of  primary 
interest  quickly  revealed  that  all  of  the  structures  involved  were  thin, 
uniform,  swept,  cantilevered  lifting  surfaces  with  the  length  significantly 
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greater  than  the  width,  (minimum  l/b  ratio  of  8).  These  facts  had  been 
previously  noted.  In  particular,  the  probable  importance  of  sweep  had 
been  previously  stated  or  implied  by  the  author  and  several  other  in¬ 
vestigators.  (See  e.g.  ,  Ref.  5.)  Landahl,  et  al  (Ref.  26)  had,  in  fact, 
used  a  specific  example  to  show  the  probable  dependence  of  low  ju  flutter 
on  sweep.  However,  in  spite  of  these  facts,  the  mechanism  of  low  mass- 
density  ratio  flutter  remained  unknown.  Hence,  the  extrapolation  of  re¬ 
sults  from  the  specific  to  the  general  could  not  be  done  with  confidence. 

The  observation  that  certain  rather  special  structural  features 
might  be  an  important  part  of  the  low  fJ.  flutter  mechanism  certainly  in¬ 
dicates  that  these  features  must  be  retained  in  the  mathematical  model 
to  be  used  to  study  the  mechanism.  Conversely,  a  truly  basic  under¬ 
standing  of  the  flutter  phenomenon  requires  that  unnecessary  complexity 
be  avoided.  Just  what  features  could  be  expected  to  be  important  may  be 
reasoned  rather  simply  by  considering  the  inter-relationships  of  the 
structures  of  interest  and  the  hydrodynamic  forces  acting  upon  them. 

In  the  case  of  a  two-dimensional  swept  lifting  surface,  the  hydro - 
dynamic  forces  can,  in  general,  be  very  well  approximated  by  consider¬ 
ing  only  the  axis -normal  component  U  cos  A  °f  the  free  stream.  This 
leads  to  the  simple  result  that  if  a  two-dimensional  lifting  surface  is 
assumed  to  be  swept  by  rotating  the  surface  relative  to  the  free  stream, 
the  critical  flutter  speed  may  be  expected  to  vary  approximately  as 
the  secant  of  the  sweep  angle.  According  to  this  simplest  of  approxima¬ 
tions,  the  only  effect  of  sweep  upon  the  flutter  behavior  of  a  system  is 
to  increase  the  critical  flutter  speed  Thus,  it  is  clear  that  this  most 
basic  effect  of  sweep  cannot  be  responsible  for  the  occurrence  of  low 
mass -density  ratio  flutter.  Hence,  the  fact  that  the  structures  of 
interest  had  been  flexible  cantilevers  has  to  be  regarded  as  significant. 
Asa  result,  it  is  obvious  that  a  satisfactory  mathematical  model  must 
retain  the  salient  features  of  both  the  flexible  cantilever  rigidly  sup¬ 
ported  at  the  root  and  the  angle  of  sweepback  of  the  elastic  axis  without 
becoming  so  complex  that  cause  and  effect  relationships  are  lost. 
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Accordingly,  Part  II  of  this  investigation  is  based  upon  a  simple  binary 
representation  of  a  uniform  cantilevered  lifting  surface  which  includes 
the  provision  for  incorporating  an  angle  of  sweep  and  retains  a  spanwise 
variation  in  the  pitch  to  bending  ratio. 

II -B  The  Axis -Normal  and  Streamwise  Formulations 

The  problem  of  obtaining  a  simplified  mathematical  model  capable 
of  accurately  representing  an  oscillating,  uniform,  cantilevered  lifting 
surface  may  be  approached  in  two  separate  but  related  ways.  One  method 
may  be  considered  to  be  a  logical  refinement  of  the  simplest  effect  of 
sweep,  which  has  been  described  in  the  previous  section.  This  method 
arises  from  the  consideration  that  a  flexible  cantilevered  lifting  surface 
experiences  spanwise  deformations  and  that  in  consequence  the  spanwise 
component  of  the  free  stream  velocity  -  U  sin  A  •  cannot  be  ignored  with 
respect  to  the  hydrodynamic  forces  produced.  Hence,  in  the  case  of  the 
oscillating  cantilever,  the  spanwise  component  of  the  free  stream  also 
produces  downwash.  The  full  downwash  on  the  surface  is  thus  considered 
to  be  the  sum  of  the  two  downwash  components.  The  technique  has  then 
been  to  start  from  the  full  downwash  but  to  apply  to  it  only  the  axis -normal 
component  of  the  free  stream.  Consequently,  axis -normal  strips  form 
the  basis  of  this  analysis.  The  result  of  following  an  approach  of  this  type 
has  been  what  may  be  termed  a  first  order  approximation  to  the  swept 
cantilever  problem  which  has  become  known  as  the  "axis -normal"  method. 
This  method  was  developed  and  used  by  Barnby,  Cunningham,  and  Garrick 
(Ref.  19).  As  may  be  noted  from  Ref.  19,  the  refinements  upon  the  two- 
dimensional  analysis  introduced  by  this  approach  primarily  take  the  form 
of  additional  unsteady  hydrodynamic  terms,  the  importance  of  which  is 
largely  determined  by  the  "sweep  parameter"  )i  which  appears  as  a 
multiplying  factor  with  each  of  the  "new"  terms.  The  sweep  parameter 
}{  is  a  consequence  of  the  local  bending  slope  of  the  structure  and  hence 
incorporates  both  the  sweepback  and  the  cantilevered  aspects  of  the 
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structure.  It  may  be  mentioned  that  Herr  (Ref.  5)  used  the  method  of 
Ref,  1  9  in  his  calculations  and  was  thus  able  to  produce  some  of  the  first 
theoretical  results  indicative  of  the  possible  important  role  of  sweep  in 
low  mass -density  ratio  flutter. 

A  second  approach  has  been  to  consider  the  full  free  stream 
velocity  and  apply  it  to  streamwise  strips.  This  method,  which  has  be¬ 
come  known  as  the  "streamwise  approach11  must  lead  to  the  simplest 
approximation  as  the  swept  lifting  surface  approaches  an  infinite  aspect 
ratio,  just  as  the  "axis-normal"  approach  does.  Hence,  it  follows  that 
the  hydrodynamic  forces  to  be  used  on  a  streamwise  strip  which  is  part 
of  a  swept  structure  are  cos  ^  times  the  forces  that  would  be  calculated 
for  the  same  strip,  and  for  its  identical  motion,  in  the  unswept  case. 

The  details  of  combining  strip  theory  and  the  streamwise  approach  with 
a  cantilevered  lifting  surface  were  included  by  Jordan  in  Ref.  20.  As 
with  the  axis -normal  method,  the  result  of  following  the  approach  just 
described  is  an  approximation  to  the  swept  cantilever  problem  which  may 
again  be  termed  as  of  first  order.  However,  unlike  the  axis -normal 
approach,  the  refinements  do  not  introduce  any  unsteady  hydrodynamic 
derivatives  other  than  those  of  the  usual  unswept  or  ’’typical  section” 
formulations  discussed  in  Part  I.  Rather,  the  effects  of  sweep  are  intro¬ 
duced  structurally.  In  particular,  the  angle  of  sweep  of  the  elastic  axis 
results  in  a  skewness  between  the  elastic  axis  -  which  is  assumed  to  be 
straight  and  held  rigidly  at  the  root  -  and  the  streamwise  strip.  As  a 
result,  the  entire  streamwise  strip  is  not  a ssociated  with  one  point  on 
the  elastic  axis  but  rather  with  a  range  of  points  and,  therefore,  a  range 
of  bending  values.  Hence,  a  bending  of  the  elastic  axis  produces  both 
translation  and  pitch  of  the  streamwise  strip;  the  pitching  motion  being 
proportional  to  the  local  bending  slope  and  thus,  in  general,  a  function 
of  the  spanwise  location  of  the  strip  under  consideration.  In  addition, 
the  moment  of  the  hydrodynamic  forces  acting  on  the  strip  produces 
both  a  torsional  moment  and  a  bending  moment  on  the  elastic  axis.  The 
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details  of  the  equation  that  arises  from  the  streamwise  formulation  of  the 
problem  are  shown  in  Ref.  20  (see  also  Appendix  B). 

Obviously,  the  two  methods  are  not  entirely  similar,  particularly 
for  cantilevered  lifting  surfaces  of  modest  aspect  ratio.  The  streamwise 
approach  seems  to  have  an  advantage  in  that  the  entire  flow  is  considered. 
In  addition,  none  of  the  strips  end  at  the  side  boundaries  of  the  surface  as 
is  the  case  with  the  axis -normal  approach  at  the  tip  and  the  root.  A  dis¬ 
advantage  of  the  streamwise  approach  arises  if  the  structure  has  axis- 
normal  ribs,  or  if  the  structure  is  solid  but  an  appreciable  sweep  angle 
is  involved,  since  in  these  cases  a  camber  of  the  streamwise  strip  may 
result.  This  latter  difficulty  can  also  arise  with  the  axis -normal  approach 
if  rigid  streamwise  ribs  are  used.  To  examine  the  inherent  differences 
between  the  two  methods  any  more  carefully  is  quite  pointless  since  ex¬ 
perience  has  shown  that  both  methods  lead  to  practically  identical  flutter 
speeds  if  the  span  to  chord  ratio  of  the  lifting  surface  is  large  enough  to 
justify  the  strip  theory  involved  in  either  method  and  if  the  sweep  angle  is 
considered  to  be  moderate.  Indeed,  for  swept  surfaces  with  a  low  span 
to  chord  ratio,  there  is  little  doubt  but  what  the  use  of  strip  theory  and  the 
assumption  of  a  straight  elastic  axis  represent  a  degree  of  approximation 
considerably  greater  than  any  differences  resulting  from  the  choice  of 
either  the  axis -normal  or  the  streamwise  approach.  Thus  it  is  apparent 
that  either  approach  could  form  the  basis  for  this  investigation.  This 
latter  statement  is  further  substantiated  when  it  is  realized  that  the  flutter 
cases  which  form  the  background  for  this  investigation  are  characterized 
by  structures  having  relatively  high  l/b  ratios  and  low  sweep  angles  -  and 
hence  low  values  of  K.  However,  it  is  felt  that  the  approach  incorporating 
the  streamwise  strips  is  more  easily  visualized  and  hence  better  suited  to 
an  investigation  seeking  a  physical  interpretation  of  the  basic  low  mass- 
density  ratio  flutter  mechanism.  The  cause  and  effect  relationships  are 
easier  to  isolate  and  interpret  if  the  forms  of  the  unsteady  hydrodynamic 
terms  remain  unchanged  by  the  sweeping  process. 
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Furthermore,  the  streamwise  method  is  closely  related  to  the 
multi-panel  point  modal  analyses  commonly  used  in  industrial  production 
work  and  known,  by  virtue  of  experimental  evidence,  to  be  capable  of 
reasonably  accurate  low  mass -density  ratio  flutter  predictions.  Hence, 
the  streamwise  approach  forms  the  basis  of  this  investigation. 

II -C  The  Salient  Features  of  the  Streamwise  Approach  Applied  to  the 
Uniform  Cantilevered  Lifting  Surface 

The  combination  of  incompressible,  unsteady,  strip  theory  and  the 
streamwise  approach  as  applied  to  the  uniform  cantilevered  lifting  sur¬ 
face  forms  the  theoretical  basis  of  Part  II  of  this  investigation.  The  fact 
that  the  unsteady  hydrodynamic  theory  is  as  stated  implies  a  number  of 
assumptions,  such  as  potential  flow,  which  are  stated  in  virtually  all 
modern  flutter  texts  and  need  not  be  repeated  herein,  (see  e.  g.  ,  Ref.  14). 
In  addition,  assumptions  must  be  made  with  regard  to  the  structure  and 
how  the  swept  and  unswept  configurations  are  assumed  to  be  related.  For 
this  study,  the  unswept  reference  cantilevered  strut  and  the  swept  strut 
are  related  as  shown  in  Figure  ZZ  and  as  stated  below  after  the  fashion 
proposed  by  Jordan  in  Ref.  18. 

The  strut  in  Figure  ZZ  is  assumed  to  have  a  straight  elastic  axis 
which  is  rigidly  held  at  the  root.  As  the  transition  from  the  unswept  to 
the  swept  takes  place,  each  spanwise  strip  keeps  its  length  and  mass 
composition.  While  the  strip  rotates  about  its  root  point,  this  point 
slides  along  the  root  chord.  In  addition,  the  structure  is  assumed  to  be 
of  such  a  nature  that  the  streamwise  strips  remain  rigid.  By  sweeping 
in  this  manner,  the  halfchord  at  a  given  spanwise  station  of  the  swept 
strut  is  related  to  the  halfchord  of  the  straight  strut  counterpart  as  fol¬ 
lows  : 


b  =  b  sec  A 

A 


(II-l) 
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The  details  of  the  integro-differential  equation  that  is  the  result  of 
applying  strip  theory  to  the  streamwise  strip  of  a  cantilever  which  is 
swept  as  noted  above,  and  then  integrating  over  the  span  are  shown  in 
Ref.  20.  The  flutter  determinant  may  then  be  obtained  by  using  the 
Rayleigh -Ritz  method  to  approximate  the  vibratory  character  of  the  struc¬ 
ture  in  question.  Inasmuch  as  modal  analyses  using  as  few  as  two  modes 
have  successfully  predicted  the  existence  of  low  mass -density  ratio  flut¬ 
ter,  the  mathematical  model  chosen  for  this  investigation  is  the  binary 
flutter  determinant  of  a  uniform,  swept  strut;  the  assumed  flutter  modes 
being  the  fundamental  bending  and  torsional  modes  of  the  uniform  canti¬ 
lever.  The  salient  features  of  the  flutter  determinant,  based  on  the 
above  representation,  can  be  summed  up  briefly  as  follows.  Again,  the 
details  are  contained  in  Ref.  20  with  an  abbreviated  development  of  the 
determinant  included  herein  as  Appendix  B. 

As  noted  in  Part  I  and  Appendix  A,  the  flutter  determinant  of  the 
strictly  two-dimensional  typical  section  may  be  written  as: 
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®h  * 
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(II -2) 


where: 

B"“xa  4  '  V*  +  a) 

c  -  «  ^  -  x-h(i  +  ») 

D  =  |i  r0=  +  -  (X.^  +  +  a)  +  L^(J  +  a  |» 

(See  Appendix  A  for  definition  of  unsteady  derivatives  L^,  etc.  ) 
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Now  the  counterpart  of  Eq.  II-2  for  the  uniform  and  unswept 
cantilever  strut  may  be  considered.  The  flutter  determinant  for  this 
system  is  obtained  if  certain  spanwise  integrals  over  the  bending  and 
torsional  modes  are  introduced  in  Eq.  II-2.  For  convenience,  the  re¬ 
sultant  equations  may  be  normalized  to  yield: 
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M— > 
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(H-3) 


This  determinant  is  seen  to  be  identical  to  Eq.  (II-2)  except  for 
the  coefficients  of  the  off-diagonal  terms.  The  nearness  of  .  96  to  1 . 0 
has  in  the  past  been  regarded  as  evidence  that  the  two-dimensional 
typical  section  and  the  cantilever  flutter  problems  are  closely  related. 
Indeed,  it  may  be  stated  that  the  practical  usefulness  of  the  simplified 
typical  section  analysis  before  high  speed  digital  computers  became 
available,  was  largely  due  to  the  fact  that  it  could  be  used  to  adequately 
represent  a  cantilevered  aircraft  wing.  However,  it  will  be  shown  later 
in  this  study  that  the  close  relationship  between  the  two-dimensional  and 
the  flexible  cantilever  flutter  problems  does  not  hold  in  the  low  mass- 
density  ratio  region  due  to  an  extreme  sensitivity  of  the  flutter  behavior 
of  a  lifting  surface  to  its  vibratory  character  -  and  hence  to  its  assumed 
flutter  modes.  The  exact  value  of  the  off-diagonal  coefficients  depends 
upon  the  particular  choice  of  flutter  modes. 

Now  if  Eq.  (II-3)  is  regarded  as  the  flutter  determinant  of  an  un- 
swept  strut  and  the  corresponding  swept  strut  is  assumed  to  be  formed 
as  discussed  earlier  in  this  section,  the  flutter  determinant  of  the  swept, 
cantilever  strut  may  be  written  as: 
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(II -4) 

It  may  be  noted  that  the  changes  introduced  by  this  first  order 
approximation  of  the  swept  cantilever  are  entirely  determined  by  terms 
containing  the  sweep  parameter  K.  The  sweep  parameter  K  combines 
both  the  sweep  angle  A  and  the  length  to  semichord  ratio  of  the  equivalent 
or  "reference u  unswept  strut  as  follows; 

K  =  2  tan  h.jljb  (II-5 ) 

(See  also  the  List  of  Symbols  and  Figure  22.  ) 

Conversely,  it  may  be  seen  that  the  sweep  parameter  value  may 
be  changed  by  altering  any  one  of  the  three  structural  parameters  of 
sweep,  span,  or  semichord.  However,  it  must  be  recalled  that 
Eq.  (II-4)  is  based  upon  certain  assumptions  regarding  the  structure 
and  the  two-dimensional  unsteady  hydrodynamic  theory.  Clearly,  these 
assumptions  limit  the  classes  of  structures  which  may  be  regarded  as 
being  represented  by  Eq.  (H-4)  regardless  of  the  value  of  the  sweep 
parameter  K. 

It  may  be  seen  from  the  swept  strut  determinant  under  discussion 
that  if  K  is  assumed  to  be  zero,  as  is  the  case  if  the  angle  of  sweep  is 
zero  or  the  span  infinite,  Eq.  ( II  —4 )  reduces  to  Eq.  (U-3).  Further,  the 
A,  B,  etc.  of  Eq.  (II -4)  are  the  same  as  equations  (II - 3 )  and  (II-2).  Also, 
all  of  the  non-dimensional  parameters  -  including  /u  -  are  referred  to  the 
unswept  strut  and  thus  invariant  with  respect  to  the  sweep  parameter  K. 
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Hence,  the  effects  upon  the  flutter  behavior  of  a  structure  which  are 
entirely  due  to  sweep  -  or  more  correctly,  the  introduction  of  a  non¬ 
zero  sweep  parameter  -  are  easily  studied.  However,  the  flutter  speed 
of  the  swept  strut  is  now  related  to  the  reduced  frequency  k  by  the  re¬ 
lation. 


V 


cu.b  .  (D  b  sec  A 
f  A  f 


(II— 6 ) 


Hence,  it  may  be  seen  that  the  above  approximation  includes  the 
most  fundamental  effect  of  sweep. 

As  mentioned  in  the  last  section,  the  terms  which  include  the 
sweep  parameter  K  arise  from  the  pitch  of  the  streamwise  section  that 
is  coupled  to  the  local  slope  of  the  bending  mode  and  from  the  bending 
moment  acting  on  the  elastic  axis  which  is  a  component  of  the  stream- 
wise  moment.  The  fact  that  the  same  parameter  K  governs  both  effects 
might  be  expected  as  both  are  essentially  the  result  of  one  structural 
change.  Nevertheless,  it  is  fortunate  that  any  additional  complexity 
may  be  avoided. 

A  plot  of  the  sweep  parameter  K  for  various  configurations  is 
included  herein  as  Figure  23.  A  glance  will  suffice  to  show  that  for 
practical  structures  for  which  the  fundamental  assumptions  of  strip 
theory,  a  straight  elastic  axis,  and  a  rigid  streamwise  section  apply; 

K  is  normally  a  small  number.  Further,  it  may  be  noted  that  the 
majority  of  flutter  cases  -  theoretical  and  experimental  -  which,  in 
fact,  provided  the  motivation  for  this  investigation  (  see  e.g.  ,  Ref.  7 
or  10)  involved  extremely  low  values  of  K,  (on  the  order  of  1 0  1  to  10  s). 
Use  is  made  of  these  facts  to  simplify  the  swept  strut  flutter  equations. 

II -D  Simplification  of  the  Flutter  Determinant 

Before  considering  the  development  of  a  simplified  swept  strut 
equation,  a  few  remarks  regarding  the  reasons  for  seeking  a  simplified 
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version  are  in  order.  In  this  era  of  high  speed  digital  computers,  the 
flutter  determinant  represented  by  Eq.  (II-4),  could  obviously  be  easily 
handled  without  eliminating  any  of  the  terms.  Further,  it  might  be 
argued  that  any  trends  determined  as  a  result  of  a  parametric  study  of 
the  unsimplified  determinant  would  inherently  be  more  accurate  than 
what  might  be  expected  from  any  simplified  version.  However,  it  must 
be  remembered  that  the  purpose  of  this  study  is  not  to  obtain  numerical 
answers  to  specific  problems.  It  is  a  known  fact  that  reasonably  accu¬ 
rate  answers  may  be  obtained  by  production  type  model  analyses.  How¬ 
ever,  the  very  completeness  and  complexity  that  produces  the  degree  of 
accuracy  associated  with  the  production  analyses  makes  a  study  of  the 
general  cause  and  effect  relationships  difficult.  Rather,  an  insight  into 
the  fundamental  mechanism  of  low  mass -density  ratio  flutter  is  sought. 
To  this  end,  the  simplest  -  or  nearly  cimplest  -  mathematical  model 
which  retains  the  salient  features  of  the  mechanism  is  required.  For¬ 
tunately,  the  pitfall  of  oversimplification  presents  no  problem  as  the 
results  of  more  elaborate  techniques  -  and  experiments  -  are  available 
for  comparison.  Indeed,  it  is  the  very  fact  that  known  results  have 
shown  that  low  values  of  the  sweep  parameter  can  be  associated  with 
known  flutter  occurrences  that  suggests  that  a  simplified  version  of 
Eq.  (11-4)  may  be  used. 

If  Eq.  (II-4)  is  expanded,  certain  of  the  terms  may  be  combined 
to  produce  the  following  result: 
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It  may  be  seen  from  Eq.  ( II  —  7 )  that  the  coefficients  of  the  last  two 
terms  are  two  orders  of  magnitude  less  than  the  coefficients  of  the  other 
terms.  This  fact,  combined  with  the  fact  that  low  mass -density  ratio 
flutter  has  been  predicted  and  experimentally  verified  for  struts  having 

—  1  —  3 

sweep  parameter  values  on  the  order  of  10  to  10  ,  serves  to  strongly 

indicate  that  the  last  two  terms  may  be  entirely  eliminated  without  the 
character  of  the  problem  being  lost.  Similar  reasoning  leads  to  the  use 
of  1. 0  as  the  coefficient  of  the  fourth  term.  In  addition  to  strictly  numeri¬ 
cal  reasons,  a  physical  justification  to  the  above  stated  simplification  is 
seen  to  exist  when  it  is  realized  that  the  exact  values  of  the  numerical 
coefficients  in  Eq.  (II  —  7 )  are  functions  of  the  assumed  flutter  modes. 

Thus,  the  terms  neglected  can  be  viewed  as  the  results  of  the  effects  of 
sweep  upon  the  mode  shapes  and  the  hydrodynamic  force  changes  result¬ 
ing  therefrom.  Clearly,  these  effects  could  well  be  expected  to  be 
negligible  for  small  angles  of  sweep.  As  will  be  shown  in  the  next  sec¬ 
tion  of  this  investigation,  the  above  reasoning  is  fully  justified  by  the 
accuracy  of  the  results  obtained  from  the  simplified  equation. 

By  making  the  simplifications  discussed  above,  the  expansion  of 
the  flutter  determinant  may  be  written  as  follows: 

a)  8  a)  3 

(A  -  n  n  — )(D  -Mr3  —  )  -  .  96*  BC 
8  O'  a 


(B  +  C+^D)  =  0 


(II-8a) 


In  effect,  the  above  is  similar  to  Eq.  (II -3 ),  but  with  a  term  due 
to  sweep  added  and  thus  is  seen  to  contain  all  of  the  essential  features 
of  the  swept  cantilever. 
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It  may  be  seen  from  Eq.  (II-8a)  that  the  effect  of  the  spanwise 
integrals  over  the  assumed  mode  shapes  appears  only  as  the  coefficient 
of  the  second  term  -  or  more  explicitly  as  the  coupling  coefficient  of  the 
off-diagonal  terms  of  the  determinant  of  the  unswept  cantilever.  Since 
production  type  modal  analyses  have  shown  a  strong  sensitivity  to  the 
number  and  character  of  the  assumed  flutter  modes,  it  is  desired  to  set 
up  the  equations  so  that  the  effect  of  varying  the  off-diagonal  coupling 
constant  may  be  studied.  Accordingly,  the  second  term  of  Eq.  (II-8a)  is 
expressed  as  a  sum  and  difference  as  follows: 
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(II -8b) 


By  substituting  the  appropriate  unsteady  aerodynamic  terms  in 
the  above  equation  as  indicated  by  the  definitions  following  Eq.  (II-2)  and 
Appendix  A,  and  separating  the  resultant  complex  equation  into  its  real 
and  imaginary  parts,  the  equations  shown  on  the  next  page  may  be  ob¬ 
tained,  (Eqs .  II-9a  and  9b). 

As  in  Part  I,  the  following  substitutions  are  made: 
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,  the  flutter  frequency  coefficient  squared;  and 
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X  =  ( - )  the  flutter  speed  coefficient  squared. 
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Equations  (II-9a)  and  (II-9b) 
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In  both  of  these  equations,  the  terms  introduced  by  the  inclusion 
of  the  sweep  parameter  are  readily  identifiable.  The  terms  which  ex¬ 
press  the  difference  of  the  off-diagonal  coupling  constants  from  unity  are 
in  the  numerator  in  both  cases,  and  are  preceded  by  the  decimal  values. 
As  with  the  equations  obtained  for  the  typical  section  and  discussed  in 
Part  I,  Eqs.  (H-9a)  and  (II-9b)  are  transcendental  and  require  a  trial 
and  error  solution.  Some  results  of  solving  Eqs.  (II-9a)  and  (II  —  9b )  are 
discussed  in  the  following  section. 

II -E  Some  Typical  Results  Based  Upon  the  Simplified  Theory 

Since  previous  work  had  shown  that  production  type  modal  analyses 
are  capable  of  predicting  the  low  mass -density  ratio  flutter  of  uniform, 
cantilever  struts,  and  further  since  at  least  part  of  the  prior  theoretical 
work  had  been  experimentally  verified,  results  are  available  which  may 
be  used  to  check  the  validity  of  the  simplified  approach  used  to  obtain 
Eqs.  (II~9a)  and  (H-9b).  Accordingly,  several  of  the  models  for  which 
results  exist  have  been  investigated  using  the  simplified  approach  out¬ 
lined  in  the  previous  two  sections.  In  addition  to  checking  the  simplified 
flutter  equation,  these  investigations  show  the  character  of  the  flutter 
solutions  as  a  function  of  the  two  parameters  of  primary  interest,  namely 
the  sweep  parameter  K  and  the  mass-density  ratio  fJ.  .  Specifically,  the 
structural  non-dimensional  parametric  values  used  for  a  major  portion 
of  the  preliminary  numerical  work  outlined  herein  are  those  of  Model  #2 
used  by  Baird,  et  al,  (Ref.  10)  and  Model  #2  used  by  Herr  (Ref.  5).  In 
addition,  certain  of  the  conclusions  of  this  investigation,  which  will  be 
discussed  in  the  sections  to  follow,  are  compared  to  some  as  yet  un¬ 
published  results  which  have  been  obtained  by  C.  E.  Squires,  Jr.  ,  of 
Grumman  Aircraft  Engineering  Corporation  for  a  third  model.  These 
models  are  referred  to  herein  as  Baird's  Model  #2,  Herr’s  Model  #2, 
and  Squires'  Model,  respectively.  That  Eq.  (H-8b)  represents  a 
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satisfactory  mathematical  representation  of  the  low  mass -density  flutter 
problem  was  verified  when  the  resultant  solutions  were  found  to  agree 
quite  well  with  the  results  of  the  more  elaborate  techniques  used  by  the 
above  named  investigators.  It  naturally  follows  that  the  agreement  with 
experimental  results  was  also  found  to  be  acceptably  good. 

Figure  24  shows  a  comparison  of  the  results  obtained  by  Herr 

(Ref.  5),  using  the  theory  of  Ref.  19,  with  results  obtained  using  the 

simplified  analysis  for  the  unswept  Herr  #2  model.  In  both  cases,  no 

low  mass -density  ratio  flutter  (/i  <  )  is  predicted,  while  a  range  of 

a 

flutter  speeds  is  predicted  for  the  higher  mass -density  ratios.  It  can  be 
seen  that  the  agreement  is  quite  good,  particularly  with  respect  to  the 
lower  boundary  which  is  the  one  of  greatest  practical  interest  Since  the 
two  underlying  unsteady  hydrodynamic  theories  approach  each  other  as 
the  sweep  angle  approaches  zero,  it  is  very  likely  that  the  slight  dif¬ 
ferences  shown  in  Figure  24  are  due  to  1)  the  slightly  different  off- 
diagonal  coupling  constants  used,  and  2)  the  fact  that  the  frequencies  used 
in  the  simplified  analysis  were  not  the  fundamental  bending  and  torsional 
frequencies  of  the  unswept  strut  but  rather  the  swept  measured  frequen¬ 
cies  of  the  first  and  third  coupled  modes,  which  were  predominantly 
fundamental  bending  and  fundamental  torsion  respectively.  It  might  be 
mentioned  that  Figure  24  is  typical  in  that  no  finite  low  mass -density 
ratio  solutions  could  be  obtained  with  a  sweep  parameter,  K,  of  zero  for 
any  of  the  many  combinations  of  system  parameters  investigated.  This 
was  in  keeping  with  what  other  investigators  had  suggested,  (see  e.g.  , 
Ref.  26)  and  was  a  factor  which  contributed  to  the  formulation  of  the 
physical  interpretation  of  the  low  mass -density  ratio  flutter  mechanism 
which  is  presented  in  a  section  to  follow. 

Figure  25  shows  that  both  theories  predict  a  low  mass -density 
ratio  flutter  boundary  for  the  Herr  Model  #2  with  a  10  degree  sweep 
angle.  As  may  be  seen,  the  agreement  is  quite  good.  The  character 
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of  the  low  mass -density  ratio  solutions  shown  in  Figure  25  may  be  re¬ 
garded  as  typical  for  parametric  combinations  possessing  such  solutions. 
In  particular,  Figure  25  illustrates  the  local  maximum  of  the  flutter  speed 
that  occurs  for  some  mass -density  ratios  near  the  fi  of  the  system.  This 

cL 

local  maximum  may  be  said  to  divide  the  MusualM  from  the  low  mass- 
density  ratio  flutter  boundaries,  (see  also  Ref.  11).  It  might  be  mentioned 
that  Herr  tested  this  model  with  a  sweep  angle  of  15  degrees  and  a  mass- 
density  ratio  jU  of  .  13  and  did  not  obtain  flutter.  However,  the  maximum 
speed  obtained  during  the  test  was  less  than  one  half  of  the  predicted 
flutter  speed  and  thus  the  absence  of  experimental  verification  is  under¬ 
standable.  It  will  be  shown  in  a  section  to  follow  that  the  comparatively 
high  flutter  speed  predicted  and  quite  possibly  physically  existent  flutter 
speed  was  at  least  partially  due  to  the  unfortunate  choice  of  a  forward 
elastic  axis  location. 

Figure  26  shows  a  comparison  of  solutions  for  the  Baird  Model  #2 
with  15  degrees  of  sweep.  The  solutions  obtained  by  Squires  and  Baird 
were  obtained  by  using  production  type  modal  analyses  incorporating  the 
lowest  three,  four,  or  five  coupled  modes  of  the  strut.  Again,  the  agree¬ 
ment  was  quite  good,  particularly  with  regard  to  the  low  mass -density 
ratio  quartic  -  or  four  mode  -  analysis.  As  may  be  seen,  the  agreement 
with  the  two  experimental  points  was  excellent.  In  addition,  Figure  26 
demonstrates  the  "unusual"  sensitivity  of  the  low  mass -density  ratio 
modal  solutions  to  the  modal  representation  of  the  structure.  The  sensi¬ 
tivity  has  been  termed  "unusual”  because  of  the  comparatively  weak  de¬ 
pendence  of  the  "usual"  (fJ.  >  fJ.  )  flutter  boundary  upon  the  assumed  flutter 

3, 

modes.  This  characteristic  of  low  mass -density  ratio  flutter  will  be  dis¬ 
cussed  in  considerable  detail  in  the  sections  to  follow.  It  will  be  shown 
that  such  a  behavior  can  be  readily  explained. 

In  addition  to  the  flutter  speed  or  flutter  speed  coefficient  compari¬ 
sons  discussed  above,  the  flutter  frequency  solutions  were  also  checked 
at  the  discrete  points  for  which  comparable  solutions  were  available. 

Aga  in,  the  agreement  was  good. 
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The  results  which  have  been  discussed  represent  a  typical  samp¬ 
ling  of  the  calculations  performed  to  assure  the  validity  of  Eq.  (II -8b ) 
and  have  been  presented  to  illustrate  why  it  was  concluded  that  the  simpli 
fied  equation  could  be  used  to  study  the  fundamental  mechanism  of  low 
mass -density  ratio  flutter.  Having  established  the  validity  of  the  simpli¬ 
fied  approach,  certain  of  the  specific  characteristics  of  the  solutions  re¬ 
sulting  from  the  use  of  the  approach  will  be  discussed. 

Figure  27  is  a  typical  plot  of  the  flutter  speed  coefficient  as  a 
function  of  the  mass -density  ratio  for  several  sweep  parameter  values. 
This  plot  shows  that  although  the  physically  meaningful  flutter  boundary 
is  only  slightly  affected  by  small  changes  in  the  value  of  the  sweep  para¬ 
meter  when  mass -density  ratios  appreciably  greater  than  p  are  in- 

a 

volved,  (the  \x  of  the  system  depicted  is  .76),  a  radical  change  in  the 
a 

low  mass -density  ratio  behavior  occurs  as  the  sweep  parameter  is  in¬ 
creased  beyond  what  might  be  termed  the  "critical*1  K  of  the  system  in 
that  a  finite  valued  flutter  boundary  is  seen  to  appear.  As  the  sweep 
parameter  is  further  increased  -  at  least  within  the  physically  meaning¬ 
ful  range  -  the  low  mass -density  ratio  boundary  is  seen  to  remain  es¬ 
sentially  unchanged.  Thus  the  flutter  speed,  for  configurations  possess¬ 
ing  a  sweep  parameter  greater  than  the  critical  value,  is  predicted  to 
rise  approximately  as  a  function  of  the  secant  of  the  sweep  angle.  The 
results  of  more  elaborate  analytical  techniques  and,  in  fact,  experi¬ 
mental  results  have  shown  the  rise  in  the  flutter  speed  as  a  function  of 
sweep  to  be  somewhat  greater.  However,  it  must  be  recalled  that  both 
the  simplified  and  the  complete  binary  flutter  determinants  involve  as¬ 
sumptions  which  deteriorate  in  validity  as  the  angle  of  sweep  is  increased 
These  assumptions  have  been  discussed  in  previous  sections. 

Another  feature  of  the  low  mass -density  ratio  solutions  which  is 
typified  by  Figure  27  is  that  if  the  mass -density  ratio  can  be  assumed  to 
be  varied  by  varying  the  fluid  density,  the  low  mass -density  ratio 
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boundary  is  seen  to  be  a  constant  dynamic  pressure  phenomenon.  It 
should  be  noted  that  certain  of  these  characteristics  have  been  noted  by 
other  investigators  for  this  and  other  models,  (see  e.g.,  Ref.  12). 

Although  Figure  Z7  may  be  regarded  as  typical  for  those  para¬ 
metric  combinations  which  possess  a  valid  low  mass -density  ratio  boun¬ 
dary,  many  combinations  do  not  possess  low  mass -density  ratio  solu¬ 
tions  (and  thus  might  be  presumed  to  be  flutter  free  in  the  low  p  region) 
for  any  reasonable  sweep  parameter  values.  Qualitative  and  quantitative 
reasons  for  this  behavior  are  discussed  in  sections  to  follow.  In  all 

cases,  the  effect  of  varying  the  sweep  parameter  upon  the  nusualM  (p  >  p  ) 

a 

flutter  boundary  will  be  shown  to  be  slight. 

Figure  28  is  a  magnification  of  the  low  p  portion  of  Figure  27,  and 
shows  the  behavior  as  the  zero  mass -density  ratio  is  approached.  In  ad¬ 
dition,  the  numerical  values  give  some  indication  of  the  range  of  reduced 
frequencies  (k)  typically  encountered  with  the  low  mass -density  ratio 
flutter  solutions.  It  may  be  noted  that  for  any  given  low  mass -density 
ratio,  the  reduced  frequency  rises  as  the  sweep  parameter  is  increased. 

Figure  29  is  the  flutter  frequency  counterpart  of  Figure  27  and 
may  again  be  regarded  as  typical.  As  can  be  seen,  the  flutter  frequency 
for  those  parametric  combinations  possessing  a  low  p  boundary  ap¬ 
proaches  zero  as  the  mass -density  ratio  approaches  zero.  Further,  the 
flutter  frequency  at  any  given  low  mass-density  ratio  decreases  as  the 
sweep  parameter  decreases.  Figure  30  is  the  flutter  frequency  coefficient 
counterpart  of  Figure  28. 

With  regard  to  the  characteristics  discussed  above,  by  far  the 
most  interesting  fact  is  that  the  flutter  frequency  due  to  any  given  para¬ 
metric  combination  is  seen  to  approach  zero  as  the  sweep  parameter  ap¬ 
proaches  the  ’'critical"  value.  Hence,  the  flutter  is  seen  to  degenerate  to 
a  pure  divergence.  This  behavior  is  illustrated  by  Figure  31  which  has 
been  computed  for  the  Baird  #2  model  at  a  mass -density  ratio  of  .  1 .  As 
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may  be  seen  from  the  figure,  the  "critical11  sweep  parameter  value  is 
approximately  .  03.  This  corresponds  to  approximately  6.8  degrees  of 
sweep  which  agrees  very  well  with  the  statement  in  Ref.  10,  that,  "a 
sweep  of  7. 5  degrees  was  sufficient  to  assure  freedom  from  static  diver¬ 
gence."  Figure  32  is  the  corresponding  flutter  speed  coefficient  plot 
which  is  seen  to  be  a  very  weak  function  of  the  sweep  parameter  for  values 
above  the  "critical"  value. 

In  an  effort  to  obtain  some  idea  of  what  value  the  approximate 
quasi -steady  (k  -♦  0)  and  the  quasi -unsteady  (k  -♦<»)  hydrodynamic  theories 
might  be,  a  few  sample  calculations  were  performed.  The  procedure 
followed  .was  simply  to  utilize  Eqs.  (H-9a)  and  (II-9b)  with  the  appropri¬ 
ate  F  and  G  values  inserted.  (The  quasi-steady  and  unsteady  approxima¬ 
tions  were  discussed  in  considerable  detail  in  Part  I.  )  Figure  33  shows 
a  comparison  of  the  flutter  speed  coefficients  predicted  by  the  "exact"  and 
by  both  approximate  theories  for  the  Baird  #2  model  at  zero  degrees  of 
sweep.  As  can  be  seen,  the  qua  si -unsteady  theory  (k  «)  agrees  reason¬ 
ably  well  with  the  "exact"  theory  while  the  quasi-steady  theory  (k  -♦  0) 
actually  predicts  a  non-existent  set  of  flutter  points  in  the  low  mass- 
density  ratio  region.  It  may  be  seen  that  the  reduced  frequencies  asso¬ 
ciated  with  the  low  mass -density  ratio  solutions  suggest  that  the  quasi- 
steady  assumption  might  not  be  valid.  Figure  34  is  the  flutter  frequency 
counterpart  of  Figure  33. 

Figure  35  shows  the  same  comparison  but  for  a  sweep  parameter 
value  greater  than  the  "critical"  value  of  the  system.  Once  again,  the 
quasi -unsteady  theory  (k  -♦  m)  is  seen  to  be  qualitatively  correct  over  the 
entire  mass -density  ratio  range  while  the  quasi -steady  theory  is  valid 
only  in  the  low  mass -density  ratio  range  and  even  there  only  if  careful 
attention  is  paid  to  the  reduced  frequencies  associated  with  the  predicted 
flutter  points.  It  might  be  mentioned  that  the  quasi-steady  boundaries 
form  a  loop  which  meets  at  the  origin.  Because  of  the  typically  low  re¬ 
duced  frequencies  normally  associated  with  low /i  flutter,  the  quasi- 
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steady  branch  associated  with  the  low  reduced  frequencies  is  very  nearly 
exact.  Figure  36  is  the  flutter  frequency  counterpart. 

Because  of  the  limited  number  of  calculations  performed  utilizing 
the  approximate  theories,  it  cannot  be  stated  with  certainty  that  the  fig¬ 
ures  discussed  above  are  typical  of  all  parametric  combinations.  In  fact, 
it  is  quite  likely  that  the  inconsistencies  encountered  in  conjunction  with 
the  quasi-steady  theory  are  at  least  partially  due  to  the  combination  of  an 
aft  center  of  gravity  and  a  very  low  frequency  ratio.  Indeed,  analyses 
discussed  in  Part  I  of  this  study  showed  that  such  a  behavior  could  be  ex¬ 
pected  with  a  Mtypical  section'1  having  such  properties.  Nevertheless,  it 
is  interesting  to  note  that  once  again  the  quasi -unsteady  theory  (k  -►  oo)  was 
found  to  be  qualitatively  reliable  while  the  quasi -steady  theory  could  only 
be  used  with  considerable  caution. 

When  all  of  the  results  which  have  been  discussed  thus  far  are 
viewed  as  a  whole,  one  fact  quickly  emerges.  Namely,  the  rather  radical 
changes  in  the  flutter  behavior  which  occur  as  a  function  of  the  mass- 
density  ratio  and  the  sweep  parameter  suggest  that  low  mass -density 
ratio  flutter  does  not  have  the  same  physical  origin  as  the  "usual" 
coupled  bending -torsion  flutter  familiar  to  aircraft  flutter  analysts. 

II -F  Investigation  of  the  Sensitivity  of  Solutions  to  the  Assumed 
Flutter  Modes 

The  simplified  two -degree -of -freedom,  uniform  cantilever,  flutter 
equation  (Eq.  H-8a)  was  developed  by  ignoring  the  effect  of  sweep  upon 
the  assumed  flutter  mode  shapes.  As  a  result,  only  the  direct  effect  of 
the  spanwise  integrals  over  the  unswept  mode  shapes  appears  in  Eq.  ( II  —  8a ) ; 
this  being  the  value  taken  by  the  numerical  coefficient  of  the  second  group 
of  terms.  A  glance  at  Eq,  (U-8a)  reveals  that  if  the  numerical  coefficient 
of  (.  96)  is  replaced  by  unity,  the  first  two  terms  of  the  equation  are 
identically  the  mathematical  model  of  the  familiar  "typical  section.  M  The 
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equation  which  is  the  result  of  altering  Eq.  (II-8a)  by  replacing  (.  96)  by 
1. 0  may  then  be  viewed  as  the  representation  of  a  simple  pitch-plunge 
section  with  the  major  sweep  effect  of  pitch-bending  added. 

The  altered  equation  discussed  above  was  investigated;  the  re¬ 
sults  being  that  no  low  mass -density  ratio  solutions  could  be  obtained  re¬ 
gardless  of  the  sweep  parameter  assumed  although  the  high  mass -density 
ratio  results  were  virtually  unaltered  from  what  could  be  obtained  by  using 
Eq.  (II-8a).  The  nature  of  the  solutions  is  shown  by  Figure  37,  which  may 
again  be  regarded  as  typical.  As  can  be  seen  from  Figure  37,  the  only 
effect  of  increasing  the  sweep  parameter,  from  a  low  mass -density  ratio 

point  of  view,  was  to  lower  the  of  the  system  an  almost  insignificant 

a 

amount.  The  obvious  implication  of  the  nature  of  the  results  depicted  by 

Figure  3  7  is  that  the  pitch-plunge  "typical  section"  analysis,  even  with 

the  sweep  effect  of  pitch-bending  included,  is  incapable  of  predicting  the 

low  mass -density  ratio  {fx  <  fx  )  flutter  behavior  of  a  swept,  flexible 

a 

cantilevered  strut.  This  conclusion  is  not  entirely  unexpected  as  it  has 
been  known  that  single  panel  point  modal  analyses  in  which  the  sweep 
condition  has  been  simulated  by  introducing  a  streamwise  angle  of  attack 
into  the  bending  mode  have  also  failed  to  predict  the  existence  of  low 
mass -density  ratio  flutter  for  configurations  known  to  be  susceptible  to 
the  phenomenon.  Note  that  this  result  implies  that  experimental  tests 
conducted  by  using  semi-rigid  models  are  of  questionable  value  in  pre¬ 
dicting  the  low  mass -density  ratio  flutter  behavior  of  flexible  structures. 

Figure  38  shows  a  typical  set  of  flutter  boundaries  obtained  by 
varying  the  off-diagonal  coupling  constant  of  Eq.  (II-8a)  from  .9*  to  1.0. 

It  may  be  seen  that  although  the  high  jx  flutter  boundary  remains  un¬ 
changed,  the  low  M  flutter  boundary  is  extremely  sensitive  to  the  coup¬ 
ling  constant  utilized.  Quite  obviously,  the  low  fx  boundary  trend  is  such 
that  an  infinite  flutter  speed  coefficient  is  rapidly  approached  as  the 
coupling  constant  approaches  unity.  Hence,  a  low  mass-density  ratio 
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(ju  <  fi  )  flutter  boundary  is  non-existent  for  a  coupling  of  1.0.  When  it  is 
a 

recalled  that  the  value  of  the  coupling  constant  arises  from  spanwise  inte¬ 
grals  over  the  assumed  mode  shapes,  it  is  clear  that  low  flutter  is  ex¬ 
tremely  “mode  shape  sensitive.  11  Indeed,  Figure  26  illustrates  that  the 
results  of  production  type  model  analyses  have  exhibited  a  strong  sensi¬ 
tivity  to  the  number  and  type  of  assumed  modes  used  in  the  flutter  analy¬ 
sis.  From  the  standpoint  of  the  production  analyses,  the  extreme  sensi¬ 
tivity  of  the  solutions  to  the  coupling  constant  used  implies  that  not  only 
must  a  considerable  number  of  modes  be  assumed,  but  the  number  of 
panel  points  must  be  sufficient  to  insure  accurate  and  representative 
spanwise  numerical  integration.  In  short,  the  vibratory  character  of  the 
structure  must  be  representative  to  a  degree  far  in  excess  of  what  has 
come  to  be  regarded  as  acceptable  for  aircraft  work.  The  problems  in¬ 
herent  with  this  extreme  sensitivity  when  complex  structures  are  in¬ 
volved  are  obvious. 

II -G  The  Zero  Mass -Density  Ratio  Analysis 

All  discussion  thus  far  in  Part  II  has  been  based  upon  the  physical 
assumption  that  a  mass -density  ratio  variation  is  achieved  by  considering 
a  variable  fluid  density.  Indeed,  the  statement  that  the  low  mass -density 
ratio  solutions  exhibited  a  constant  critical  dynamic  pressure  implied  that 
such  must  be  the  case.  If  the  above  reasoning  is  extended  to  the  limiting 
case,  it  is  clear  that  a  zero  mass -density  ratio  implies  a  fluid  of  infinite 
density.  As  one  might  intuitively  expect,  the  flutter  behavior  as  the  zero 
mass -density  case  is  assumed  to  be  approached  in  the  above  stated  man¬ 
ner  is  for  the  flutter  speed  and  flutter  frequency  to  both  approach  zero. 
This  results  in  an  indeterminate  solution  for  the  reduced  frequency. 

These  tendencies  are  both  very  well  illustrated  by  Figures  28  and  29 
which  are,  of  course,  based  upon  parametric  studies  which  assumed  con¬ 
stant  structural  mass  and  stiffness  properties  throughout  the  mass -density 
ratio  range. 
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However,  a  zero  mass -density  ratio  can  also  be  approached  by  de¬ 
creasing  the  section  mass  while  holding  the  fluid  density  constant.  The 
details  of  these  two  physically  different  approaches  were  discussed  in 
some  detail  in  Part  I  of  this  investigation  and  the  advantages  of  each  were 
pointed  out  at  that  time.  Briefly,  a  meaningful  investigation  of  the  limit¬ 
ing  zero  mass -density  case  is  possible  by  writing  the  flutter  equations  in 
terms  of  "in-fluid11  frequencies  and  then  using  the  concept  of  a  "decreas¬ 
ing  section  mass.  "  The  physical  interpretation  of  finite  "in  fluid"  fre¬ 
quency  ratios  coupled  with  a  zero  section  mass  is  quite  simply  that  the 
strut  is  assumed  to  retain  its  stiffness  properties  as  the  section  mass 
approaches  zero. 

Using  the  relationships  shown  herein  as  Eqs.  (1-7),  Eq.  (H-9a) 
may  be  converted  to  its  "in  fluid"  equivalent  which  may  be  written  as 
shown  on  the  following  page  as  Eq.  (Il-lOa).  Similarly,  Eq.  (U-9b)  may 
be  written  as  Eq.  (II-10b). 

It  is  obvious  from  the  simple  substitutions  which  result  in 
Eqs.  (II -10a)  and  (II-10b)  that  the  two  sets  of  equations  are  mathematically 
equivalent.  However,  a  glance  at  Eqs.  (II-9a)  and  (II-9b)  will  suffice  to 
show  that  the  equations  become  worthless,  in  a  practical  sense,  as  a  zero 
mass -density  ratio  is  approached  inasmuch  as  the  reduced  frequency  k  is 
indeterminate  regardless  of  the  choice  of  structural  parameters.  How¬ 
ever,  Eqs.  (II-10a)  and  (II-10b)  may  yield  finite  and  determinate  reduced 
frequency  k  solutions  if  the  zero  mass -density  ratio  is  considered  to  be 
approached  by  decreasing  the  section  mass  while  retaining  the  stiffness 
properties;  i.e.,  ^  while  II  *  remains  real  and  finite. 

Assuming  the  limiting  case  of  a  zero  mass -density  ratio  to  be 
approached  in  the  above  manner  results  in  the  Eqs.  (II-lla)  and  (U-llb) 
which  are  shown  on  page  111. 


k[2  - F  +n  2(1  -F -2a+4a®F+8aFx0 -2xa+4Frtta)] +4G[l  -n  (1  +2a  +2^  )] 
3136{/ix0(4aFk  -  k  -2G)+|-(1  +  2a)(Ga  - F3  )  +  2FG(  1  -4aa) 

+  G(1  +  4a)  +  ak(l  -  4aF)j  _ _ _ 

8Fk<?+a"  +/i  ^■8)+kn*<1  +*i  K2  ~2F  -4a  +  8a*F)  -4G(2a  +  1)(1  +fi  )n  a 
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Equations  (II-10a)  and  (II-10b) 
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Hence,  a  simple  mathematical  interpretation  of  what  has  been 
shown  is  that  the  "in  vacuum"  approach  causes  all  solutions  to  be  col¬ 
lapsed  onto  the  origin  while  the  Min  fluid”  approach  stretches  the  solu¬ 
tions  from  a  point  to  a  line.  From  a  practical  standpoint,  the  latter  ap¬ 
proach  offers  two  significant  advantages.  First,  all  flutter  frequencies 
do  not  approach  zero  as  the  mass-density  ratio  approaches  zero.  Rather, 
even  for  the  limiting  zero  mass -density  case,  only  the  "critical"  sweep 
parameter  values  yield  zero  flutter  frequency  solutions.  (This  point  is 
illustrated  by  Figure  39;  Figure  40  being  the  flutter  speed  counterpart.  ) 
As  may  be  seen  from  Figure  39,  real  finite  flutter  frequencies  may  be 
calculated  for  sweep  parameter  values  greater  than  the  critical  value 
while  no  solutions  exist  for  lesser  values.  Secondly,  the  structural  non- 

dimensional  parameters  IIs,  r  a,  and  x  do  not  appear  independently  in 

c  m 

the  zero  mass -density  equations,  but  rather  only  appear  in  combined 
form  as  the  "in  fluid”  frequency  ratio  11^.  (To  avoid  any  misunderstand¬ 
ing  it  is  recalled  that  11^  is  a  form  of  stiffness  ratio.  )  For  the  limiting 
case  of  interest,  the  relationship  may  be  written  as: 
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(11-12) 


As  was  shown  in  Part  I,  this  reduces  the  task  of  conducting  a 
parametric  study  by  several  orders  of  magnitude.  In  particular,  the 
above  stated  fact  allows  the  complete  zero  mass-density  ratio  solution 
pattern  for  uniform  cantilever  struts  to  be  discerned  by  considering  only 
the  three  parameters  11^,  a>  an<^  This  fact,  coupled  with  the  use  of 
high  speed  digital  computers,  allows  a  complete  parametric  study  cover¬ 
ing  the  physically  useful  ranges  of  parameters  to  be  conducted  quite 
easily  and  efficiently.  The  result  is  a  compact  representation  of  what 
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parametric  combinations  possess  tero  mass -density  ratio  solutions  and, 

of  particular  interest,  the  minimum  sweep  parameter  value  required  for 

a  valid  flutter  solution.  This  information  is  included  herein  as  Figure  41, 

which  shows  the  critical  sweep  parameter  value  as  a  function  of  the  elastic 

axis  location  for  various  values  of  II  *  . 

* 

If  Figure  41  were  only  valid  for  the  limiting  zero  mass-density 
ratio  case,  it  would  only  be  of  academic  interest.  However,  it  may  be 
recalled  that  the  zero  mass -density  ratio  case  is  merely  the  extreme 
limit  of  the  low  mass -density  ratio  region  and  as  such,  any  physical  ef¬ 
fects  due  to  the  decreasing  of  the  mass -density  ratio  are  retained  and, 
in  fact,  magnified  as  the  limiting  condition  is  approached.  Indeed,  it  is 
intuitively  obvious  and  is,  in  fact,  illustrated  by  many  of  the  figures  in¬ 
cluded  herein  that  the  low  mass-density  ratio  flutter  boundaries  (fJ.  <  ) 

a 

terminate  at  the  zero  mass -density  ratio.  (Literally  hundreds  of  solu¬ 
tions  by  the  author  and  others,  in  addition  to  what  has  been  included  here¬ 
in,  furnish  the  numerical  proof  that  the  above  statements  are  correct.  ) 
Hence,  the  non-existence  of  a  zero  mass -density  ratio  flutter  solution 
implies  the  non-existence  of  low  mass -density  ratio  flutter  for  the  con¬ 
figuration  in  question.  Therefore,  Figure  41  is  seen  to  be  a  completely 
practical  plot  of  which  configurations  may  be  expected  to  have  a  low 
mass -density  ratio  flutter  problem  for  some  real,  finite,  free  stream 
velocity;  and  further,  the  minimum  sweep  parameter  value  required  for 
a  dynamic  instability  to  develop.  At  this  point  it  is  again  emphasized 
that  Figure  41  is  based  upon  the  assumptions  of  a  uniform,  cantilevered, 
fully  wetted,  lifting  surface  and  further  that  all  the  assumptions  in¬ 
volved  in  two-dimensional,  unsteady,  incompressible,  potential  flow 
are  inherent  in  the  development.  Further,  the  simplified  equation  used 
to  obtain  Figure  41  inherently  decreases  in  accuracy  as  the  "critical11 
sweep  parameter  increases.  However,  it  may  be  noted  that  Figure  41 
checked  very  well  with  all  applicable  cases  found  in  the  literature,  in¬ 
cluding  the  three  models  specifically  mentioned  in  an  earlier  section. 
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As  an  illustration  of  the  applicability  of  Figure  41  to  a  low  but  non¬ 
zero  mass -density  ratio  case,  Figure  42  in  included  herein.  Figure  42 
shows  the  flutter  frequency  as  a  function  of  the  sweep  parameter  for  a  set 
of  structural  parameters  equivalent  to  those  of  Figure  39.  Two  of  the 
structural  parameters  used  to  construct  Figure  42  were  chosen  at  random 
as,  in  fact,  were  the  two  mass -density  ratios  plotted.  This  figure  may  be 
regarded  as  typical.  Figure  43  shows  the  corresponding  flutter  speed  co¬ 
efficients  . 

Returning  to  Figure  41,  several  interesting  facts  may  be  noted. 
First,  no  low  mass -density  ratio  flutter  is  predicted  if  the  elastic  axis  is 
forward  of  the  quarter  chord.  Second,  it  may  be  noted  that  the  structural 
configurations  most  likely  to  encounter  a  low  mass -density  ratio  flutter 
problem  (i.  e.  ,  for  minimal  values  of  K  and  the  largest  range  of  elastic 
axis  (a)  locations)  are  those  much  less  stiff  in  bending  than  in  torsion, 
i.  e.  ,  low  II*  values.  Conversely,  if  the  elastic  axis  is  just  aft  of  the 
quarter  chord,  virtually  all  combinations  of  stiffness  may  be  expected  to 
yield  valid  flutter  solutions  for  comparatively  low  values  of  the  sweep 
parameter.  However,  Figure  44  shows  that  the  flutter  speed  coefficient 
asymptotically  approaches  infinity  as  the  elastic  axis  is  assumed  to  move 
forward  toward  the  quarter  chord.  To  obtain  Figure  44,  the  sweep  para¬ 
meter  value  is  always  assumed  to  be  the  minimum  required  for  a  valid 
flutter  solution.  This  plot  illustrates  why  the  choice  of  a  forward  elastic 
axis  by  Herr  (Ref.  5)  was  previously  termed  unfortunate  from  the  stand¬ 
point  of  an  experimental  verification  of  his  theoretical  results.  Figure  45 
is  a  similar  plot  for  a  strut  assumed  to  retain  the  mass  properties  of  the 
Baird  model  #2  with  the  mass -density  ratio  assumed  to  be  .  1.  It  can  be 
seen  that  there  is  little  difference  between  the  general  zero  mass -density 
ratio  plot  and  the  specific  plot  of  Figure  45.  However,  as  the  mass- 
density  ratio  is  assumed  to  increase,  the  influence  of  the  location  of  the 
center  of  gravity  of  the  section  mass  may  be  expected  to  increase.  Ac¬ 
cordingly,  the  shape  of  the  curve  may  be  expected  to  be  altered  somewhat. 
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As  has  been  previously  noted,  the  zero  mass -density  ratio  analy¬ 
sis  shows  the  ••critical11  sweep  parameter  value  to  be  associated  with  a 
zero  flutter  frequency.  This  is  not  unexpected  as  a  similar  trend  has 
been  noted  and  has  been  discussed  for  the  several  specific  non -zero 
mass -density  ratio  cases  included  herein.  For  the  special  case  of  a 
zero  flutter  frequency,  the  HexactM  unsteady  hydrodynamic  theory  de¬ 
generates  to  the  simpler  quasi-steady  (k  0)  theory.  Hence,  it  is 
apparent  that  the  ••critical11  sweep  parameter  values  may  be  determined 
by  considering  the  quasi-steady  zero  mass -density  ratio  flutter  equation. 

If  the  quasi -steady  values  of  F  and  G  are  substituted  into 
Eqs.  (II-  11a)  and  (U-llb),  the  following  may  be  obtained. 

1  +  .3136(2(1  +  2a)X*n*  -  a(l  -  4a) 

«  3  =  -  (II-13a) 

1  +  8a3  +  II3  (8a3  ^  4a) 

-  K(^-  +  a3  ){l6a  -  4  +  K[4a  -  8a3]] 

O 


It  may  be  recalled  that  a  similar  simplification  of  the  flutter  equa¬ 
tions  of  the  "typical  sectionH(see  Part  I)  immediately  resulted  in  a  set  of 
uncoupled  equations.  However,  such  is  not  now  the  case  as  Eq.  (II-l3a) 
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contains  an  X*  term.  Thus  the  independent  determination  of  the  flutter 
frequency  coefficient  requires  that  Eq.  (II-13b)  be  substituted  into 
Eq.  (II-13a).  Performing  this  substitution  results  in  a  quadratic  in  £2* 
which  may  be  written  as  follows: 


f2*{4(l  +2a)njz  -K(l  +  8a3)  [l  -K(a+|)]Z  -.  3136(1  +  2a)  n*( l  +  8a3  )j 

+  f2f  [4Z  -  .6272(1  +  2a  -2aa)Z  3136(1  +  2a)(l  +  8a3  )[K[2a -K(^-+aa)3 

-(1+nJ)]  -[l  +  .3136a(4a  -1)][4(1  +  2a)IIa  -K(l  +  8aa)[l  -X(a+j)]]] 

-  .  3136(1  +  2a)(.6272aa  +  1)  -  [l  +  .  3136  a(4a  -  1)][4  -  .  6272(1  +2a  -2a3)]  =  0 

(11-14) 


where  Z,  the  denominator  of  Eq.  (II-13a),  is  written  as: 

Z  =  [1  +  8aa  +  na(8a3  -  4a)]  -  K{(^-  +  aa)(l6a  -  4)} 

V  Q 

-K*{(|-  +  a3)(4a  -  8a3  )J 

Eq.  (11-14)  may  be  easily  solved  by  the  quadratic  formula.  Since 
the  critical  sweep  parameter  value  results  in  a  zero  flutter  frequency, 
the  critical  value  may  be  determined  by  letting  the  quadratic  solution  of 
Eq.  (II- 14 )  approach  infinity.  For  finite  values  of  the  structural  para¬ 
meters,  an  infinite  solution  is  approached  as  the  coefficient  of  the  £2* 
term  approaches  zero.  This  coefficient  may  be  written  as  follows: 
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X4{(1  +  8aa)(2a  -  l)(i-+a3)(4a  -  8a3)} 

O 

+  K3  {(1  +  8a3  )(2a  -  1  )(i- +  a3 )( 1 6a  -  4)  +  (1  +  8a3  )(i-+a3  )(2a  -  1  )8a} 

O  O 

-  K3  {( 1  +  8a3)(2a  -  1)(1  +  8a3  +  4aII3[2a  -  l]) 

+  (1  +8a3)(i-+a3)(4a  -  1)8  -  (1  +  8a3  )(1  +2a)(l  -2a)4aIIa} 

+  X{(1  +8a3)(l  +2a)(4a  -1)4II3  +2(1  +8aa)(l  +  8aa  +  4aII3  [2a  -  1  ] )} 

-  8(1 +2a)(l +8a3  +  4an3[2a -1JT3  +  .6272(1  +  2a)(l  +  8aa)na  =  0 

(H-15) 


Hence,  the  above  may  be  solved  for  K  to  yield  Figure  41. 

It  may  be  seen  from  Figure  23  that,  in  general,  practical  sweep 
parameter  values  are  less  than  unity  while  Figure  41  shows  that  many 
critical  sweep  parameter  values  are  an  order  of  magnitude  less  than 
unity.  Accordingly,  the  quartic  and  cubic  terms  of  Eq.  (11-15)  may  be 
expected  to  be  negligible  for  many  cases  of  practical  interest.  This 
leaves  a  comparatively  simple  quadratic  expression.  The  remaining 
quadratic  expression  may  be  solved  for  K,  the  results  being  included 
herein  as  Figure  46.  As  might  be  expected.  Figures  46  and  41  agree 
very  well  for  those  combinations  possessing  a  low  K  critical  while  the 
agreement  deteriorates  for  those  combinations  of  II*  and  (a)  which  have 
a  comparatively  large  K  critical.  Extending  the  comparison  leads  to  the 
conclusion  that  the  simplified  flutter  equation  upon  which  this  investiga¬ 
tion  is  based  may  be  expected  to  predict  somewhat  erroneous  results  for 
configurations  possessing  an  aft  elastic  axis  and  an  equivalent  xero  mass- 
density  ratio  Min  fluid11  frequency  ratio  squared,  II*,  on  the  order  of  one 
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tenth.  However,  it  is  evident  that  the  fundamental  features  of  the 
phenomenon  of  low  mass -density  ratio  flutter  are  retained  even  in  the 
much  simplified  equation  used  to  obtain  Figure  46. 

II -H  A  Physical  Interpretation  of  the  Low  Mass -Density  Ratio 
Flutter  Mechanism 

If  the  results  of  the  parameter  studies  discussed  in  Part  II  of  this 
investigation  are  considered  as  a  whole,  an  explanation  of  the  physical 
origin  of  the  coupled  bending -torsion,  low  mass -density  ratio  flutter  of 
uniform,  cantilevered,  fully  wetted  lifting  surfaces  may  be  formulated. 
The  most  important  features  of  the  low  mass -density  ratio  solutions 
which  lead  to  the  explanation  to  follow  are: 

1.  A  positive,  non -zero  sweep  parameter  value  is  required 
for  a  low  mass -density  ratio  dynamic  instability  to  occur. 

This  implies  a  structure  of  finite  length,  the  elastic  axis 
of  which  is  swept  back. 

2.  The  structures  most  likely  to  encounter  a  low  mass -density 
ratio  flutter  problem  are  those  much  less  stiff  in  bending 
than  in  torsion. 

3.  The  structure  must  be  flexible  in  the  spanwise  direction 
and  properly  represented  as  such.  This  is  implied  by  the 
fact  that  an  off-diagonal  coupling  constant  less  than  unity 
is  required  for  the  proper  prediction  of  low  mass -density 
ratio  flutter,  and, 

4.  The  minimum  sweep  parameter  value  for  which  a  dynamic 
instability  can  occur  yields  a  zero  flutter  frequency.  Thus 
the  flutter  instability  degenerates  to  a  pure  divergence  as 
the  sweep  parameter  is  decreased. 

Item  1  suggests  that  the  streamwise  angle  of  attack  which  results 
from  the  bending  of  a  swept  strut  is  an  important  part  of  the  mechanism. 
Item  2  supplements  the  suggestion  of  Item  1  since  a  strut  much  less  stiff 
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in  bending  than  in  torsion  would  have  the  largest  ratio  of  streamwise 
angle  due  to  bending  to  streamwise  angle  due  to  torsion.  Item  3  suggests 
that  the  spanwise  pitch  to  bending  ratio  is  important,  while  Item  4  indi¬ 
cates  that  the  dynamic  instability  is  a  direct  outgrowth  of  a  tendency  of 
the  strut  to  diverge.  The  fact  that  no  low  mass -density  ratio  flutter  may 
be  predicted  for  struts  having  an  elastic  axis  forward  of  the  quarter  chord 
further  strengthens  the  latter  indication.  Hence,  it  is  suggested  that  the 
flutter  condition  occurs  in  the  following  manner.  First,  a  perturbation  in 
angle  of  attack  causes  a  moment  about  the  elastic  axis  which  tends  to 
further  increase  the  angle  of  attack  of  the  strut.  This  in  turn  results  in 
an  incremental  increase  in  lift  which  tends  to  bend  the  strut,  while  the 
translation  resulting  from  the  lift  force  tends  to  produce  an  opposing 
streamwise  angle  of  attack.  Since  the  relieving  angle  of  attack  is  a 
function  of  the  local  bending  slope  and  the  translational  velocity,  it  is 
not  uniform  over  the  span.  Hence,  the  instantaneous  spanwise  leading 
is  also  not  uniform  -  and  may  even  be  non -unidirectional.  Therefore, 
it  is  suggested  that  the  flutter  condition  develops  -  at  least  for  the  low 
flutter  frequency  cases  associated  with  sweep  parameter  values  near 
the  critical  value  -  when  the  positive  hydrodynamic  spring  on  the  outer 
end  of  a  strut  overcomes  the  negative  hydrodynamic  spring  of  the  fixed 
end  and  causes  an  oscillatory  motion  to  develop.  Obviously,  the  direc¬ 
tion  of  the  hydrodynamic  spring  is  a  function  of  the  streamwise  angle  of 
attack.  In  all  probability,  the  oscillation  is  sustained  by  the  inertia  ef¬ 
fects  of  the  structural  and  virtual  masses  and  the  negative  hydrodynamic 
spring  which  exists  over  a  portion  of  the  strut.  Naturally  the  balance  of 
forces  must  be  such  that  a  complete  cycle  results  in  positive  work  being 
performed  on  the  strut. 

It  may  be  noted  that  the  rather  delicate  balance  of  loading  re¬ 
quired  to  produce  the  condition  described  above  is  a  direct  function  of 
the  local  spanwise  pitch  to  bending  ratio.  This  would  account  for  the 
extreme  sensitivity  of  the  analyses  to  the  vibratory  mode  shapes  assumed 
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in  the  representation  of  the  structure.  In  addition,  the  proposed  neces¬ 
sity  of  having  both  a  positive  and  a  negative  hydrodynamic  spring  acting 
on  the  strut  would  seem  to  explain  why  the  production  modal  analyses 
such  as  those  discussed  in  Ref  10,  have  always  shown  the  second  bend¬ 
ing  mode  branch  to  be  the  critical  branch  leading  to  the  flutter  condition. 

If  the  inertia  load  of  the  translating  virtual  mass  which  is  cen¬ 
tered  at  the  midchord  is  considered  in  even  the  most  rudimentary  fashion, 
it  may  be  seen  that  the  proposed  mechanism  accounts  for  ther  shape  of 
the  curves  shown  in  Figure  41.  It  may  be  seen  from  Figure  41  that  the 
highest  critical  sweep  parameter  values  are  associated  with  parametric 
combinations  having  an  elastic  axis  location  near  the  midchord  position. 
Clearly,  the  effect  of  the  inertia  force  of  the  virtual  mass  is  such  that 
it  tends  to  increase  the  diverging  angle  of  attack  for  forward  elastic 
axis  locations  while  it  has  a  relieving  effect  for  the  aft  axis  configura¬ 
tion.  Hence,  the  aft  axis  configurations  require  that  less  of  the  total 
relieving  angle  be  provided  by  the  pitch-bending  effect. 

It  must  be  once  again  mentioned  that  the  proposed  flutter  mech¬ 
anism  has  been  formulated  along  the  lines  of  the  streamwise  approach 
followed  throughout  this  investigation.  However,  there  is  little  doubt 
that  an  entirely  equivalent  explanation  can  be  formulated  from  the  ”axis- 
normal"  approach.  Further,  it  must  be  stated  that  the  above  explanation 
does  not  consider  the  unsteady  fluid  forces  which  are  the  natural  result 
of  an  oscillatory  motion.  Hence,  it  is  possible  that  the  mechanism  may 
be  somewhat  more  complex  as  the  sweep  parameter  is  increased  sub¬ 
stantially  beyond  the  critical  value  as  the  increase  is  accompanied  by  a 
rise  in  the  flutter  frequency. 

Figure  47  shows  the  calculated  divergence  or  flutter  speeds  as 
a  function  of  the  sweep  angle  for  Squires’  strut  model.  It  may  be  noted 
that  the  divergence  speed  of  the  unswept  strut  is  considerably  less  than 
the  minimum  flutter  speed,  which  is  as sociated  with  K  critical.  Landahl, 
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et  al,  showed  a  similar  plot  in  Ref.  26,  for  the  Herr  #2  model.  When 
these  results  are  combined  with  the  proposed  theory,  one  is  drawn  to 
the  conclusion  that  an  unswept,  uniform,  cantilevered  strut  which  will 
not  diverge  at  some  given  free  stream  velocity  will  not,  for  any  reason¬ 
able  sweep  angle,  experience  a  low  mass -density  ratio  flutter  problem 
for  at  least  as  great  a  free  stream  velocity.  This  statement  is  made  be¬ 
cause  the  flutter  condition  appears  to  develop  directly  from  an  initial 
attempt  of  the  structure  to  diverge. 

II -J  Conclusions  Resulting  from  Part  II 

Part  II  of  this  investigation  has  dealt  specifically  with  fully  wetted, 
uniform,  cantilevered  lifting  surfaces.  Further,  it  has  been  assumed 
that  the  surfaces  could  be  considered  rigid  in  the  streamwise  direction 
with  the  span  to  chord  ratio  large  enough  to  justify  the  use  of  two- 
dimensional  strip  theory.  Accordingly,  the  specific  conclusions  to  be 
stated  below  may  be  considered  to  be  entirely  valid  only  when  dealing 
with  such  structures. 

However,  this  last  statement  is  not  meant  to  imply  that  a  low 
mass-density  ratio  flutter  mechanism  such  as  that  discussed  in  the  pre¬ 
vious  section  cannot  exist  for  any  structures  other  than  those  supported 
in  a  cantilevered  fashion.  Clearly,  other  types  of  structures  could  ex¬ 
perience  a  non-unidirectional  spanwise  force  distribution;  the  only  re¬ 
quirement  being  a  spanwise  variation  in  the  pitch  to  plunge  ratio. 

Rather,  the  reason  for  qualifying  the  conclusions  stated  in  this  section 
is  simply  that  the  sensitivity  of  the  flutter  mechanism  to  small  para¬ 
metric  changes  could  make  certain  extrapolations  of  the  results  pre¬ 
sented  in  Part  II  of  questionable  value.  Some  extrapolations  and  gene¬ 
ralisations  of  the  results  which  are  believed  to  be  valid  will  be  treated 
in  the  next  section  which  deals  with  the  conclusions  resulting  from  the 
consideration  of  this  entire  investigation  in  combination  with  the  results 
obtained  by  other  investigators. 
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As  a  direct  result  of  Part  II  of  this  investigation,  the  following 
conclusions  may  be  drawn: 

1  A  comparatively  simple  modal  analysis  using  the  funda¬ 
mental  bending  and  torsion  modes  as  the  assumed  flutter 
modes  may  be  expected  to  predict  the  existence  of  a  poten¬ 
tial  low  mass -density  ratio  flutter  problem  provided  the 
spanwise  variation  in  coupling  between  the  pitch  and  bend¬ 
ing  modes  is  considered  and  further  provided  that  any 
sweep  of  the  elastic  axis  is  considered.  Conversely,  the 
prediction  will  not  be  reliable  if  the  system  is  reduced  to 
the  so-called  "typical  section,  "  as  one  or  more  features 
important  to  the  fundamental  flutter  mechanism  will  be 
lost.  In  particular,  the  local  spanwise  pitch  to  bending 
ratio  is  an  important  parameter.  It  further  follows  that 
semi-rigid  strut  models  are  not  susceptible  to  low  mass- 
density  ratio  dynamic  instabilities  of  the  type  discussed 
herein,  and  hence,  tests  performed  with  semi-rigid 
models  are  of  questionable  value. 

2.  Although  the  simplified  analysis  mentioned  above  is  suf¬ 
ficient  to  predict  the  existence  of  a  potential  problem,  the 
sensitivity  of  the  flutter  speed  prediction  to  small  changes 
in  the  assumed  vibratory  character  of  the  structure  im¬ 
plies  that  an  accurate  prediction  of  flutter  speed  may  re¬ 
quire  a  relatively  complex  analysis  involving  a  considerable 
number  of  modes  and  panel  points.  Again,  this  is  the  result 
of  the  extreme  dependence  of  the  mechanism  upon  the  local 
spanwise  pitch  to  bending  ratio. 

3.  No  low  mass -density  ratio  flutter  may  be  predicted  -  and 
none  is  believed  to  be  physically  possible  -  unless  the  struc¬ 
ture  has  a  positive,  non -zero  sweep  parameter  value.  For 
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physically  plausible  structures,  this  implies  that  the 
elastic  axis  must  be  swept  back  relative  to  the  normal 
to  the  free  stream.  The  minimum  sweep  parameter 
value  for  which  a  dynamic  instability  can  occur  yields 
a  zero  flutter  frequency. 

4.  The  streamwise  relieving  angle  of  attack  of  a  swept 
cantilever  subjected  to  bending  is  an  important  part 

of  the  flutter  mechanism.  Accordingly,  the  structures 
most  likely  to  encounter  a  low  mass -density  ratio  flut¬ 
ter  problem  are  those  much  less  stiff  in  bending  than 
in  torsion. 

5.  No  low  mass -density  ratio  flutter  may  be  predicted  - 
and  none  is  believed  to  be  physically  possible  -  for 
structures  having  an  elastic  axis  forward  of  the  quarter 
chord. 

6.  A  structure  which  may  experience  a  low  mass-density 
ratio  flutter  problem  will  encounter  the  instability  at 
approximately  the  same  dynamic  pressure  throughout 
the  low  mass -density  ratio  range. 

It  is  further  concluded  that  the  analysis  of  the  limiting  zero 
mass -density  ratio  case  which  has  been  shown  herein  allows  a  rapid 
determination  of  whether  a  structure  can  experience  a  low  mass -density 
ratio  dynamic  instability  at  any  finite  free  stream  velocity  and  further, 
that  parametric  trends  obtained  from  the  analysis  of  the  limiting  case 
are  qualitatively  valid  throughout  the  low  mass -density  ratio  range. 

In  short,  it  may  be  stated  that  the  low  mass -density  ratio  flutter 
cases  which  motivated  this  study  (see  e.g.  ,  Ref.  10)  can  be  said  to  be 
both  qualitatively  and  quantitatively  understood  as  a  direct  result  of  this 
investigation. 
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GENERAL  CONCLUSIONS  AND  RECOMMENDATIONS 

If  the  results  of  this  investigation  and  the  previously  documented 
experiences  of  other  investigators  are  considered  as  a  whole,  one  is 
drawn  to  the  conclusion  that,  in  principle,  the  existing  unsteady  hydro- 
dynamic  and  flutter  theories  may  be  used  to  predict  the  low  mass -density 
ratio  flutter  behavior  of  fully  wetted  lifting  surfaces.  However,  the 
physical  origin  of  the  low  mass -density  ratio  flutter  phenomenon  is  such 
that  the  hydro -elastic  behavior  of  a  lifting  surface  is  extraordinarily  de¬ 
pendent  upon  its  vibratory  character.  Hence  an  accurate  prediction  re¬ 
quires  a  much  more  accurate  structural  representation  that  what  might 
be  expected  on  the  basis  of  well  documented  aircraft  experience.  Con¬ 
versely,  it  may  be  stated  that  the  early  low  mass -density  ratio  investiga¬ 
tions  which  were  based  upon  the  mathematical  model  of  a  ’’typical  section” 
were  not  fundamentally  wrong.  However,  the  conclusions  which  resulted 
from  such  studies  were  extended,  on  the  basis  of  aircraft  experience,  to 
such  a  degree  that  completely  erroneous  notions  followed.  It  follows  that 
the  most  important  lesson  to  be  learned  from  all  of  the  low  mass -density 
ratio  work  performed  to  date  is  that  the  region  is  typified  by  the  extreme¬ 
ly  sensitive  nature  of  the  hydro -ela  Stic  behavior  of  a  system  to  small 
changes  in  the  system  parameters.  Thus,  it  is  a  region  for  which  it  is 
dangerous  to  generalize.  It  is  the  very  sensitivity  of  the  dynamic  be¬ 
havior  of  lifting  surfaces  operating  in  low  mass -density  ratio  region  that 
suggests  some  areas  which  might  be  fruitfully  investigated  in  the  future. 

Both  Part  I  and  Part  II  of  this  investigation  showed  qualitatively 
the  sensitivity  of  /i  or  the  /i  peak  to  small  changes  in  the  system  para- 

cL 

meters.  It  is  suggested  that  a  quantitative  investigation  of  this  phenome¬ 
non  might  be  valuable.  If  reliable  theoretical  techniques  could  be  de¬ 
veloped,  it  would  be  possible,  in  certain  cases,  to  make  a  slight  struc¬ 
tural  change  and  greatly  improve  the  hydro-elastic  behavior  of  a  surface. 


-  125  - 


It  is  likely  that  an  investigation  such  as  the  one  suggested  would  require 
that  considerable  effort  be  expended  upon  the  probable  improvements  in 
theory  made  possible  by  the  use  of  three  dimensional  unsteady  deriva¬ 
tives. 

As  an  extension  of  the  work  presented  in  Part  II,  it  is  suggested 
that  a  cantilevered  structure  with  an  aft  mounted  control  surface  be  con¬ 
sidered,  as  many  aft  mounted  hydrofoil  struts  are  of  this  design.  In 
addition,  the  transient  response  of  typical  hydrofoil  designs  could  be  in¬ 
vestigated  to  determine  whether  fatigue  problems  might  exist  even  when 
an  actual  instability  does  not  occur. 

All  of  the  work  performed  to  date  has  considered  the  structure  to 
be  rigid  in  either  the  streamwise  or  axis -normal  directions.  In  view  of 
the  demonstrated  sensitivity  of  the  predicted  hydro -ela Stic  behavior  of  a 
structure  to  its  assumed  vibratory  character,  a  study  considering  the 
possible  effects  of  a  completely  elastic  structure  would  seem  to  be  in 
order. 

Lastly,  it  must  be  once  again  stated  that  almost  all  of  the  work 
performed  to  date  has  been  concerned  with  fully  wetted  surfaces.  As 
higher  speeds  are  attained,  it  will  be  necessary  to  investigate  the  effects 
of  cavitation  upon  the  hydro -elastic  behavior  of  a  system. 
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APPENDIX  A 


An  Abbreviated  Development  of  the  Complex  MFlutter  Pete rminant"  of 
a  Two-Degree  -of  -Freedom  Section  Oscillating  in  Incompressible  Flow 


The  equations  of  motion  of  the  representative  section  of  Figure  2, 
may  be  written  in  the  following  form: 

mh  +  S  a  +  m.tu^h  =  -  L» 
a  h 

(A-l) 

3.  h  +  I  a  +  I  uu  3a  =  M 

a  a  a  a 


The  determination  of  the  aerodynamic  forces  which  serve  as  the 
forcing  functions  of  the  above  set  of  equations  is  a  lengthy  one.  However, 
it  has  been  well  documented  in  a  number  of  references  including  Refs.  1, 
13,  14,  15,  1 6 ,  and  hence  need  not  be  considered  here  in  any  detail.  Quali¬ 
tatively  it  may  be  mentioned  that  for  the  case  of  unsteady  motion  the  lift 
on  a  foil  section  is  no  longer  a  simple  function  of  the  circulation.  For 
one  thing,  hydrodynamic  inertia  forces  and  moments  can  result  due  to 
the  acceleration  and  deceleration  of  the  fluid  moving  with  the  foil  -  even 
at  zero  forward  velocity.  Secondly,  there  is  a  contribution  known  as  the 
quasi -steady  lift.  This  is  the  lift  which  would  be  produced  by  the  motion 
of  the  foil  if  the  wake  were  assumed  to  have  no  effect.  The  third  con¬ 
tribution  may  be  regarded  as  a  correction  to  the  quasi -steady  lift  result¬ 
ing  from  considering  the  normal  velocities  in  the  neighborhood  of  the  foil 
which  are  produced  by  the  trailing  vortices.  Hence,  the  lift  and  moment 
are  a  function  of  the  motion  and  the  configuration  of  the  system. 

Theodor  sen  and  others  considered  the  motion  of  the  foil  at  the  in¬ 
ception  of  flutter  to  be  simple  harmonic  and  after  some  rather  involved 
mathematics  showed  that  the  aerodynamic  forces  and  moments  can  be 
expressed  as: 
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L  =  irpb*[L  +  U  o  -baa]  +  2wpU  bC(k)[h  +  U  a  +b(y-a)a] 

I  I  I  w 

M  =  lrpb8[bah  -  U,b(|-  -  a)a  -  b8  (|-  +  a8  )a  ] 

1  L  o 

+  277pUfb8(a+j)C(k)[h  +  U  a-  +  b(j  -  a)a] 

(A -2) 


where  C(k)  is  the  well  known  Theodorsen  function  expressible  as  the 
couplex  number  F  +  j G.  The  Theodorsen  function  is  a  function  of  only 
the  reduced  frequency  k.  F  and  G  may  both  be  expressed  as  a  ratio  of 
Bessel  functions  of  the  first  and  second  kinds,  of  argument  k.  C(k)  can 
also  be  expressed  in  terms  of  Hankel  functions.  The  dependence  of  C(k) 
upon  k  can  most  easily  be  shown  by  a  complex  polar  plot. 


.  5  F<k> 


Figure  A-l 
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Eq.  (A -2)  may  now  be  combined  with  Eq.  (A-l)  and  the  harmonic 
motion  assumption  again  introduced  to  yield  the  following: 

-  u^mh  -  tu^mh  =  -  TT.p  ba  [  -o>®  h  +  jU^UM?  +  tt)^ba] 

+  2JrpUfbC(k)[ja>fh  +  UfQ’  +  jwfb(j  -  a)a  ] 


-a)fsStth-«»f*I +%SIaQ'  =  tfpbs[-<i>*bah-ju>fUfb(|--a)tt  +  <oaba(|-+aa  )a 

+  ff2pUftf(a  +~)C(k)[ji»fh  +  UfO’  +  j«fb(l--a)or] 

(A -3) 

The  lift  and  moment  expressions  may  be  put  into  the  notation  of 
Ref.  16  by  the  introduction  of  the  unsteady  derivatives,  this  yields  the 
following  set  of  expressions* 

- -■  mh  -  »,*S  a  +  oo  amh  =  ffpb3®/  { L.  ~  +  [L  -  LJv  +  a)]^  } 
f  f  or  h  fhb  a  h  2 

-  u>  aS  h-0)aI  O'  +  <u  a I  a  =TT  pb*tu/  {[M,-  LJ^  +  a)] 
f  a  {  a  a  a  f  h  h  2  b 
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Each  of  the  above  equations  may  now  be  suitably  non -dimens  ion  - 

alized  by  dividing  the  first  by  -  7T  pb3(i)f8  and  the  second  by  -7rpb4cuf3  and 

introducing  the  non-dimensional  expressions  x  and  r  3. 

a  a 


m  h  m  m  h  h 

ff pb*  b  irpb*  *C,C‘  'up b*  »•  b 


-Lhr-[L»-Lh(r+a|]“ 


m  h,  m  9  O'  m  8  r,,  .  ,1  .  ,n  h 

-  x  —  + -  r  a - r  or  =  -  [M,  -  L  (—  +  a)J  — 

_  ,  s  o  b  3  a  8  a  O'  h  h  2  b 

If pb  IT  pb  (ju*  ttpb 


-[M  -(L  +M  )(j+a)+L  (i+a)8]^ 
o'  h  2  h  2 


(A -5) 


Transposing,  collecting  coefficients  of  the  same  motion  coordi- 
m 

nates,  and  letting  [X  =  .  “  the  mass -density  ratio  -  will  yield  the  fol¬ 

lowing  set  of  homogeneous  algebraic  equations: 


«.»  o>  8 


CuCl  -  —  •  —  +  L,  )  £  +  fax  +  L  -  L,  4  +  a)}o-  = 
3  8  h  b  a  ot  h  2 


0)  U)„ 

or  f 


m 

-Jr'i  +  M„ 


•(L«+Mhl(r+a)  +  Lh(r+a)*1‘' =  0 


(A -6) 


For  a  meaningful  solution,  the  determinant  of  the  above  set  must 
vanish.  This  determinant  is  the  familiar  "flutter  determinant"  of 
Eq.  (1-4). 
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APPENDIX  B 


An  Abbreviated  Development  of  a  Binary  Flutter  Determinant  for  a 
Uniform,  Swept,  Cantilever  Lifting  Surface 

The  calculation  of  the  steady  lift  distribution  of  a  uniform,  swept, 
cantilever  lifting  surface  of  large  aspect  ratio  presents  little  difficulty 
apart  from  the  root  and  tip  since,  if  the  flow  is  assumed  to  be  friction  - 
less,  the  spanwise  component  of  the  free  stream  has  little  influence. 
However,  such  is  not  the  case  when  the  unsteady  hydrodynamic  forces 
acting  on  such  surfaces  are  of  interest  as  the  downwash  will,  in  general, 
vary  along  the  span  of  an  oscillating  surface  even  if  the  root  and  tip 
effects  are  neglected. 

Ideally,  an  unsteady  hydrodynamic  theory  applicable  to  a  swept 
surface  should  include  all  of  the  effects  mentioned  above.  However,  an 
approach  which  is  both  simplifying  and,  in  general,  conservative  is  to 
neglect  the  root  and  tip  effects  when  dealing  with  surfaces  which  have  a 
moderate  to  large  aspect  ratio  Thus,  a  two-dimensional  "strip"  theory 
is  assumed  to  apply.  This  strip  theory  is  similar  to  what  is  commonly 
used  in  the  flutter  calculations  of  unswept  surfaces. 

When  the  strip  method  is  to  be  applied  to  a  swept  surface,  the 
question  arises  as  to  whether  the  strips  should  be  taken  normal  to  the 
elastic  axis  or  in  the  streamwise  direction  For  reasons  which  have 
been  discussed  in  the  text,  the  latter  approach  is  used  in  the  abbreviated 
development  to  follow.  Ref.  20  may  be  consulted  if  more  details  are 
desired  , 

With  regard  to  the  surface  being  considered,  the  effects  of  sweep¬ 
ing  a  surface  on  the  relation  between  the  forces  and  the  flutter  speed  can 
be  considered  in  a  convenient  manner  if  the  axial  length  (unswept  span) 
and  the  normal  chord  are  taken  to  be  invarient  during  the  sweeping  pro¬ 
cess.  The  retention  of  these  dimensions  can  be  accomplished  if  the 
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surface  is  swept  by  rotation  in  the  manner  described  in  the  text  and 
shown  in  Figure  22. 

For  this  development,  the  point  of  origin  of  the  axis  system  is 

taken  to  be  the  root  of  the  elastic  axis,  i.  e.  ,  point  0  in  Figure  22.  x 

0 

and  y  are  non-dimensional  coordinates.  The  length  xb  sec  A  is  meas¬ 
ured  in  the  streamwise  direction  from  the  elastic  axis,  positive  back¬ 
ward.  The  length  y£  is  measured  along  the  elastic  axis,  b  and  t  are 
both  taken  as  defined  in  the  general  list  of  symbols,  i.  e.,  invariant 
with  sweep.  The  lift  per  unit  span  L»(  y  )  is  taken  to  be  opposite  to  the 
positive  normal  displacement  z  and  a  nnose-upM  twist  is  regarded  as 
positive  while  the  moment  M(  y  )  is  positive  ’’nose -down.  "  It  is  to  be 
understood  that  the  forces  involved  imply  a  harmonic  motion  of  the 
structure.  Hence,  any  reference  to  the  harmonic  factor  is  omitted  in 
the  subsequent  treatment. 

Now  if  the  lift  per  unit  span  is  L.(  y  )  and  the  torsional  moment 
per  unit  span  is  M(  y  ),  the  force  and  moment  acting  on  the  elastic  axis 
per  unit  length  of  the  elastic  axis  are: 

L»(  y  )  =  L»(  y  )  cos  A 

and  (B-l) 

M(  y  )  =  M(  y  )  cos8  A 

The  lift  and  moment  are  expressed  in  terms  of  non-dimensional 
load  functions  as: 

L(  y  )  =  bjWfx,  y  )dx  =  pcuf*b3Zo(  y  ) 

and  (B-2) 

M(  y  )  =  baJx1F(x,  y  )dx  =  p<t)*b*Zi  (y  ) 

The  subscripts  o  and  i  refer  to  translation  and  pitch  or  bending 
and  torsion  respectively  as  indicated  in  the  supplementary  list  of  symbols. 

See  supplementary  list  of  symbols  at  end  of  Appendix  B. 
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The  deformations  of  the  elastic  axis  are  represented  by  the  non- 

dimensional  "modes'1  z  (  y  ). 

s 

Vertical  Translation  =  z  (y’)b 

o 

Nose  up  Angle  of  Twist  =zx(y)  (B-3) 

The  local  stiffnesses  are  represented  by  the  function  a  and  o 

01 

of  dimensions  time  as  follows: 

Bending  Stiffness  =  El  =  pb3£4  of 
Tor sional  Stiffness  =  GJ  =  pb4£3o"  (B-4) 

By  combining  expressions  B  -  2,3,4,  the  elastic  equations  re¬ 
lating  the  deformations  of  the  elastic  axis  to  the  external  forces  may  be 
written.  These  equations  are  written  by  considering  the  fact  that  the 
bending  moment  M(  y  )  tan  A  acts  on  the  axis  in  addition  to  the  bending 
force  L»(  y  ).  Further,  the  equations  below  are  written  by  considering 
the  phase  angle  of  the  structural  damping  to  be  zero  and  without  regard 
to  the  boundary  conditions  at  the  point  y  =  0. 

J*  [(n-y)Z  (ti)  +  |z  frpldt,  =  -  -2- - 2 - 

-  o  2  1  a  d_3 

y  f  y 

(B-5) 

o,  d[zl(y)] 

U}fS  dy 

In  the  above,  r)  is  used  as  a  dummy  variable  along  the  y  axis. 
Equations  (B-5)  refer  to  quantities  measured  in  the  axis -normal 
system.  However,  the  loads  Z  (y  )  are  referred  to  the  external  forces 
on  the  streamwise  strip  through  Eq.  (B-2). 


J*  Z1  (  T|  Ml)  ■  - 

y 
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The  coefficient  X  which  appears  above  has  been  termed  the  "sweep 
parameter"  in  the  text.  It  is  written  as: 

K  =  2  tan  A  (B-6) 

A  relationship  between  the  displacement  z(x,  y  )b  along  the  stream- 
wise  strip  and  the  modes  zg(y  )  of  the  elastic  axis  may  be  written  as  fol¬ 
lows: 

d[zo(7)b]  xb 

z(x,  y  )b  =  z  (  y  )b  +  xb  tan  A - = - + - 7  cos  A  z,  (  y  ) 

0  d(  y  l)  cos  A  1 

Dividing  by  b  and  noting  the  second  term  of  the  right  hand  side 
contains  the  "sweep  parameter,  "  the  above  may  be  written: 

v  d[«  (y)] 

z(x,y)  =  z  (  y  )  +  x[-  -  +z  (  y  )]  (B-7) 

0  *  dy  1 


Having  obtained  the  elastic  and  the  kinematic  relationships,  the 
determination  of  the  external  forces  on  the  streamwise  strip  which  are 
caused  by  a  given  displacement  of  the  strip  remains. 

The  external  forces  consist  of  inertia  forces  and  hydrodynamic 
forces.  All  of  these  forces  may  be  combined  by  writing  the  local  ex¬ 
ternal  force  as: 


F(x,  y)  =  -  w*m(x,  y)z(x,  y)b  (B-8) 

'V 

Here  m(x,y)  is  not  the  local  structural  mass  but  rather  that 
"complex  mass"  which,  if  fixed  to  the  structure  at  the  point  (x,y),  would 
produce  the  inertia  force  F(x,  y).  Hence,  m(x,  y)  is  the  sum  of  the  struc¬ 
tural  mass,  virtual  hydrodynamic  mass,  and  further  components  corres¬ 
ponding  to  hydrodynamic  stiffness  and  hydrodynamic  damping,  (see 
Appendix  A). 
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The  displacement  z(x,y")  of  the  streamwise  strip  is  now  split  up 
into  elementary  deformations,  x*1,  assuming  a  rigid  streamwise  section. 


z(x,  y)  =  L  z  (  y  )xn 
n  =0  n 


(B-9) 


Note  that  in  the  above  functions  z  (  y  )  are  non-dimensional  ampli¬ 
tudes  of  the  streamwise  strips  and  not  the  "modes"  of  the  elastic  axis. 

In  integrating  F(x,y)  along  the  streamwise  strip,  use  maybe  made 
of  Eq.  (B-2).  Thus  the  forces  on  the  elastic  axis  are: 

~T.  E. 

P(U  az  (  y  )  =  J’xrF(x,  y)dx  =  -03  8b2  z  (y)  |  m(x,  y  )xF+ndx 

f  r  f  n  n 


L.  E. 


Now  from  the  above  and  by  replacing  the  "complex  mass"  by  the 

coefficients  A  ,  the  following  expression  may  be  obtained 
rn 

XT.  E 

-pbZ  (  y  )  =  pb  2  A  z  (y)=£z  (y)  \  m(x,  y)xr+Rdx  (B-10) 

r  n  rn  n  n  n  \ 


E.E. 

The  "complex  mass"  not  only  depends  upon  the  reduced  frequency 

k  but  also  on  the  motion  x*1  being  considered.  Hence,  in  general,  A  ^A 

rn  •  n: 

However,  under  certain  conditions,  (see  Ref.  20)  the  "complex  mass"  is 
invariant  with  sweep  just  as  is  the  structural  mass.  The  physical  explana¬ 
tion  for  the  invariability  of  the  hydrodynamic  reactions  is  that  if  the  sweep 
ing  is  accomplished  by  rotation,  the  cos  A  term  which  is  the  result  of  the 
most  fundamental  effect  of  sweep  is  just  cancelled  by  the  sec  Afactor  in 

the  semichord,  i.e.  ,  bA  =  b  sec  A. 

A 
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By  combining  the  results  which  have  been  shown;  the  equation  re¬ 
lating  the  modes  z  (  y  )  to  the  load  functions  Z  (  y  )  may  be  obtained.  For 
s  r 

streamwise  strips,  these  may  be  written: 


_  _  _  _  v  dCz  (n)] 

-Zr«”  >  *  »  +  Ari  1 1  — ft—  +  *1  < 


(B-ll) 


The  flutter  equations  result  if  Eqs.  (B-ll)  are  combined  with 
Eqs.  (B -5 ). 

The  complete  set  of  flutter  equations,  for  a  given  reduced  fre¬ 
quency,  is  a  complex  eigenvalve  problem  with  the  unknown  eigen-solu¬ 
tions  being  the  complex  modes  z  (y  ).  In  most  flutter  calculations  these 

s 

modes  are  replaced  by  a  set  of  dynamic  coordinates. 

For  the  binary  case  being  considered  here,  the  modes  are  re¬ 
placed  by  one  coordinate  each  in  bending  and  torsion. 


zG(y)  =  qQh(y ) 
\  ( y")  =  sx a  ( y ) 


(B-12) 


Now  the  amplitudes  q  and  q^  become  the  unknowns  in  the 
problem. 

For  simplicity,  and  because  the  case  discussed  in  the  text  was 
that  of  a  rigidly  supported  cantilever,  the  boundary  conditions  which 
are  assumed  are: 


h(o)  =  — ^  (o)  =  a(o)  =  0 
dy 


The  Lagrangian  equations  can  now  be  obtained  if  the  elastic 

equations,  (Eqs.  B-5),  are  multiplied  by  — —  and  f  respectively, 

dy*  dy 

and  then  integrated  along  the  span  0  *  y  *  1. 
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Integration  by  parts,  using  the  boundary  conditions  previously 
stated,  results  in  the  following. 

i  i  8 

J*  fczj  y)+77z.(y  )^dy  +—  J'  (^-^)  dy  =  o 

o  £  dv  1  2  O  ,  —9  ' 

o  y  mod  v 


0)  o  d  y 


( B  —  1 3 ) 


l  l 

q  g 

J’  aZ((  y  )dy  +  —  J*  a  (  “ *  )  dy  =  0 
o  u^8  o  dy 


The  integrals  on  the  right  side  of  the  above  set  represent  the 
stiffness  and  can  be  expressed  by  the  uncoupled  "in  vacuum"  natural 
frequencies  a)^  and  u>^  .  This  results  in  the  following: 


J[hZo(y)+!  ^Z^ytfdy  +q<>-^-  V  J  Hh8dy  = 


8  1 


tt>£  O 


(B  -14) 


9  1 


I  ttZJ  y  )dy  +  qx  V  J*  u  r^a^dy  = 


(1)^  o 


Eqs.  (B-14)  are  the  binary  flutter  equations.  For  a  uniform 
listing  surface,  Eqs.  (B-14)  and  Eqs.  (B-ll)  can  be  combined  to  yield 
the  following  flutter  determinant. 


V  v* 

A  +(C  A  +c  A  )^+c  A  -A - — u 

oo'  1  01  9  lor2  3  1  1  4  9 


c  A  +  c  A  — 
4  ol  s  1  1  2 


c .  A  +  c.  A  — 
4  1 o  6112 


A  oc  8 

A - Llr 


(B-15) 
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In  Eq.  (B-15)  the  ^Qetc.  are  the  "complex  mass"  coefficients, 
or  more  precisely,  the  individual  structural  mass,  virtual  mass,  and 
hydrodynamic  reaction  terms  of  the  "typical  section"  or  unswept  analy¬ 
sis,  i.e.,  +  L  ,  etc.  (See  Appendix  A.  )  The  coefficients  c  are 

h  1 

integrals  taken  over  the  mode  shapes  and  can  be  calculated  numeri¬ 
cally  for  a  given  design.  Jordan  (Ref.  20)  calculated  the  values  for  a 
uniform  surface  by  using  as  the  normal  modes  h  and  a  the  fundamental 
modes  of  a  uniform  beam.  This  yielded  the  following  values. 


c  =4.648 

3 

c4  =  .959 

=  “  2.146 

It  may  be  noted  that  if  the  above  coefficients  are  used  Eq.  (B-15) 
becomes  the  fundamental  determinant  used  in  Part  II,  (see  Eq.  II-4). 
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Additional  Symbols  Used  in  Appendix  B 

A  "Complex  mass"  coefficient,  see  Eq.  (B-15) 

rn 

F(x,  y)  External  force  per  unit  area 

L  Lift  per  unit  span 

M  Moment  per  unit  span 

m(x,  y)  Structural  mass  per  unit  area 

m(x,  y)  "Complex  mass"  per  unit  area,  see  Eq.  (B-8) 

x,  y,  z  Non-dimensional  coordinates 

z^(  y  )  Bending  mode,  non-dimensional 

zi  (  y  )  Torsional  mode,  non -dimensional 

Zo(y)  Non-dimensional  lift  per  unit  length  of  elastic  axis 

(  y  )  Non-dimensional  moment  per  unit  length  of  elastic  axis 

T|  Dummy  variable  equal  to  y  in  integrals 

Bending  stiffness  function,  invariant  for  uniform  structure 

a  Torsional  stiffness  function,  invariant  for  uniform  structure 

l 

Indices 


i  General  index  (i  =0,1,2....) 

r  Denotes  force  (o  lift;  1  moment) 

s  Denotes  elastic  axis  deformation  (o  bending;  1  torsion) 

n  Denotes  streamwise  strip  deformation  (o  translation;  1  pitch) 


A  superimposed  bar  denotes  quantities  referred  to  "axis -normal"  direc¬ 
tions  except  for  "axis -normal"  or  unswept  length  t  and  semi -chord  b. 
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Figures  11a -llg  show  how  the  "critical" 
mass -density  ratio  varies  as  a  function 
of  the  location  of  the  elastic  axis  for 
various  center  of  gravity  positions  and 
a  non-dimensional  radius  of  gyration 
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fi^uke  11  au 


/6^-jl  VS-  cL  FOR 


life 


FIGURE.  11c 


1 5  i 


VS.  2L  For  VARIOUS  If2- 


Xm  -  .  Z& 
r|  =  •  25 


FUjUSE  lldL 


-154- 


vs.  X  Foe 


.1 1- 

-4 


-3 


-  2. 


va?ioo$  'TT2' 


5o 


Figures  13a-13i  show  what  combinations  of 
the  frequency  ratio  and  the  elastic  axis  loca¬ 
tion  will  yield  a  finite,  positive  value  for  the 
"critical"  mass-density  ratio.  These  plots 
are  for  various  locations  of  the  center  of 
gravity  and  for  a  non-dimensional  radius  of 
gyration  squared  of  .25. 
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Figures  l6a-l6o  show  the  quasi -steady  and 
qua  si -unsteady  "critical"  mass-density  ratios 
and  the  qua  si -unsteady  "zero"  mass -density 
ratio  as  a  function  of  the  location  of  the  axis 
of  zero  static  coupling.  The  plots  are  for 
various  locations  of  the  center  of  gravity  and 
for  various  frequency  ratios.  In  all  cases, 
the  non-dimensional  radius  of  gyration  about 
the  c  g.  was  considered  to  be  .  50. 
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Figures  18a-18h  show  the  relative  sensi¬ 
tivity  of  various  flutter  boundaries  to 
changes  in  the  non-dimensional  system 
parameters.  These  plots  are  based  upon 


quasi -unsteady  theory. 
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CANTI  LBVtuKBD  FOI  I — , 


More:  TT*  =  [^Vfc+3^]  T1/ r^-*- 2_ 


FIGURE  4.1 
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-THE  inverse  of  the  FLUTTER  FF?CO*U^KOV 

COEFFICIENT  AS  A  FUNCTION  OF  THE  SWEEP 
PARAMETER  FOR  A  CONFIGURATION  FD^ESJlNG 

A  LOW  /C  SOUNDARV. 

PARAMETERS—  rj  —  .Z5 

Xm  -  .X 

nr  z  -  .0411 
a.  =  — .  z_ 


FIGURE  AZ 


o 9'  os’  ov  or  ©z'  ©v  «,cf  ©q-  «9o'  so'  -be 


A  PLOT  OP  THE  IN- FLUID  FLUTTER.  SPEED 
COEFFICIENT  FoR.  THE  MINIMUM  5 WEEP 
PARAMETER  VALUE  FOE  WHICH  A  VALID  Z.ERO 
MASS-DENSITY  Solution  may  e> E  found. 


Note.: 

FLUTTER  5 PEEP  COEFFlClg 
APPROACHES  INFINITY  FOE 
ELASTIC  AXIS  AT  QUARTER 
CHORD.  ,  ii,, 


H 


a 


7r:*=  -is. 


FIGURE  A-A 


FLUTTER  speed  coefficient  &S  A  FUNCTION 
OF  THE  LOCATION  OF  THE  ELASTIC  AXIS. 

MODEL- MASS  AUD  STIFFNESS  OF  BAIRD #"2. 
PARAMETERS-  I —  .IG» 


Xm3  .1&C« 
TTl  “  oil 


-  Ill  - 

A  PLOT  5HOWIN^  "THE-  CRITICAL-  SWEEP  PARAMETER. 
VALUE  FOR  VARIOUS  AND  3l 

Not*.  : 

Tiu9  PLOT  19  BKStD  UPON  “THE.  Q)O^DCA>TlC. 

►.P'^irox. INACTION  PISC0556D  IN  T«l.  lew.  HENC* 
tub,  VAIIDITV  DS«B*sSaS  \£  c^ir  mCRSMCS . 

see  F’ic^ovze.  f-ote.  <nu**te*nc.  tc«.9vl-ts. 


-  l.o  -fe  -.<&  -.A  -•  Z  o  .2 


a 


A>  .»  lo  13.  1,4 
Fi^UWE.  4<2> 
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A  PLOT  or  THE  FLUTTER  SPEED  AMD  FLUTTER  FREQUENCY 
A*  A  FUNCTION  OF  THE  AN^LE  OF  SWEEP  SAC X  POE 
^QUIRES  MODEL.  . 

PARAME.TE.RS—  «  .  25  1=  "08 

Xrn  —  .1^4  -C.  ■*  5 Co  IK. 


U).f  =  CPS 


Jl-  PE^WEtS 


FIGURE  4"l 


